Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,897 C→T 5.8% M13I (ATG→ATA) ‡ yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8970CT 5.8% 409.1 / 10.8 121M13I (ATG→ATA) ‡yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (54/60);  new base T (3/4);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 121 total aligned reads displayed.

TTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATA  >  NC_000913/383752‑384044
                                                                                                                                                 |                                                                                                                                                   
  aCCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                                                                >  1:135447/1‑149 (MQ=255)
      ggCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATg                                                                                                                                            <  2:200898/149‑1 (MQ=255)
        cAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                                                                          <  3:1449555/149‑1 (MQ=255)
            aaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATa                                                                                                                                      <  4:343772/149‑1 (MQ=255)
                  ggTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATgg                                                                                                                                           >  3:1268853/1‑138 (MQ=255)
                  ggTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATgg                                                                                                                                           <  4:1268853/138‑1 (MQ=255)
                      aCGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGGCTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                                                                            >  4:1189145/1‑133 (MQ=255)
                      aCGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                                                                            <  3:1189145/133‑1 (MQ=255)
                        gCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATc                                                                                                                          <  1:263671/149‑1 (MQ=255)
                           gTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCt                                                                                                                                               >  2:481885/1‑125 (MQ=255)
                            tAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCa                                                                                                                        >  4:32237/1‑147 (MQ=255)
                            tAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCa                                                                                                                        <  3:32237/147‑1 (MQ=255)
                            tAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAgg                                                                                                                      >  4:1049674/1‑149 (MQ=255)
                              gATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                                                                          >  1:627919/1‑127 (MQ=255)
                                  aTCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                                                                <  4:513261/149‑1 (MQ=255)
                                  aTCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCa                                                                                                                <  2:135447/149‑1 (MQ=255)
                                       gTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCaa                                                                                                           >  3:357493/1‑149 (MQ=255)
                                                ccAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaa                                                                                                  <  4:357493/149‑1 (MQ=255)
                                                 cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTc                                                                                                                 >  3:903310/1‑133 (MQ=255)
                                                 cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGAATGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCATGTTTTCATGCAATTTGCTCaaa                                                                                                 >  1:1253557/1‑149 (MQ=255)
                                                  aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                                                                >  1:820462/1‑149 (MQ=255)
                                                      acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCt                                                                                                                                               <  1:1340464/98‑1 (MQ=255)
                                                      acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCt                                                                                                                                               >  2:1340464/1‑98 (MQ=255)
                                                      acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATc                                                                                            >  1:1350651/1‑149 (MQ=255)
                                                        gggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                                                                          >  4:459561/1‑149 (MQ=255)
                                                        gggCCTTCGTTATACGCAGGTATAATGGCGTCAACACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAt                                                                                          >  1:534423/1‑149 (MQ=255)
                                                         ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtt                                                                                         <  2:1350651/149‑1 (MQ=255)
                                                         ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtt                                                                                         <  2:820462/149‑1 (MQ=255)
                                                          gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                                                                        <  4:721023/149‑1 (MQ=255)
                                                           ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                       >  2:1374839/1‑149 (MQ=255)
                                                             ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                                                                         >  3:382544/1‑97 (MQ=255)
                                                                gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGaa                                                                                  <  1:1421992/149‑1 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                                                                  <  4:576585/84‑1 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACAtt                                                                                                                                                  >  3:576585/1‑84 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                       >  1:535258/1‑111 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                       <  2:535258/111‑1 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                                                                   <  3:1148757/115‑1 (MQ=255)
                                                                 ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                                                                   >  4:1148757/1‑115 (MQ=255)
                                                                       gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                  <  3:1444972/94‑1 (MQ=255)
                                                                       gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGa                                                                                                                                  >  4:1444972/1‑94 (MQ=255)
                                                                       gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                                                                           >  2:1108763/1‑149 (MQ=255)
                                                                       gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                                                                           <  1:1374839/149‑1 (MQ=255)
                                                                                aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                                                                >  4:282830/1‑119 (MQ=255)
                                                                                aTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                                                                <  3:282830/119‑1 (MQ=255)
                                                                                   gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                 <  3:629813/115‑1 (MQ=255)
                                                                                   gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                 >  4:629813/1‑115 (MQ=255)
                                                                                     gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                       >  1:1060895/1‑123 (MQ=255)
                                                                                     gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                                                                       >  2:1001801/1‑139 (MQ=255)
                                                                                     gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCt                                                                       <  1:1001801/139‑1 (MQ=255)
                                                                                         aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                       <  3:105678/87‑1 (MQ=255)
                                                                                         aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt                                                         >  2:1220751/1‑149 (MQ=255)
                                                                                           aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                                                                           <  3:1107945/65‑1 (MQ=255)
                                                                                           aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                                                                           >  4:1107945/1‑65 (MQ=255)
                                                                                           aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                       <  2:1172480/117‑1 (MQ=255)
                                                                                           aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                                       >  1:1172480/1‑117 (MQ=255)
                                                                                               ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                       <  2:334942/81‑1 (MQ=255)
                                                                                               ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAg                                                                                                                       >  1:334942/1‑81 (MQ=255)
                                                                                               ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                                                                 <  1:100821/119‑1 (MQ=255)
                                                                                               ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                                                                 >  2:100821/1‑119 (MQ=255)
                                                                                               ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAc                                                   >  2:1303390/1‑149 (MQ=255)
                                                                                                ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTACATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                                                  >  4:498479/1‑149 (MQ=255)
                                                                                                  ttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtg                                                <  2:534423/149‑1 (MQ=255)
                                                                                                     ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgttg                                             >  1:615421/1‑149 (MQ=255)
                                                                                                            tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                                                                   <  3:128739/104‑1 (MQ=255)
                                                                                                            tGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                                                                   >  4:128739/1‑104 (MQ=255)
                                                                                                               gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAt                                   <  1:1220751/149‑1 (MQ=255)
                                                                                                                    cGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGc                              <  2:615421/149‑1 (MQ=255)
                                                                                                                        ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                                                                     <  4:82407/74‑1 (MQ=255)
                                                                                                                        ttAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                                                                     >  3:82407/1‑74 (MQ=255)
                                                                                                                              tACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                 >  2:240196/1‑72 (MQ=255)
                                                                                                                              tACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                                                 <  1:240196/72‑1 (MQ=255)
                                                                                                                               aCTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                                                                        >  3:1004447/1‑80 (MQ=255)
                                                                                                                               aCTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                                                                        <  4:1004447/80‑1 (MQ=255)
                                                                                                                                cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAg                  <  1:1108763/149‑1 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGc                                                                                                      >  1:1378783/1‑64 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                                                                           >  3:468767/1‑91 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACaaa                                                                           <  4:468767/91‑1 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt                                                      >  1:861002/1‑112 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt                                                      <  2:861002/112‑1 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt                                                       >  2:1466161/1‑111 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg                                                     >  4:1145570/1‑113 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg                                                     <  3:1145570/113‑1 (MQ=255)
                                                                                                                                 ttAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAgg                 >  4:1374285/1‑149 (MQ=255)
                                                                                                                                  tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCCTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGc                <  2:1253557/149‑1 (MQ=255)
                                                                                                                                    aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgtt                                              >  1:289425/1‑117 (MQ=255)
                                                                                                                                    aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGc                              <  2:822203/133‑1 (MQ=255)
                                                                                                                                    aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGc                              >  1:822203/1‑133 (MQ=255)
                                                                                                                                     aaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGAt             <  1:1303390/149‑1 (MQ=255)
                                                                                                                                      aaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTcc                                          <  1:1287879/119‑1 (MQ=255)
                                                                                                                                      aaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGAt             >  3:602164/1‑148 (MQ=255)
                                                                                                                                        ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                                                                     <  2:264419/74‑1 (MQ=255)
                                                                                                                                        ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                                                                     >  1:264419/1‑74 (MQ=255)
                                                                                                                                         cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                                                                     <  3:859095/73‑1 (MQ=255)
                                                                                                                                         cATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                                                                     >  4:859095/1‑73 (MQ=255)
                                                                                                                                           tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTACATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGAc       <  3:498479/149‑1 (MQ=255)
                                                                                                                                            aGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGaaa                                    <  4:355592/119‑1 (MQ=255)
                                                                                                                                            aGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGaaa                                    >  3:355592/1‑119 (MQ=255)
                                                                                                                                            aGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATtgt                                <  2:729821/123‑1 (MQ=255)
                                                                                                                                               aaCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATc            >  2:1215021/1‑140 (MQ=255)
                                                                                                                                                aCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCata  <  4:254008/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCACGCAGGCGATCATGACGCATA  >  NC_000913/383752‑384044

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: