Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,935,100 G→T 100% intergenic (‑413/‑68) yieP ← / → rrsC putative transcriptional regulator/16S ribosomal RNA of rrnC operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,935,1000GT95.1% 212.7 / 6.5 102intergenic (‑413/‑68)yieP/rrsCputative transcriptional regulator/16S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base T (48/49);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
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 102 total aligned reads displayed.

ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAA  >  CP009273/3934955‑3935161
                                                                                                                                                 |                                                             
atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                             >  1:2933709/1‑148 (MQ=21)
 tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAAt                                                            >  2:2497742/1‑148 (MQ=21)
 tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAAt                                                            >  2:2497739/1‑148 (MQ=21)
  gCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATa                                                           <  2:1036796/148‑1 (MQ=21)
      gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACCATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCaa                                                       >  2:1075776/1‑148 (MQ=21)
      gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCaa                                                       <  1:1378212/148‑1 (MQ=21)
       aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                      >  1:2491309/1‑148 (MQ=21)
       aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                      >  2:2995054/1‑148 (MQ=21)
       aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAg                                                      >  1:3546032/1‑148 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGAAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                     >  2:3087010/1‑148 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGAAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                     >  2:3086701/1‑148 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGACAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAATATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                     <  1:2723156/148‑1 (MQ=255)
         tctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtc                                                    <  2:52527/148‑1 (MQ=21)
         tctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtc                                                    >  2:377281/1‑148 (MQ=21)
          ctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtct                                                   <  1:113215/148‑1 (MQ=21)
            gTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCa                                                 >  1:1295361/1‑148 (MQ=21)
            gTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCa                                                 <  2:1443844/148‑1 (MQ=21)
              agcgttaaaGCGTATTATGCACACCCCGGGCCGCTGATAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                               <  1:221224/142‑1 (MQ=11)
              aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAag                                               <  1:221295/148‑1 (MQ=34)
                cGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAagag                                             >  2:2530107/1‑148 (MQ=21)
                  ggAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTg                                           <  1:2004795/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGATAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  1:2494777/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  1:1271047/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  1:1191962/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  1:1271034/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          >  2:239376/1‑148 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  2:3516282/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                          <  2:2826383/148‑1 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATATATCATACAATCTCTGTGGGCACTCTAAGATACGTATTCTTAACGTCGCAAGACGAAAAATTAATAACAATTCTCAAGAGTGa                                          >  2:239462/1‑148 (MQ=11)
                    aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGaa                                         >  2:683094/1‑148 (MQ=21)
                     aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCACAAGAGTGAac                                        <  1:569536/148‑1 (MQ=16)
                     aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                        <  1:2057453/148‑1 (MQ=18)
                     aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                        >  1:641394/1‑148 (MQ=18)
                      aGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAaca                                       >  2:842477/1‑148 (MQ=17)
                      aGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAaca                                       >  2:842474/1‑148 (MQ=17)
                      aGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAaca                                       >  1:1045017/1‑148 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                     <  1:3517019/148‑1 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                     <  2:2872250/148‑1 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                     <  1:1439584/148‑1 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                     <  1:3518198/148‑1 (MQ=17)
                          tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                   <  2:1707457/148‑1 (MQ=17)
                          tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGACACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                   >  2:184604/1‑148 (MQ=17)
                          tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGACACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                   >  2:184579/1‑148 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                  <  1:2810995/148‑1 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACCATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCTAGTCGAAAAATTAATACCAAGTCTCAAGAGTGAACACGtta                                  >  1:584418/1‑146 (MQ=255)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                  >  2:1182360/1‑148 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                  >  2:1182414/1‑148 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                  >  1:1260708/1‑148 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTaa                                  >  2:376863/1‑148 (MQ=32)
                            ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                 <  2:1772102/148‑1 (MQ=17)
                            ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                 <  2:2665258/148‑1 (MQ=17)
                            ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                 <  2:1771572/148‑1 (MQ=17)
                            ttatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                 <  2:3412460/148‑1 (MQ=17)
                             tatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                                <  2:1910329/148‑1 (MQ=17)
                                gCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAt                             <  2:2718722/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                            <  1:3364730/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                            <  1:875472/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                            <  1:875447/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                            <  1:354557/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                            <  1:607041/148‑1 (MQ=32)
                                  acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                           <  2:2833917/148‑1 (MQ=17)
                                   cacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTAc                          >  1:1905654/1‑148 (MQ=17)
                                    acCCCGCGCCTCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACg                         <  1:1969698/148‑1 (MQ=17)
                                     ccccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGa                        <  1:3321845/148‑1 (MQ=17)
                                      cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                       <  2:3478927/148‑1 (MQ=32)
                                       ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                      >  2:2773004/1‑148 (MQ=17)
                                       ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                      <  2:486848/148‑1 (MQ=17)
                                       ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                      <  2:486777/148‑1 (MQ=17)
                                       ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                      >  2:1737288/1‑148 (MQ=17)
                                         gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGtt                    <  2:2736376/148‑1 (MQ=17)
                                         gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGtt                    <  2:2736519/148‑1 (MQ=17)
                                         gcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGtt                    >  1:1490433/1‑148 (MQ=17)
                                           gccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTa                  <  1:1746776/148‑1 (MQ=17)
                                               cTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc              >  1:3170638/1‑148 (MQ=17)
                                               cTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc              <  2:3001779/148‑1 (MQ=17)
                                               cTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc              >  1:3170586/1‑148 (MQ=17)
                                                tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt             <  1:859829/148‑1 (MQ=17)
                                                tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt             <  2:3301206/148‑1 (MQ=17)
                                                tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt             >  2:312006/1‑148 (MQ=17)
                                                tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt             >  2:311980/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:3259977/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:1491042/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:3457965/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:373276/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:373311/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt           >  2:373338/1‑148 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg          >  2:1062410/1‑148 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg          >  2:1062980/1‑148 (MQ=17)
                                                    aaaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGa         >  1:2963292/1‑148 (MQ=17)
                                                    aaaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGa         >  2:377480/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg        <  2:1466491/148‑1 (MQ=17)
                                                      aaagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGc       <  2:3227782/148‑1 (MQ=17)
                                                        agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCAt     >  1:495817/1‑148 (MQ=17)
                                                         gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATc    >  1:1701781/1‑148 (MQ=17)
                                                         gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATc    >  1:1702064/1‑148 (MQ=17)
                                                          caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa   >  1:912687/1‑148 (MQ=14)
                                                          caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa   >  1:1137919/1‑148 (MQ=14)
                                                           aaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaa  <  2:2998558/148‑1 (MQ=11)
                                                           aaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaa  <  2:729524/148‑1 (MQ=11)
                                                           aaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCaa  <  2:729529/148‑1 (MQ=11)
                                                                                                                                                 |                                                             
ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAA  >  CP009273/3934955‑3935161

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