Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,915,998 C→A 100% A116S (GCT→TCT)  yqcA ← short‑chain flavodoxin, FMN‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,915,9980CA96.5% 615.6 / 12.1 200A116S (GCT→TCT) yqcAshort‑chain flavodoxin, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (93/100);  total (97/103)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 201 total aligned reads displayed.

TTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAA  >  CP009273/2915853‑2916141
                                                                                                                                                 |                                                                                                                                               
ttcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAgagc                                                                                                                                               >  1:351417/1‑148 (MQ=255)
    gCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTg                                                                                                                                           >  2:3156529/1‑148 (MQ=255)
       ccGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTc                                                                                                                                        <  2:1679917/148‑1 (MQ=255)
        cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTcc                                                                                                                                       >  2:344562/1‑148 (MQ=255)
        cGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTcc                                                                                                                                       >  1:627174/1‑148 (MQ=255)
            aaaaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGca                                                                                                                                   <  1:2232414/148‑1 (MQ=255)
             aaaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGcaa                                                                                                                                  >  2:1714739/1‑148 (MQ=255)
                 gACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAg                                                                                                                              <  2:2245015/148‑1 (MQ=255)
                 gACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAACATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAg                                                                                                                              >  1:2371843/1‑148 (MQ=255)
                   cATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGc                                                                                                                            >  2:727779/1‑148 (MQ=255)
                     tttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCAt                                                                                                                          >  1:519137/1‑148 (MQ=255)
                      ttACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATc                                                                                                                         >  2:3462431/1‑148 (MQ=255)
                         ccGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGaa                                                                                                                      <  1:1431956/148‑1 (MQ=255)
                         ccGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGaa                                                                                                                      <  1:1457114/148‑1 (MQ=255)
                           gctcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACt                                                                                                                    <  2:2163815/148‑1 (MQ=255)
                            ctcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTg                                                                                                                   <  2:817193/148‑1 (MQ=255)
                             tcgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGt                                                                                                                  <  1:1538659/148‑1 (MQ=255)
                              cgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                 <  2:106420/148‑1 (MQ=255)
                              cgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                 <  1:3131902/148‑1 (MQ=255)
                              cgctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGtt                                                                                                                 >  1:2074228/1‑148 (MQ=255)
                                ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTg                                                                                                               <  1:1751013/148‑1 (MQ=255)
                                ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTg                                                                                                               >  2:1622813/1‑148 (MQ=255)
                                ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTg                                                                                                               >  1:341291/1‑148 (MQ=255)
                                ctGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATACCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAAACCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTCTTTg                                                                                                               >  1:2501295/1‑148 (MQ=255)
                                 tGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgc                                                                                                              <  1:3156529/148‑1 (MQ=255)
                                    aCATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgc                                                                                                           <  1:1884956/148‑1 (MQ=255)
                                    aCATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgc                                                                                                           <  1:1027527/148‑1 (MQ=255)
                                     cATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgcc                                                                                                          >  2:2722081/1‑148 (MQ=255)
                                     cATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTgccgcc                                                                                                          >  2:938228/1‑148 (MQ=255)
                                      atatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                         <  2:2352323/148‑1 (MQ=255)
                                      atatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                         <  1:2192789/148‑1 (MQ=255)
                                      atatCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCa                                                                                                         <  2:560759/148‑1 (MQ=255)
                                         tCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTg                                                                                                      <  1:2557706/148‑1 (MQ=255)
                                         tCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTg                                                                                                      >  1:933320/1‑148 (MQ=255)
                                         tCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTg                                                                                                      >  2:335887/1‑148 (MQ=255)
                                           aGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCa                                                                                                    <  2:519137/148‑1 (MQ=255)
                                             gACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGa                                                                                                  >  1:2851163/1‑148 (MQ=255)
                                             gACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGa                                                                                                  >  1:1070872/1‑148 (MQ=255)
                                               caacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                <  1:3269661/148‑1 (MQ=255)
                                               caacaGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGaaa                                                                                                <  2:1511739/148‑1 (MQ=255)
                                                  caGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTc                                                                                             <  1:3493880/148‑1 (MQ=255)
                                                    gCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGTATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTcac                                                                                           <  1:1243337/148‑1 (MQ=255)
                                                    gCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTcac                                                                                           <  2:356405/148‑1 (MQ=255)
                                                    gCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTcac                                                                                           <  1:956155/148‑1 (MQ=255)
                                                    gCGTGCCCCACTGTTCGACCCACGGATTTGAGTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTcac                                                                                           <  1:772986/148‑1 (MQ=255)
                                                     cGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAACTTcaca                                                                                          >  2:3053603/1‑148 (MQ=255)
                                                          cccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAact                                                                                     <  2:3135938/148‑1 (MQ=255)
                                                           ccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAacta                                                                                    <  1:148898/148‑1 (MQ=255)
                                                           ccACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAacta                                                                                    <  1:1949233/148‑1 (MQ=255)
                                                            cACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAactac                                                                                   >  1:2965727/1‑148 (MQ=255)
                                                              cTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                 <  1:2280531/148‑1 (MQ=255)
                                                              cTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTAAAGGTAGTCAATCAAAATCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTg                                                                                 >  1:3497602/1‑148 (MQ=255)
                                                                 ttCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCg                                                                              >  1:984823/1‑148 (MQ=255)
                                                                 ttCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCAGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCTCATAACTACTGTCg                                                                              >  1:2088428/1‑148 (MQ=255)
                                                                     aCCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGa                                                                          >  1:1969454/1‑148 (MQ=255)
                                                                          cGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                     <  2:2052777/148‑1 (MQ=255)
                                                                          cGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                     <  2:358125/148‑1 (MQ=255)
                                                                          cGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCa                                                                     <  2:701614/148‑1 (MQ=255)
                                                                             aTTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATc                                                                  <  1:1080075/148‑1 (MQ=255)
                                                                             aTTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAAATCACATAACTACTGTCGCCGAGCGCAATc                                                                  >  2:1793592/1‑148 (MQ=255)
                                                                               ttGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCTCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCAc                                                                <  2:1357278/148‑1 (MQ=255)
                                                                                 gATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGc                                                              <  2:215546/148‑1 (MQ=255)
                                                                                      cGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATaa                                                         <  1:2591508/148‑1 (MQ=255)
                                                                                      cGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATaa                                                         <  1:899702/148‑1 (MQ=255)
                                                                                      cGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATaa                                                         <  1:807292/148‑1 (MQ=255)
                                                                                        tttCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACg                                                       >  1:1032170/1‑148 (MQ=255)
                                                                                         ttccggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGc                                                      >  2:2603876/1‑148 (MQ=255)
                                                                                           ccggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAg                                                    <  1:2954194/148‑1 (MQ=255)
                                                                                            cggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGa                                                   <  2:3544784/148‑1 (MQ=255)
                                                                                             ggttccggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGAt                                                  >  1:459865/1‑148 (MQ=255)
                                                                                                  cggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGc                                             <  1:3331210/148‑1 (MQ=255)
                                                                                                  cggGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGc                                             <  1:1678665/148‑1 (MQ=255)
                                                                                                      ttCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGa                                         <  1:1317478/148‑1 (MQ=255)
                                                                                                          cACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                     >  2:1150146/1‑148 (MQ=255)
                                                                                                          cACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                     <  2:872071/148‑1 (MQ=255)
                                                                                                          cACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                     <  2:2072591/148‑1 (MQ=255)
                                                                                                          cACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGcc                                     >  2:1150043/1‑148 (MQ=255)
                                                                                                            cTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCa                                   >  1:2040161/1‑148 (MQ=255)
                                                                                                                cGTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTACCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACt                               >  1:2045168/1‑148 (MQ=255)
                                                                                                                 gTCAATCAACAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTa                              <  2:2437497/148‑1 (MQ=255)
                                                                                                                   caatcaaCAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATc                            >  1:1850950/1‑148 (MQ=255)
                                                                                                                    aatcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCt                           <  1:1714739/148‑1 (MQ=255)
                                                                                                                    aatcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCt                           >  2:1901311/1‑148 (MQ=255)
                                                                                                                     atcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCtt                          >  2:547726/1‑148 (MQ=255)
                                                                                                                     atcaaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCtt                          >  1:3195329/1‑148 (MQ=255)
                                                                                                                       caaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTg                        >  1:2253839/1‑148 (MQ=255)
                                                                                                                        aaCAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGa                       <  1:2817598/148‑1 (MQ=255)
                                                                                                                          cAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAtt                     <  1:2722081/148‑1 (MQ=255)
                                                                                                                          cAGCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGAtt                     <  1:938228/148‑1 (MQ=255)
                                                                                                                            gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                   <  2:351417/148‑1 (MQ=255)
                                                                                                                            gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                   >  1:273599/1‑148 (MQ=255)
                                                                                                                            gCATTTCGCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTcc                   >  1:400887/1‑148 (MQ=255)
                                                                                                                                   gCCAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAag            >  2:2571613/1‑148 (MQ=255)
                                                                                                                                    ccAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAaga           >  1:3018821/1‑148 (MQ=255)
                                                                                                                                     cAACCCGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAagag          <  2:2198231/148‑1 (MQ=255)
                                                                                                                                        cccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTgg       >  1:1183294/1‑148 (MQ=255)
                                                                                                                                        cccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTgg       >  2:724085/1‑148 (MQ=255)
                                                                                                                                         ccGCTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATCACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGc      <  1:1260254/148‑1 (MQ=255)
                                                                                                                                            cTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGcaca   <  2:2549621/148‑1 (MQ=255)
                                                                                                                                            cTGAGAGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGcaca   <  2:1255677/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAACTGTTTGCCGCCATTGCAGAAATTCACATAACTACTGTCGCCGAGCGCAATCACGCCATAACGCAGATTCGGCTGGAAGCCCAGACTATCTTTGATTCCCTGAAAGAGTGGCACAA  >  CP009273/2915853‑2916141

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: ‑