Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,263,409 C→T 100% intergenic (+448/+166) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,263,4090CT95.1% 253.0 / 5.5 109intergenic (+448/+166)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (43/59);  total (46/63)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 112 total aligned reads displayed.

TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552
                                                                                                                                                |                                                                                                                                               
tttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                              <  2:2425087/148‑1 (MQ=37)
tttGGGTATAGTGTCGTGGACAGTCAGTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATtttt                                                                                                                                              <  2:2740795/148‑1 (MQ=255)
 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                             <  2:3706332/148‑1 (MQ=255)
 ttGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                             <  2:3171720/148‑1 (MQ=37)
 gtGGGTATAGTGTCGTGGACACTCATTCATCTTACGGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTc                                                                                                                                             <  2:3706348/147‑1 (MQ=255)
       ataGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCGTTTGGGTATAGTGTCGTGGGCAGGGCTTGATTTTTCTGgccc                                                                                                                                       >  1:810820/1‑148 (MQ=11)
         aGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                     <  2:3474747/148‑1 (MQ=37)
         aGTGTCGTGGACACTCATTCATCTTTCTGCCACTCAAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTc                                                                                                                                     <  2:952236/148‑1 (MQ=255)
              cGTGGACAGTGATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCaaaa                                                                                                                                <  1:1737506/148‑1 (MQ=255)
               gttgACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAg                                                                                                                               <  2:2070834/145‑1 (MQ=255)
                       tcattcatCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       <  1:1629292/148‑1 (MQ=255)
                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       >  1:1127106/1‑148 (MQ=37)
                       tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       <  2:2043591/148‑1 (MQ=37)
                          ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTTACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGTCATTAAATTTTCTGCCCCACCCAAAACAAAAAccccc                                                                                                                    >  1:1984573/1‑146 (MQ=11)
                               cTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                               <  1:1853155/148‑1 (MQ=37)
                               cTTTCGGCCCCTCCAAAAGTAAAAACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                                                                                               <  1:1853140/148‑1 (MQ=255)
                                 ttCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTTCGCCTCCAAAATCAAAAACCCCCCGCAGCg                                                                                                             >  1:2604916/1‑148 (MQ=11)
                                 ttCTGCCACTCCAAAAGTAAAAACCCGCCGAATCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                             <  2:2724236/148‑1 (MQ=255)
                                  tccgCCCCTCCACAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                            <  1:3140147/145‑1 (MQ=255)
                                   cTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                           <  2:1917577/148‑1 (MQ=37)
                                    tGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                          <  1:2926967/148‑1 (MQ=37)
                                    ggCCCCTCCACAAGTAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGt                                                                                                          <  1:2926898/147‑1 (MQ=255)
                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAACAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  2:2553797/148‑1 (MQ=255)
                                      ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttt                                                                                                        <  2:1678778/148‑1 (MQ=37)
                                       cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGtttt                                                                                                       <  1:2166141/148‑1 (MQ=37)
                                        ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCCAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                      >  1:2383787/1‑148 (MQ=255)
                                        ccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGttttt                                                                                                      <  2:3160023/148‑1 (MQ=37)
                                              caaGTAAAAACCCGCCGCAGCGGGTTTTTACGTAACACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaa                                                                                                <  2:800635/147‑1 (MQ=255)
                                                attAAAAACCCGCCGCATCGGTTTATTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAc                                                                                              <  2:2212720/146‑1 (MQ=11)
                                                aGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                              >  1:1577087/1‑147 (MQ=37)
                                                 gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                             <  1:3194708/148‑1 (MQ=38)
                                                 ctAAAAACCCGCCGAAGCGGTTTTTTACGTAAAACAGTTGAAACTGAACTATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                                                                                             <  1:1035645/147‑1 (MQ=11)
                                                  tAAACACCCGCCGCAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                                                                                            <  2:2311735/148‑1 (MQ=255)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                           <  1:2369311/148‑1 (MQ=255)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                           <  1:2369388/148‑1 (MQ=37)
                                                   aaaaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAg                                                                                           <  1:2369313/148‑1 (MQ=37)
                                                       aCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGa                                                                                       >  2:508824/1‑148 (MQ=37)
                                                         ccgccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaaa                                                                                     <  1:3662725/148‑1 (MQ=37)
                                                             cgcatcGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                 <  2:3617639/143‑1 (MQ=255)
                                                             cgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa                                                                                 <  2:2882003/148‑1 (MQ=37)
                                                               aaGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc                                                                               >  2:129536/1‑148 (MQ=37)
                                                                aGCGGGTTAGTACGTACAACAGGTCAAACTGACCGATACGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCCCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCg                                                                              <  2:2120234/148‑1 (MQ=11)
                                                                 gCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGa                                                                             <  2:2383787/148‑1 (MQ=37)
                                                                      tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGc                                                                        >  2:3370264/1‑148 (MQ=12)
                                                                      tttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCAATTTTGAGCCCCACCAAAAGCAAAAACCCGCCGAAGCAGGTTTTTACGGAAAACAAGTGAAACTGACCGATaatc                                                                        >  1:1516519/1‑146 (MQ=11)
                                                                       ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGGCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACAGACCCATAAGcc                                                                       >  1:2071045/1‑148 (MQ=11)
                                                                       ttttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTATTCAGCCCCGCCAAAATCAAAAACCCGCCGAAGCGGGTTTTTAAGGAAAACAGGTGAAACTGACCGATAAGcc                                                                       >  1:193241/1‑148 (MQ=11)
                                                                        tttACGTAAAACAGGTGAAACTGACCGATAATCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCTTTCCTTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGTAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCg                                                                      <  2:915975/148‑1 (MQ=11)
                                                                         ttACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGc                                                                     <  2:2455521/148‑1 (MQ=255)
                                                                          tACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGct                                                                    <  2:3132442/148‑1 (MQ=38)
                                                                             gTACAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACTGTCATTATTTTTTCTGCCCCTACATAAGCACAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                 <  1:2812703/148‑1 (MQ=11)
                                                                             gTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                 <  1:3300769/148‑1 (MQ=255)
                                                                             gTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttc                                                                 <  1:3300603/148‑1 (MQ=255)
                                                                              tcaaaCACGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGctttct                                                                <  2:2347018/146‑1 (MQ=255)
                                                                              tAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTAAGTTAAACAGGTGAAACTGACCGATAAGCCGcgttct                                                                >  1:3472528/1‑148 (MQ=11)
                                                                                aaaCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGctttcttt                                                              <  2:3431997/148‑1 (MQ=255)
                                                                                  aCAGGTTAAACTGACCGATAAGCCTCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                            <  2:1300749/148‑1 (MQ=2)
                                                                                  aCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGGCAGGCATTAATTTTTCTGCCCCTCCAAAAGCAAAAAACCCCCCAAGCGGGTTTTTAAGTAAAACAGGTGAAACTGACCGATAATCCGCGTTGTTTTg                                                            >  1:370874/1‑148 (MQ=11)
                                                                                   cAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGcgttctgtggg                                                           >  2:3332646/1‑143 (MQ=255)
                                                                                     ggTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgtgt                                                         >  2:3250697/1‑148 (MQ=255)
                                                                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                      <  2:1127106/148‑1 (MQ=25)
                                                                                        gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                      <  1:70591/148‑1 (MQ=25)
                                                                                        caaaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGtata                                                      <  1:1686914/147‑1 (MQ=21)
                                                                                           aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCCTTCATTTTTAAGCCCCGCCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAATCCGCTTTCGTTTGGGGATAGcg                                                   >  1:2661909/1‑148 (MQ=11)
                                                                                           aCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcg                                                   <  2:2068861/148‑1 (MQ=25)
                                                                                            cTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTCTGCCCCTCCAAAACCAAAAAACCCCCGAACCGGGGTTTTACGTAAACCAGGGGAAAATGACCGATAACCCGCGTTCGGTTGGGTGGAGcgg                                                  >  1:466655/1‑147 (MQ=11)
                                                                                              gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAAcc                                                                                                                       >  1:3501796/1‑77 (MQ=21)
                                                                                              gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcg                                                <  1:2631264/148‑1 (MQ=37)
                                                                                              gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtcg                                                >  1:3501825/1‑148 (MQ=37)
                                                                                               aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGcgtatt                                               >  2:2468469/1‑145 (MQ=255)
                                                                                               aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGaa                                                                                                                >  2:2271753/1‑83 (MQ=33)
                                                                                                    gaaGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACa                                          <  2:3027211/147‑1 (MQ=37)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAATCCGCTTTCTTTTGGGTATAGAGTCGTGGCCAGtg                                       >  2:2082866/1‑147 (MQ=255)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTAAGTAAACCAGGTGAAACTGACCGATAAACCGCTTTCTTTTGGGGATAGCCTCGGGGACAGtc                                       >  2:3232846/1‑148 (MQ=11)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCAATTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTATTTTGGGTATAGCGTCGTGGACAGtc                                       >  1:935281/1‑148 (MQ=255)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCAATCATTTTTCTGCCCCACCAAAATCAAAAAACCGCCCAACCGGGTTTTTAAGTAAACCACGTGAAACTGACCCATAAGCCGCGTTCTTTTGGGTATAGCGTAGTGGGCACtc                                       >  1:1136448/1‑148 (MQ=11)
                                                                                                       gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTAATTCATTTTTGTGCCCCTCCAAAAGCAAATACCCGCCGGAGCGGGGTTTTACGTAAACCAGGTGTAACTGTCCGATAAGCCGGTTTGTTTTGGGGATAGCGTCGTGGACAGtg                                       >  1:935275/1‑147 (MQ=11)
                                                                                                         cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCGTTCTGTTGGGGATAGCGTGGTGGACAGtcat                                     >  2:2972287/1‑148 (MQ=11)
                                                                                                         cGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTAATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCCAAGCGGGTTTTTACGTAAAACAGGTGAAACAGACCGATAAGCCGCTTTCTTTTGGGTAGAGCGTAGTGGTCCGtcat                                     >  2:1064940/1‑148 (MQ=11)
                                                                                                         cGCCTTCTTTTGGTTATCGTCTCGTCGACAGTGCTTCAAATTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGTTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                     <  2:1766702/148‑1 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACCGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAACCCGGTTTTTAAGTAAACCAGCTGAAACTGACCGAAAAGCCGCTTTCTTTTGGGTATAGCGTAGTGGCCACtcaag                                    >  1:2005551/1‑146 (MQ=11)
                                                                                                          gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCCGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcatt                                    >  2:891426/1‑148 (MQ=255)
                                                                                                            tttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                  <  1:3313792/148‑1 (MQ=37)
                                                                                                            tttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                                                                                             >  2:919859/1‑73 (MQ=34)
                                                                                                             ttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattcct                                 >  2:3035312/1‑146 (MQ=37)
                                                                                                             ttctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCCCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcaataat                                 >  2:3217079/1‑143 (MQ=255)
                                                                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCt                               <  1:1229652/148‑1 (MQ=17)
                                                                                                               ctttTGGGTATAGTGTCGTGGACAGTCATTCATTTTTCTGCCCCTACAAAAGCAAAAAACCCCCGAAGCGGGTTTTTACGTAAACCACGTGAAACTGAACGATAAGCCCCGTTCTTTTGGGGATAGGGTCGAGGAAAGtcaatcaaaa                               >  1:3285120/1‑145 (MQ=11)
                                                                                                                 tttGTGTATAGTGTCGTGTACAGTCATTCTTATTTCTGCCCCTCCAAAAGAACATACCCGTCGATGGGTTTTTTTACGTATCCCATGTGTAACTGACCGATTAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTAATCttt                             <  1:2133468/148‑1 (MQ=11)
                                                                                                                            tgtCGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGCTAAGCCGCTTTCTTTTGGGTTTAGCGTCGTGGACAGGCAATCATATTTATTCCCCGCCa                  >  1:3750504/1‑148 (MQ=11)
                                                                                                                               cGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGGATAGCGTCGTGGACAGTCATTCATCATTATGCACCTCCaaaa               >  2:2293433/1‑148 (MQ=11)
                                                                                                                               cGTGGACAGTCATTCATTTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa               <  2:2584874/148‑1 (MQ=21)
                                                                                                                                   gACAGTCATTCATTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCaa           >  2:2695833/1‑148 (MQ=255)
                                                                                                                                        tcattcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc      <  2:3501825/148‑1 (MQ=37)
                                                                                                                                        tcattcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc      >  2:596973/1‑148 (MQ=37)
                                                                                                                                        tcattcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGGCATGCATCTGTCTGCCCCACCAAAAGCACTAAcc      >  2:596450/1‑148 (MQ=11)
                                                                                                                                        tcattcatTTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAAcc      <  2:1612506/148‑1 (MQ=21)
                                                                                                                                         cattcattTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTCTAGCGTCGTGGACAGTCAATCATCTTTGTGCCCCTGCAAAAGCAAAAAccc     >  1:3389849/1‑148 (MQ=11)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc  >  1:3169761/1‑148 (MQ=37)
                                                                                                                                            tcatTTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTGTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCTGCCCCTCCAAAATCAAAAACccgcc  >  2:2837417/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCC  >  CP009273/3263265‑3263552

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