Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,711,755 C→A 100% D45Y (GAC→TAC)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,711,7550CA98.1% 652.7 / 2.2 204D45Y (GAC→TAC) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (96/104);  total (100/104)
Fisher's exact test for biased strand distribution p-value = 5.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 204 total aligned reads displayed.

CTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCAC  >  CP009273/2711610‑2711899
                                                                                                                                                 |                                                                                                                                                
 ttttCCCCCTGAGCCCCCGGTGAGCGGCGTTTCCCGGCGACTCAACACCCGGCGCGCGGGAACCGGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGc                                                                                                                                               <  2:3286899/148‑1 (MQ=255)
      ccACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGGGGGAAACGGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACgag                                                                                                                                          <  2:3813365/148‑1 (MQ=255)
      ccACTCAGGCCACCGTGTTCGGCTTTGCCCGGGGTCACACCCCCCGGCGCGCGGAAACCCGTCCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACgag                                                                                                                                          <  2:2561484/148‑1 (MQ=255)
         cTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGAtt                                                                                                                                       >  1:3083284/1‑148 (MQ=255)
          tGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGAttt                                                                                                                                      >  1:3726193/1‑148 (MQ=255)
          tGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGAttt                                                                                                                                      >  2:2537679/1‑148 (MQ=255)
             tGCCACGGTGAGCGGCCTTTCCCGGCGACACACCCCCCGGCGCGGGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTcc                                                                                                                                   <  1:2826609/148‑1 (MQ=255)
             gccccccGTGAGCGGCCTTGCCCGGCGTCACATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTcc                                                                                                                                   <  2:2478450/143‑1 (MQ=255)
              gCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACACTGTTTGTCATCGCGGGCGGGGCGACGGACACCGTAAGCACGAGATTTTCCt                                                                                                                                  >  2:160964/1‑148 (MQ=255)
              gCCACCGTGACCGGCCTTTCCCGGCGTCTCATCACCCGGCGCGCGGGAACCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCt                                                                                                                                  <  2:3158175/148‑1 (MQ=255)
                cACCGTGATGGGCCTTGGCCGGCGTCACATCCCCGGGCGCGCGGTAACCCGGGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTgg                                                                                                                                <  1:3317921/148‑1 (MQ=255)
                cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTgg                                                                                                                                >  2:1925733/1‑148 (MQ=255)
                cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTgg                                                                                                                                >  1:1830952/1‑148 (MQ=255)
                cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTgg                                                                                                                                >  1:2411912/1‑148 (MQ=255)
                 aCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGa                                                                                                                               >  2:356767/1‑148 (MQ=255)
                    gTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgc                                                                                                                            >  1:26272/1‑148 (MQ=255)
                     tGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgcc                                                                                                                           >  2:1790326/1‑148 (MQ=255)
                      gATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgccc                                                                                                                          >  1:1030633/1‑148 (MQ=255)
                       aTCGGCCTTTGCCGGCTTCTCACCACCCGGCGCGCGGTAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGAcgcccg                                                                                                                         <  2:2784726/148‑1 (MQ=255)
                           gCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGcccc                                                                                                                     >  1:1061105/1‑148 (MQ=255)
                            ccttggccGGCGTCACAACACCCGGCGCGGGGGAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGccccc                                                                                                                    <  2:2614940/142‑1 (MQ=255)
                             ttgggccGGCGCCACACCCCCGGGCGCGCGGCACCCGGGCCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCt                                                                                                                   <  1:162289/143‑1 (MQ=255)
                             cTTTGCCGGCGTCTCACCCCCGGGGGCGCGTCAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCt                                                                                                                   <  1:2091943/148‑1 (MQ=255)
                                     gCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGt                                                                                                           >  2:1434756/1‑148 (MQ=255)
                                           catcaCCCGGCGCGCGGTAAACGGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCgg                                                                                                     <  2:1161509/148‑1 (MQ=255)
                                             tcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGt                                                                                                   >  1:14864/1‑148 (MQ=255)
                                             tcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGt                                                                                                   >  1:2156721/1‑148 (MQ=255)
                                              caCCCGGCGCGCGGCAACCGGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGtt                                                                                                  <  2:446704/148‑1 (MQ=255)
                                                ccgggggcgcgGGACCCCGGCCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtgt                                                                                                <  1:732391/142‑1 (MQ=255)
                                                 ccGGCGCGCGGTAAACCGGGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTtgtg                                                                                               <  1:935404/148‑1 (MQ=255)
                                                    gcgcgcGGGAACCCGGCCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGt                                                                                            <  2:1830898/148‑1 (MQ=255)
                                                       cgcgGGAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACg                                                                                         <  2:2135489/148‑1 (MQ=255)
                                                        gcgGCAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGt                                                                                        <  1:2132267/148‑1 (MQ=255)
                                                             aaaCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATc                                                                                   >  1:647901/1‑148 (MQ=255)
                                                               aCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAg                                                                                 >  2:269431/1‑148 (MQ=255)
                                                               aCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAg                                                                                 >  2:310278/1‑148 (MQ=255)
                                                               aCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAg                                                                                 >  2:733914/1‑148 (MQ=255)
                                                                cggtcccccACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                                                                <  1:1110910/141‑1 (MQ=255)
                                                                ccGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                                                                >  2:1550854/1‑148 (MQ=255)
                                                                ccGTGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                                                                <  2:1619351/148‑1 (MQ=255)
                                                                ccGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                                                                <  1:3576687/148‑1 (MQ=255)
                                                                 cGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGaa                                                                               >  2:2871863/1‑148 (MQ=255)
                                                                   cgccccACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAatc                                                                             <  1:18229/144‑1 (MQ=255)
                                                                          cGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTGTGTCAGCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAAtcatcatca                                                                      <  2:3230103/148‑1 (MQ=255)
                                                                          cGGCGAATCGTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAAtcatcatca                                                                      <  2:144321/148‑1 (MQ=255)
                                                                           ggCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCAc                                                                     <  1:2008054/148‑1 (MQ=255)
                                                                             cGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTg                                                                   >  1:3803178/1‑148 (MQ=255)
                                                                               aaTCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCg                                                                 >  1:2936006/1‑148 (MQ=255)
                                                                               aaTCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCg                                                                 >  1:3102197/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  2:515344/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  1:1752814/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  2:65823/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  1:3533118/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  1:600181/1‑148 (MQ=255)
                                                                                    tCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacat                                                            >  1:2993725/1‑148 (MQ=255)
                                                                                          ttGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCa                                                      <  2:188677/148‑1 (MQ=255)
                                                                                           tGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCAc                                                     <  2:1571238/148‑1 (MQ=255)
                                                                                            gCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACt                                                    >  1:3380726/1‑148 (MQ=255)
                                                                                            gCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACt                                                    <  2:1340072/148‑1 (MQ=255)
                                                                                            gCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACt                                                    <  1:1002941/148‑1 (MQ=255)
                                                                                               gTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTg                                                 >  2:3208683/1‑148 (MQ=255)
                                                                                                  aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                              <  2:1375859/148‑1 (MQ=255)
                                                                                                  aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                              >  2:3121261/1‑148 (MQ=255)
                                                                                                  aCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGAcc                                              <  1:1312461/148‑1 (MQ=255)
                                                                                                     gggCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTg                                           <  1:1283160/148‑1 (MQ=255)
                                                                                                      ggCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTCCGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGc                                          >  1:1167042/1‑148 (MQ=255)
                                                                                                          gACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCt                                      <  1:3042843/148‑1 (MQ=255)
                                                                                                          gACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCt                                      <  1:2729901/148‑1 (MQ=255)
                                                                                                           aCTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtt                                     >  2:2944154/1‑148 (MQ=255)
                                                                                                           aCTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtt                                     <  1:1109730/148‑1 (MQ=255)
                                                                                                           aCTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtt                                     <  1:526023/148‑1 (MQ=255)
                                                                                                            ctctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCttt                                    >  1:1814886/1‑148 (MQ=255)
                                                                                                             tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                   <  2:55570/148‑1 (MQ=255)
                                                                                                             tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                   >  1:3623066/1‑148 (MQ=255)
                                                                                                             tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                   <  1:1790326/148‑1 (MQ=255)
                                                                                                             tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                   <  1:3683063/148‑1 (MQ=255)
                                                                                                             tctGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCtttt                                   <  2:2446517/148‑1 (MQ=255)
                                                                                                                gTTTGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTAcc                                >  2:1983678/1‑148 (MQ=255)
                                                                                                                 tttGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCt                               <  2:1061105/148‑1 (MQ=255)
                                                                                                                  ttGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtt                              >  2:1486500/1‑148 (MQ=255)
                                                                                                                   tGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCttt                             <  2:1781136/148‑1 (MQ=255)
                                                                                                                   tGTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCttt                             <  2:3364864/148‑1 (MQ=255)
                                                                                                                    gTCATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttc                            >  1:1154002/1‑148 (MQ=255)
                                                                                                                      cATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcat                          >  1:1891871/1‑148 (MQ=255)
                                                                                                                      cATCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcat                          <  2:1254491/148‑1 (MQ=255)
                                                                                                                        tCGCGGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcattt                        >  2:539296/1‑148 (MQ=255)
                                                                                                                         cgcgGGCGGGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTAACtttcatttc                       >  2:402882/1‑148 (MQ=255)
                                                                                                                          gcgggcggGGCGACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCtttcatttca                      <  1:1369198/148‑1 (MQ=255)
                                                                                                                              gcggggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGt                  <  1:1671902/148‑1 (MQ=255)
                                                                                                                                 gggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCa               <  1:2834154/148‑1 (MQ=255)
                                                                                                                                 gggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCa               <  2:3069621/148‑1 (MQ=255)
                                                                                                                                 gggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCa               >  2:3012934/1‑148 (MQ=255)
                                                                                                                                  ggcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCAt              >  2:144780/1‑148 (MQ=255)
                                                                                                                                   gcgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATg             <  1:3714558/148‑1 (MQ=255)
                                                                                                                                    cgACGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGc            <  2:3726122/148‑1 (MQ=255)
                                                                                                                                      actgCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtt          <  1:160964/145‑1 (MQ=255)
                                                                                                                                      aCGGCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtt          <  2:17019/148‑1 (MQ=255)
                                                                                                                                        ggCCACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCtttt        <  1:3692599/148‑1 (MQ=255)
                                                                                                                                          ccACCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTGTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTtct      <  2:1154207/148‑1 (MQ=255)
                                                                                                                                            aCCGTAAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTcc    >  2:3502122/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTGACCTTGCCGCTTTTACCTTTCATTTCATCGTTCATGCTTTTCTCCAC  >  CP009273/2711610‑2711899

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