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,1000GT98.0% 226.5 / ‑1.2 100intergenic (‑413/‑68)yieP/rrsCputative transcriptional regulator/16S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (46/52);  total (48/52)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCA  >  CP009273/3934955‑3935160
                                                                                                                                                 |                                                            
atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACCCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                            <  1:2716192/148‑1 (MQ=21)
atGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTaa                                                            >  1:1199307/1‑148 (MQ=21)
 tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAATCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAAt                                                           >  1:2519525/1‑148 (MQ=21)
 tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTGAACGTCGCAAGACGATAAATTAAt                                                           >  1:2519462/1‑148 (MQ=255)
  gCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATa                                                          <  1:1083287/148‑1 (MQ=21)
   cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATAc                                                         >  1:2610358/1‑148 (MQ=21)
    ttGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCACGACGAAAAATTAATAcc                                                        >  2:2396639/1‑148 (MQ=21)
    ttGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATAcc                                                        >  1:1966405/1‑148 (MQ=21)
     tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                       <  2:2045784/148‑1 (MQ=21)
     tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                       <  2:2044880/148‑1 (MQ=21)
     tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCa                                                       <  2:2797852/148‑1 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                    >  1:992144/1‑148 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                    <  2:2717969/148‑1 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                    <  2:1572909/148‑1 (MQ=21)
        ctctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGt                                                    <  2:204934/148‑1 (MQ=21)
          ctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtct                                                  <  1:2465702/148‑1 (MQ=21)
          ctGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGtct                                                  <  1:1412411/148‑1 (MQ=21)
             tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCaa                                               >  1:3246690/1‑148 (MQ=21)
             tAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCaa                                               >  1:3246766/1‑148 (MQ=21)
              aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAag                                              <  2:1419274/148‑1 (MQ=21)
              aGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAag                                              <  2:2173919/148‑1 (MQ=21)
               gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCaa                                                      <  2:129393/139‑1 (MQ=21)
               gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCaa                                                      >  1:129393/1‑139 (MQ=21)
               gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                             <  1:1899836/148‑1 (MQ=21)
               gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                             <  2:355029/148‑1 (MQ=21)
               gCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAaga                                             >  1:3444826/1‑148 (MQ=21)
                cGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAagag                                            >  1:3331442/1‑148 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                         >  2:2261972/1‑148 (MQ=21)
                   gAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                         >  2:467872/1‑148 (MQ=21)
                   gAAAGCGTATTACGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGa                                         <  1:1919199/148‑1 (MQ=21)
                    aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGaa                                        <  2:1687085/148‑1 (MQ=21)
                    aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGaa                                        >  1:1734975/1‑148 (MQ=21)
                    aaaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGaa                                        <  2:232720/148‑1 (MQ=21)
                     aaGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAac                                       >  1:2109690/1‑148 (MQ=18)
                       gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                     >  1:2775603/1‑148 (MQ=17)
                       gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                     >  1:2778287/1‑148 (MQ=17)
                       gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAacac                                     >  2:342380/1‑148 (MQ=17)
                       gCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAacac                                     >  2:1300482/1‑148 (MQ=32)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                    >  1:709924/1‑148 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                    >  1:1171517/1‑148 (MQ=17)
                        cGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACg                                    >  1:2634326/1‑148 (MQ=17)
                          tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                  <  1:437014/148‑1 (MQ=17)
                          tattatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTa                                  <  2:2545784/148‑1 (MQ=17)
                           attatGCACCCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                 <  1:3081613/148‑1 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                 <  2:892533/148‑1 (MQ=17)
                           attatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTaa                                 <  2:723320/148‑1 (MQ=17)
                            ttatGCACACCCCGCGCCGCTGATAAAAAGCAACGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAt                                <  1:710035/148‑1 (MQ=12)
                             tatGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAAtt                               <  1:1269948/148‑1 (MQ=17)
                               tGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACTAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCa                             >  1:150208/1‑148 (MQ=17)
                               tGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCa                             <  1:2936527/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                           <  2:2974637/148‑1 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                           >  2:14075/1‑148 (MQ=17)
                                 cacaCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCAtt                           <  2:2974618/148‑1 (MQ=17)
                                  acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                          <  2:292453/148‑1 (MQ=17)
                                  acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                          <  2:2778246/148‑1 (MQ=17)
                                  acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                          <  2:291477/148‑1 (MQ=17)
                                  acacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTa                          <  2:291473/148‑1 (MQ=17)
                                     ccccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGa                       <  1:3256474/148‑1 (MQ=17)
                                     ccccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGa                       <  1:2478203/148‑1 (MQ=17)
                                      cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                      <  2:2610358/148‑1 (MQ=17)
                                      cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                      <  1:276409/148‑1 (MQ=17)
                                      cccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                      <  1:276375/148‑1 (MQ=17)
                                      ccagcgcCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGaa                      >  2:3194237/4‑148 (MQ=17)
                                       ccGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAg                     <  1:843843/148‑1 (MQ=17)
                                        cgcgCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGt                    <  2:133850/148‑1 (MQ=17)
                                          cgccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAATTCTCAAGAGTGAACACGTAATTAATTACGAAGttt                  >  1:132430/1‑148 (MQ=12)
                                           gccgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTa                 <  2:1082493/148‑1 (MQ=17)
                                             cgcTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAt               <  2:965111/148‑1 (MQ=17)
                                              gcTGATAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAAtt              <  1:1870104/148‑1 (MQ=17)
                                               cTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTc             >  2:3357839/1‑148 (MQ=17)
                                                tGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCt            >  2:229905/1‑148 (MQ=17)
                                                 gagaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCtt           >  2:1736187/1‑148 (MQ=17)
                                                 gagaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCtt           >  2:2478133/1‑148 (MQ=17)
                                                  agaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCttt          >  1:585319/1‑148 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg         <  2:999140/148‑1 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg         <  1:2031070/148‑1 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg         <  1:231659/148‑1 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg         <  1:563808/148‑1 (MQ=17)
                                                   gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTg         >  2:2783181/1‑148 (MQ=17)
                                                    aaaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGa        >  1:151893/1‑148 (MQ=17)
                                                    aaaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGa        <  1:799408/148‑1 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  1:476841/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:605544/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:518980/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:2978516/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:2699052/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  1:1557747/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:1264127/1‑148 (MQ=17)
                                                     aaaaGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAAAACGTAATTCATTACGAAGTTTAATTCTTTGAg       >  2:1264128/1‑148 (MQ=17)
                                                       aagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa     >  2:722047/1‑148 (MQ=17)
                                                       aagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa     <  1:2780374/148‑1 (MQ=17)
                                                       aagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa     <  1:2780221/148‑1 (MQ=17)
                                                       aagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCa     <  1:3213591/148‑1 (MQ=17)
                                                        agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCAt    >  1:3090905/1‑148 (MQ=17)
                                                        agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCAt    <  1:1488049/148‑1 (MQ=17)
                                                        agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCAt    <  2:1966405/148‑1 (MQ=17)
                                                         gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATc   >  2:2846247/1‑148 (MQ=17)
                                                         gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATc   >  2:2751137/1‑148 (MQ=17)
                                                         gcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATc   >  2:3301885/1‑148 (MQ=32)
                                                          caaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATTAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCa  >  2:989211/1‑148 (MQ=14)
                                                                                                                                                 |                                                            
ATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCA  >  CP009273/3934955‑3935160

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: