Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 694,721 C→T G79S (GGC→AGC)  miaB ← isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913694,7210CT98.7% 248.6 / NA 75G79S (GGC→AGC) miaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/43);  minor base G (0/1);  total (31/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGC  >  NC_000913/694474‑694963
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
cAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGAt                                                                                                                                                                                                                                                 <  2:131656/251‑1 (MQ=255)
  atatTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGCCGCTGATa                                                                                                                                                                                                                                                >  2:472034/1‑250 (MQ=255)
  atatTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATa                                                                                                                                                                                                                                                >  1:470141/1‑250 (MQ=255)
      ttattGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCTCCGACGCTGATAATCa                                                                                                                                                                                                                                            <  1:454881/250‑1 (MQ=255)
        attGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCACTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGCAGCTGATATCTACCACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGt                                                                                                                                                                                                                                         >  1:313168/1‑251 (MQ=255)
            cAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTgg                                                                                                                                                                                                                                     <  1:335302/251‑1 (MQ=255)
            cAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTgg                                                                                                                                                                                                                                     >  2:416198/1‑251 (MQ=255)
                  tCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAAGCAGGTCTGGAttct                                                                                                                                                                                                                                <  1:107647/250‑1 (MQ=255)
                   ccATGATGGAGACAAACGCGGTCGGCCCTGCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCGCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAttcttc                                                                                                                                                                                                                              <  1:551758/251‑1 (MQ=255)
                     atgaagGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTATACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGTCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTtctc                                                                                                                                                                                                                            <  2:341530/251‑1 (MQ=255)
                       gatgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCtt                                                                                                                                                                                                                          <  1:23664/251‑1 (MQ=255)
                         tgGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGCTGCCACGCAGCCACCGACGCTGATCATCCGGGCTGGATTCTTCTCtttt                                                                                                                                                                                                                        >  2:170304/1‑251 (MQ=255)
                           gagaAAAACGGGGTCGGCCAGTCGGCGAGCGGTTCCGGTAGACGGTCAAACTGCGCGATTTCCGGGAAGCTGATATCTACAACCGGGATGCGGTCGCCACGAACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCAAAAAATAATATCGACAGAGTGGGCGCGCTGGCGAATGT‑GCTCGCATTAT‑TGCGATGCCACGCAGCCACCGGCGCTGATAATCAGGTCTGGATTCTTCTCTTTTaa                                                                                                                                                                                                                      <  1:529559/251‑1 (MQ=255)
                            agaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCCTCTCTTTTAAc                                                                                                                                                                                                                     >  2:359327/1‑251 (MQ=255)
                            agaCAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGACGCCACGCACGGAGTTGATCATCTCCGGCAGACGTTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTAGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATACTCAGGTCTGGATTCTTCTCTTTTAAc                                                                                                                                                                                                                     >  1:235047/1‑251 (MQ=255)
                                 aaCGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAGGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGttt                                                                                                                                                                                                                >  1:455126/1‑251 (MQ=255)
                                   cgcgGTCGGCCCTTCGGCGCGCGGTGCCGGTAGACGGTCAAACTTCTCGATTGCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTcc                                                                                                                                                                                                              <  1:28069/251‑1 (MQ=255)
                                         cGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAc                                                                                                                                                                                                        <  1:280299/251‑1 (MQ=255)
                                          ggCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGACACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGAcc                                                                                                                                                                                                       >  2:551045/1‑251 (MQ=255)
                                           gCCCTGCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATGTCCGGGAAGCTGATATCTACAACCGGGCTGCTGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCGCTTTTAACAGTTTCCAGCGAccc                                                                                                                                                                                                      <  2:29337/251‑1 (MQ=255)
                                              cTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGCCATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAAc                                                                                                                                                                                                   >  2:217032/1‑251 (MQ=255)
                                                tCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTg                                                                                                                                                                                                 >  2:605641/1‑251 (MQ=255)
                                                    cgcgcgGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATgg                                                                                                                                                                                             >  2:526268/1‑251 (MQ=255)
                                                    cgcgcgGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATgg                                                                                                                                                                                             <  2:470141/251‑1 (MQ=255)
                                                         ggTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGAc                                                                                                                                                                                        >  1:238145/1‑251 (MQ=255)
                                                            tCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACtt                                                                                                                                                                                      <  1:342570/250‑1 (MQ=255)
                                                                     cGGTCACACTTCTCGATTTCCGGGAAGCTGATATATACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGGTTGCGGCCCAACAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCCCGCCTTCT‑TGCGATGCCACGCAGCCACCGCCGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGc                                                                                                                                                                            >  2:571121/1‑251 (MQ=255)
                                                                       gTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCt                                                                                                                                                                          <  1:449065/251‑1 (MQ=255)
                                                                           aaCTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGCCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTtctc                                                                                                                                                                      >  1:518570/1‑251 (MQ=255)
                                                                                    atgtCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGa                                                                                                                                                             <  2:86323/248‑1 (MQ=255)
                                                                                        ccGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCCCGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCCTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAg                                                                                                                                                         <  1:36590/251‑1 (MQ=255)
                                                                                         cGGGAAGCTGATATCTACAACCGGTCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAgg                                                                                                                                                        <  1:155778/251‑1 (MQ=255)
                                                                                             aaGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCCGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGCCGCTGATAATCAGGTCTGGCTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGtt                                                                                                                                                    >  2:553044/1‑251 (MQ=255)
                                                                                             aaGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGtt                                                                                                                                                    <  1:448315/251‑1 (MQ=255)
                                                                                              aGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTc                                                                                                                                                   <  2:135416/251‑1 (MQ=255)
                                                                                                   atatCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGCGCAGGTGTTcagc                                                                                                                                                >  2:137731/1‑249 (MQ=255)
                                                                                                         aCAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCagcagcagc                                                                                                                                          <  1:575683/249‑1 (MQ=255)
                                                                                                             ccGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATcc                                                                                                                                    >  2:174814/1‑251 (MQ=255)
                                                                                                                     cGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGGGCAGCGTCTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGAGTCTTCGCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCACCGCATCCGCttcttc                                                                                                                            <  2:101059/251‑1 (MQ=255)
                                                                                                                         cGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGc                                                                                                                         <  2:455126/250‑1 (MQ=255)
                                                                                                                                       ttGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTg                                                                                                          >  1:156214/1‑251 (MQ=255)
                                                                                                                                                tcCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCTGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtg                                                                                                 >  1:568847/1‑251 (MQ=255)
                                                                                                                                                 ccGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgt                                                                                                  <  1:295433/249‑1 (MQ=255)
                                                                                                                                                 ccGGAAGACGGTGCAGCGGTGGCGGCACAAAAATAGGAACGACAAAGTGGGCGCGATGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCCCGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGTAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTAAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgg                                                                                                <  2:432648/251‑1 (MQ=255)
                                                                                                                                                   ggCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGGCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCgtgggg                                                                                              >  2:449389/1‑250 (MQ=255)
                                                                                                                                                      aGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGc                                                                                           >  1:437024/1‑251 (MQ=255)
                                                                                                                                                        gcGGTGCAGCGTTGGCGGCCCAAAAAAAATATCGACTGAGGGGGCGCGCGGGCGAA‑GTCGCTCGCCT‑CTGTGCGATGCCAGGCAGCCACCGACGCGGATAATCAGGTCTGGATCCTGCTCGTCTAACAGGTCCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGCTCAGCAGCAGCGCATCCGTTTCTTCCGCCACGTCGGTCATTTGATTGCCGTGGGCGGCAt                                                                                         <  2:483224/250‑1 (MQ=255)
                                                                                                                                                         gggTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGGGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATc                                                                                        <  2:192660/250‑1 (MQ=255)
                                                                                                                                                             gCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCcagc                                                                                    <  1:108378/251‑1 (MQ=255)
                                                                                                                                                                 cGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCcagc                                                                                    >  1:35849/1‑247 (MQ=255)
                                                                                                                                                                  gTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATc                                                                               <  2:537853/251‑1 (MQ=255)
                                                                                                                                                                  gTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATc                                                                               <  1:15832/251‑1 (MQ=255)
                                                                                                                                                                        ggCCCAAAAATAATATCTACATAGTGGGCGCTCTGGCGAATGT‑GCTCACCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGTTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGGTTGAGCAGGTGTTCAGCAGCATCACATCCGCTTCTTCCGCCACGTCGGTCAGTTTGTAGCCgtgtgaggcaccgcgccggacggccat                                                                         >  2:326021/1‑227 (MQ=255)
                                                                                                                                                                             aaaaaTAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCg                                                                    <  1:472034/251‑1 (MQ=255)
                                                                                                                                                                                aataatATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATg                                                                 <  1:170304/251‑1 (MQ=255)
                                                                                                                                                                                   aataTCGACATGGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCCGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTGGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAGCTTCGATGAAt                                                              <  2:316304/251‑1 (MQ=255)
                                                                                                                                                                                   aataTCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGCTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAAt                                                              >  1:71777/1‑251 (MQ=255)
                                                                                                                                                                                        cGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATTGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGTTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTAc                                                         <  1:551045/251‑1 (MQ=255)
                                                                                                                                                                                         gACATAGTGGGCGCGCTGGCGAAGGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACt                                                        <  2:253330/251‑1 (MQ=255)
                                                                                                                                                                                           cATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGCTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTATCGCGGATTGAGCAGGTGTTCAGAAGCAGCACATCCGCTTCTTCTGCAACGTCGGTCAGTTGATATCCGTGGGTGGCATCTAGCGGATCGGCCATCTTCGCTGAATCGTACTCg                                                      >  2:337050/1‑251 (MQ=255)
                                                                                                                                                                                              aGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGtt                                                    <  1:605641/250‑1 (MQ=255)
                                                                                                                                                                                               gTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTc                                                   >  1:454519/1‑250 (MQ=255)
                                                                                                                                                                                                           gCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCCGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGCTGAATCGTACTCGTTCATCTGACAGcccc                                      >  2:482807/1‑251 (MQ=255)
                                                                                                                                                                                                           gCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGcccc                                      <  1:577576/251‑1 (MQ=255)
                                                                                                                                                                                                            cGAATGT‑GCTCCCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCATCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCa                                     <  1:279633/251‑1 (MQ=255)
                                                                                                                                                                                                                    gCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCAGCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGt                                  <  1:105772/247‑1 (MQ=255)
                                                                                                                                                                                                                    gCTCGCCTTCT‑TGCGATGCAACGCAGCCACCGAAGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAggc                                  >  2:105772/1‑246 (MQ=255)
                                                                                                                                                                                                                                tGCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAg                      <  1:174814/248‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGTCCCAACTGATGGAAGACTTTTGCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGGGGGTGGCATCCAGCAGATCGTCCATCTTCGATGAATCGTACTCGTTCATCTGACATCCCCAg                                    <  1:408764/233‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGATGACACGCAGCCACCGACGCTGATACTCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGCTTGCGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGAGGGTGGACTCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAg                                    >  2:408764/1‑233 (MQ=255)
                                                                                                                                                                                                                                   gATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGATCAGGGGTTCAGCAGCAGCACAGCCGCGTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGGTCAGCTGACAGCGCCAGGTTTTAATATGGAGtttttt                <  2:53825/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gATGCCACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCGTCTTCCGCCACGTCGGTCAGTGGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGtttttt                <  1:498486/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATTGAGTTTTTTGGTc            <  1:482807/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ccACGCAGCCACCGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTc            >  1:494228/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cccGACGCTGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCCGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTTTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTGCTCGTTCATCTGACTGCCCCAGTTTTTAATATGGAGTTTGTTCGTCATCGACTTGc  <  2:348962/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
CAATATTTATTGCAGCCTTCCATGATGGAGACAAACGCGGTCGGCCCTTCGGCGCGCGGTTCCGGTAGACGGTCAAACTTCTCGATTTCCGGGAAGCTGATATCTACAACCGGGCTGCGGTCGCCACGCACGGAGTTGATCATCTCCGGCAGACGGTGCAGCGTTTGCGGCCCAAAAATAATATCGACATAGTGGGCGCGCTGGCGAATGT‑GCTCGCCTTCT‑TGCGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGC  >  NC_000913/694474‑694963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑