Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,035,421 G→T 93.7% intergenic (‑62/‑8) mobA ← / → yihD molybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,035,4210GT93.7% 539.5 / 27.1 190intergenic (‑62/‑8)mobA/yihDmolybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base T (95/83);  total (99/91)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 195 total aligned reads displayed.

AAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAC  >  CP009273/4035276‑4035565
                                                                                                                                                 |                                                                                                                                                
aaGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCa                                                                                                                                                >  1:2213420/1‑148 (MQ=255)
aaGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCa                                                                                                                                                <  2:695996/148‑1 (MQ=255)
  gCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGatcatc                                                                                                                                              <  1:2344304/148‑1 (MQ=255)
  gCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGatcatc                                                                                                                                              >  1:556340/1‑148 (MQ=255)
   cAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAtcatca                                                                                                                                             <  1:449533/148‑1 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  1:86789/148‑1 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          >  1:3056209/1‑148 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  2:3425225/148‑1 (MQ=255)
          ttATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTCGGATTGACCCTGCTAACGTGTCTGTCTCAGGAGTAAGGATCATCACTATGaa                                                                                                                                      <  1:1173668/148‑1 (MQ=255)
            aTCTACGCCGCCCTTTCGTCTGGCTTTACCGCCTGCCAGCACAACGACTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAAt                                                                                                                                    <  2:3515495/148‑1 (MQ=255)
            aTCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAAt                                                                                                                                    <  1:1650618/148‑1 (MQ=255)
              cTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGt                                                                                                                                  >  2:3190335/1‑148 (MQ=255)
               tACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTa                                                                                                                                 <  1:3143077/148‑1 (MQ=255)
                aCGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGGTCATCACTATGAAATGTaa                                                                                                                                >  2:2663030/1‑148 (MQ=255)
                 cgccgcCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAAAGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTaaa                                                                                                                               >  1:1325017/1‑148 (MQ=255)
                    cgccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGt                                                                                                                            <  1:2023131/148‑1 (MQ=255)
                        cATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtga                                                                                                                        >  2:1188511/1‑148 (MQ=255)
                         aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaa                                                                                                                       >  1:102443/1‑148 (MQ=255)
                         aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAAAGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAAAGTCtgaa                                                                                                                       >  1:2876002/1‑148 (MQ=255)
                          ttCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaat                                                                                                                      >  1:832164/1‑148 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                     >  2:2526531/1‑148 (MQ=255)
                           tCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatg                                                                                                                     >  2:737579/1‑148 (MQ=255)
                              tCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAg                                                                                                                  >  2:433118/1‑148 (MQ=255)
                               cTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGt                                                                                                                 <  1:715814/148‑1 (MQ=255)
                                 ggCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                               >  2:1614129/1‑148 (MQ=255)
                                   cTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAtt                                                                                                             >  2:2391415/1‑148 (MQ=255)
                                     ttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGa                                                                                                           >  2:1737338/1‑148 (MQ=255)
                                     ttACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGa                                                                                                           >  2:227801/1‑148 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:616765/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:2827311/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  1:433839/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:2003100/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:2003286/148‑1 (MQ=255)
                                        ccGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAAct                                                                                                        <  1:675/148‑1 (MQ=255)
                                        ccGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAAct                                                                                                        >  2:258391/1‑148 (MQ=255)
                                         cGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATGCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActc                                                                                                       >  1:3378219/1‑148 (MQ=255)
                                            cTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcctc                                                                                                    <  1:461517/148‑1 (MQ=255)
                                                cAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagc                                                                                                <  1:251448/148‑1 (MQ=255)
                                                 aGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGGTATTGAACTCCTccagcc                                                                                               >  2:2138557/1‑148 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                            <  2:166814/148‑1 (MQ=255)
                                                        cGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTgg                                                                                        >  1:1184143/1‑148 (MQ=255)
                                                         gCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTTGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTTAATGAAGTTATTTAACTCCACCATCCAGCCAGGc                                                                                       >  1:3384521/1‑148 (MQ=255)
                                                                aTCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAaga                                                                                <  2:1482429/148‑1 (MQ=255)
                                                                        catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        >  2:2254628/1‑148 (MQ=255)
                                                                        catcaGAGTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        <  2:375821/148‑1 (MQ=255)
                                                                         atcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCt                                                                       <  1:1948781/148‑1 (MQ=255)
                                                                          tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAGGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                      >  1:1812126/1‑148 (MQ=255)
                                                                          tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                      >  1:1457642/1‑148 (MQ=255)
                                                                            aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTTTCTTTCTCAAGATTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTaa                                                                    <  1:3291546/148‑1 (MQ=255)
                                                                            aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTaa                                                                    <  2:685385/148‑1 (MQ=255)
                                                                               ttCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACct                                                                 >  2:852902/1‑148 (MQ=255)
                                                                                 cACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgc                                                               >  1:3131369/1‑148 (MQ=255)
                                                                                   ccGATATCGCCTCGTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgctg                                                             >  1:324445/1‑148 (MQ=255)
                                                                                    cGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCtgctgc                                                            >  2:770874/1‑148 (MQ=255)
                                                                                     gATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCa                                                           <  1:1066591/148‑1 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         >  2:3350808/1‑148 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         >  2:3414174/1‑148 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         >  2:1117947/1‑148 (MQ=255)
                                                                                         tCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAttt                                                       >  1:3476452/1‑148 (MQ=255)
                                                                                          cGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAtttt                                                      >  2:3449411/1‑148 (MQ=255)
                                                                                          cGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAtttt                                                      >  2:1016559/1‑148 (MQ=255)
                                                                                           gCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCAGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAttttt                                                     >  2:3279729/1‑148 (MQ=255)
                                                                                            ccTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTg                                                    >  1:716910/1‑148 (MQ=255)
                                                                                             ctctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGc                                                   >  2:2996960/1‑148 (MQ=255)
                                                                                             ctctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGc                                                   <  1:3442913/148‑1 (MQ=255)
                                                                                              tctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCa                                                  <  2:1073488/148‑1 (MQ=255)
                                                                                                 tttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaa                                               <  1:2752251/148‑1 (MQ=255)
                                                                                                  ttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaaa                                              >  2:71135/1‑148 (MQ=255)
                                                                                                    aTTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACt                                            <  1:1112589/148‑1 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                           >  2:457025/1‑148 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                           >  2:607200/1‑148 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                           <  1:1217831/148‑1 (MQ=255)
                                                                                                      tgtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                          <  2:2425178/148‑1 (MQ=255)
                                                                                                      tgtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTgg                                          <  2:1593911/148‑1 (MQ=255)
                                                                                                       gtgGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGc                                         <  2:1347310/148‑1 (MQ=255)
                                                                                                         gggATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGa                                       >  1:1595085/1‑148 (MQ=255)
                                                                                                            aTTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAag                                    >  2:305196/1‑148 (MQ=255)
                                                                                                             ttGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAaga                                   >  1:92650/1‑148 (MQ=255)
                                                                                                             ttGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAaga                                   <  2:556340/148‑1 (MQ=255)
                                                                                                                aCCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCCATTGTTGCAGAAACTGGCGCACGAGTc                                >  2:1865217/1‑148 (MQ=255)
                                                                                                                aCCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTc                                >  2:1865195/1‑148 (MQ=255)
                                                                                                                 cccTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCa                               <  2:1347988/148‑1 (MQ=255)
                                                                                                                         aCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGa                       >  1:2272732/1‑148 (MQ=255)
                                                                                                                         aCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGa                       >  1:2661212/1‑148 (MQ=255)
                                                                                                                             gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGc                   >  2:3315760/1‑148 (MQ=255)
                                                                                                                              tctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCg                  >  2:1315336/1‑148 (MQ=255)
                                                                                                                                tgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGaa                <  1:524074/148‑1 (MQ=255)
                                                                                                                                 gtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAAc               <  2:2294936/148‑1 (MQ=255)
                                                                                                                                  tctcAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACt              <  1:2549632/148‑1 (MQ=255)
                                                                                                                                     cAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGc           >  2:432096/1‑148 (MQ=255)
                                                                                                                                       agagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCgg         <  1:2286959/148‑1 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        >  2:1721738/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        >  2:2724258/1‑148 (MQ=255)
                                                                                                                                         agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAt       <  1:737579/148‑1 (MQ=255)
                                                                                                                                         agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAt       <  1:2974413/148‑1 (MQ=255)
                                                                                                                                         agTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAt       <  2:2747513/148‑1 (MQ=255)
                                                                                                                                            aaGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTg    <  2:1457642/148‑1 (MQ=255)
                                                                                                                                             aGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGa   <  1:2663030/148‑1 (MQ=255)
                                                                                                                                             aGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGa   >  2:3149593/1‑148 (MQ=255)
                                                                                                                                              ggAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAc  >  1:1552112/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAC  >  CP009273/4035276‑4035565

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: