Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,308 A→T 7.8% I289N (ATT→AAT) ‡ ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,3080AT 7.8% 407.7 / 18.2 116I289N (ATT→AAT) ‡ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (50/57);  new base T (4/5);  total (54/62)
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.77e-01

Only 100 of 123 total aligned reads displayed.

GTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCA  >  NC_000913/3066160‑3066452
                                                                                                                                                     |                                                                                                                                                 
     gggATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGc                                                                                                                                             <  2:374434/149‑1 (MQ=255)
         ttATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                         <  3:698599/149‑1 (MQ=255)
          tATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGaa                                                                                                                                         <  4:1071922/149‑1 (MQ=255)
            tCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                                                                                                      >  2:521113/1‑149 (MQ=255)
              tGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                    <  2:113689/149‑1 (MQ=255)
              tGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                    <  2:1063772/149‑1 (MQ=255)
                cATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGt                                                                                                                                  <  1:816320/149‑1 (MQ=255)
                  tCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                    <  2:1403307/145‑1 (MQ=255)
                      tATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAga                                                                                                                                       >  2:248945/1‑138 (MQ=255)
                      tATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAga                                                                                                                                       <  1:248945/139‑2 (MQ=255)
                        tCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                                            <  1:757580/131‑1 (MQ=255)
                           ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                                                                                                       <  3:1089963/149‑1 (MQ=255)
                                  gcgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGAtt                                                                                                                >  2:467931/1‑149 (MQ=255)
                                    gcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                              >  4:1338185/1‑149 (MQ=255)
                                         gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         <  2:1082802/149‑1 (MQ=255)
                                             aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                         >  4:303393/1‑113 (MQ=255)
                                             aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCACGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                                                                                                         <  3:303393/113‑1 (MQ=255)
                                               atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                                                             >  4:198081/1‑123 (MQ=37)
                                               atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                                                             <  3:198081/123‑1 (MQ=25)
                                                tatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                  <  4:423930/149‑1 (MQ=255)
                                                 atTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAg                                                                                                 >  2:1136479/1‑149 (MQ=255)
                                                     accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                    <  2:1231409/110‑1 (MQ=255)
                                                     accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATc                                                                                                                                    >  1:1231409/1‑110 (MQ=255)
                                                     accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                             >  1:488675/1‑133 (MQ=255)
                                                     accacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                             <  2:488675/133‑1 (MQ=255)
                                                     accacGATCTTAGCCACGTCCTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                             >  1:1161424/1‑149 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                       <  3:1060710/105‑1 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                                                                                                       >  4:1060710/1‑105 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACga                                                                                                                     <  3:518989/123‑1 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                           <  2:381249/149‑1 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTc                                                                                           >  1:381249/1‑149 (MQ=255)
                                                       cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAgat                                                                                                                                      >  2:619011/1‑105 (MQ=37)
                                                         cGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                                                         <  3:1338185/149‑1 (MQ=255)
                                                           aTCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                                     >  1:25617/1‑103 (MQ=255)
                                                           aTCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCAt                                                                                                                                     <  2:25617/103‑1 (MQ=255)
                                                            tCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                                      <  2:1161424/149‑1 (MQ=255)
                                                            tCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                                      <  1:467931/149‑1 (MQ=255)
                                                            tCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCa                                                                                      >  1:542740/1‑149 (MQ=255)
                                                                   cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                              >  2:517663/1‑118 (MQ=255)
                                                                   cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                              <  1:517663/118‑1 (MQ=255)
                                                                       tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGc                                                                           >  3:1399098/1‑149 (MQ=255)
                                                                         ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                             >  3:710809/1‑113 (MQ=255)
                                                                         ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                             <  4:710809/113‑1 (MQ=255)
                                                                         ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                     <  1:1072022/121‑1 (MQ=255)
                                                                         ttCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                     >  2:1072022/1‑121 (MQ=255)
                                                                             aCTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                                                     <  2:542740/149‑1 (MQ=255)
                                                                               ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                                            >  4:1459451/1‑76 (MQ=255)
                                                                               ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCt                                                                                                                                            <  3:1459451/76‑1 (MQ=255)
                                                                               ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCATATTGCTGAAga                                                                                                                                       <  2:875637/81‑2 (MQ=255)
                                                                               ttACACCGTAATTTAATTGCAACCGGCTCGATGCTGGATTAATGACTCGTCCCGTCACGATCTTCAGCATATTGCTGAAga                                                                                                                                       >  1:875637/1‑80 (MQ=255)
                                                                                   acCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                                                       >  1:476398/1‑125 (MQ=255)
                                                                                     cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATAATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                                             >  4:641592/1‑149 (MQ=255)
                                                                                     cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATAATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                                             <  3:641592/149‑1 (MQ=255)
                                                                                      gTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                                                                                                              >  2:828088/1‑99 (MQ=255)
                                                                                       taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCg                                                                                                                                   <  2:660467/77‑1 (MQ=255)
                                                                                       taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCg                                                                                                                                   >  1:660467/1‑77 (MQ=255)
                                                                                       taatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                                           <  1:1136479/149‑1 (MQ=255)
                                                                                          tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                                                         >  2:715400/1‑116 (MQ=255)
                                                                                          tttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                                                         <  1:715400/116‑1 (MQ=255)
                                                                                            taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAatat                                                      <  3:630323/149‑1 (MQ=255)
                                                                                                tGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCa                                                  <  4:1399098/149‑1 (MQ=255)
                                                                                                   aaCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                               >  3:1388526/1‑149 (MQ=255)
                                                                                                          tCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                      <  1:535862/87‑1 (MQ=255)
                                                                                                          tCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                      >  2:535862/1‑87 (MQ=255)
                                                                                                          tCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                                                       >  4:317722/1‑102 (MQ=255)
                                                                                                          tCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc                                                                                       <  3:317722/102‑1 (MQ=255)
                                                                                                            gATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                                                             <  4:297283/110‑1 (MQ=255)
                                                                                                            gATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCg                                                                             >  3:297283/1‑110 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAata                                                                                                                                                 <  2:1221337/41‑1 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAata                                                                                                                                                 >  1:1221337/1‑41 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         >  3:1121514/1‑80 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         <  4:1121514/80‑1 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         >  1:1010567/1‑80 (MQ=255)
                                                                                                              tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         <  2:1010567/80‑1 (MQ=255)
                                                                                                                  gggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCa                                <  4:205237/149‑1 (MQ=255)
                                                                                                                   ggTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAg                               >  2:39445/1‑149 (MQ=255)
                                                                                                                        aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                  >  3:116340/1‑77 (MQ=255)
                                                                                                                        aTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                  <  4:116340/77‑1 (MQ=255)
                                                                                                                         tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                              >  1:1202886/1‑96 (MQ=255)
                                                                                                                         tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATc                                                                              <  2:1202886/96‑1 (MQ=255)
                                                                                                                         tGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCa                                   <  3:278660/139‑1 (MQ=255)
                                                                                                                             tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                                                     <  2:1347737/101‑1 (MQ=255)
                                                                                                                             tCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGc                                                                     >  1:1347737/1‑101 (MQ=255)
                                                                                                                                    gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg              <  4:1388526/149‑1 (MQ=255)
                                                                                                                                       aCGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGa           >  4:1300528/1‑149 (MQ=255)
                                                                                                                                         gatgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                                           <  1:577500/95‑1 (MQ=37)
                                                                                                                                         gatgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaa                                                           >  2:577500/5‑99 (MQ=38)
                                                                                                                                         gatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                                             <  3:79465/93‑1 (MQ=25)
                                                                                                                                         gatgcTTCAGCATATTGCTGAAGATCGTGACGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATc                                                             >  4:79465/5‑97 (MQ=25)
                                                                                                                                         gATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATc         >  2:803688/1‑149 (MQ=255)
                                                                                                                                         gATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATc         <  1:803688/149‑1 (MQ=255)
                                                                                                                                          aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCa        >  2:212496/1‑149 (MQ=255)
                                                                                                                                           tgcTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg             >  4:713063/3‑143 (MQ=255)
                                                                                                                                              ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCa                                   <  1:1073548/119‑1 (MQ=255)
                                                                                                                                              ttCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg     >  2:803168/1‑149 (MQ=255)
                                                                                                                                               tCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGt    >  1:272872/1‑149 (MQ=255)
                                                                                                                                                 aGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCa  <  2:272872/149‑1 (MQ=255)
                                                                                                                                                  gCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                         >  1:573464/1‑45 (MQ=255)
                                                                                                                                                  gCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg             >  3:825279/1‑137 (MQ=255)
                                                                                                                                                    aaTATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                             <  1:152199/55‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCA  >  NC_000913/3066160‑3066452

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