Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,035,421 G→T 100% 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,4210GT97.8% 728.0 / 5.7 229intergenic (‑62/‑8)mobA/yihDmolybdopterin‑guanine dinucleotide synthase/DUF1040 protein YihD
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base T (108/116);  total (113/116)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 240 total aligned reads displayed.

AAGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAC  >  CP009273/4035276‑4035565
                                                                                                                                                 |                                                                                                                                                
aaGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCa                                                                                                                                                <  2:959896/148‑1 (MQ=255)
aaGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCa                                                                                                                                                <  1:2719183/148‑1 (MQ=255)
 aGCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCAt                                                                                                                                               >  1:2291965/1‑148 (MQ=255)
  gCACTCCTTTATCTACGCCGCCCAGTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACAGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGatcatc                                                                                                                                              <  2:2022520/148‑1 (MQ=255)
  gCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGatcatc                                                                                                                                              <  1:1973882/148‑1 (MQ=255)
  gCAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTGTTATTGTGGGATTGACCCTGCTCACGTGTCTGTCTCAAGAGTAAGGatcatc                                                                                                                                              >  2:1746508/1‑148 (MQ=255)
  gCAATCCTTTATCTACGCCGCCCAATCGTCTGGCTTTAACGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGTATTGACCCTGATAACGTGTCTTTCTCAAGAGTAAGGatcatc                                                                                                                                              >  2:3654648/1‑148 (MQ=255)
   cAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAtcatca                                                                                                                                             >  2:458406/1‑148 (MQ=255)
   cAATCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAtcatca                                                                                                                                             <  2:174843/148‑1 (MQ=255)
   cAATCCTTTATCTACGCCGCCCAGTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAtcatca                                                                                                                                             <  2:175755/148‑1 (MQ=255)
     aTCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACt                                                                                                                                           >  2:2918915/1‑148 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  1:172490/148‑1 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  1:2022737/148‑1 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  2:3603875/148‑1 (MQ=255)
      tCCTTTATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTa                                                                                                                                          <  2:2138902/148‑1 (MQ=255)
         tttATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGa                                                                                                                                       >  1:527640/1‑148 (MQ=255)
         tttATCTACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGa                                                                                                                                       >  2:2203284/1‑148 (MQ=255)
          ttATATACGCCGCCCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGaa                                                                                                                                      <  1:141230/148‑1 (MQ=255)
                    cgccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGt                                                                                                                            <  2:589314/148‑1 (MQ=255)
                     gccCATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTc                                                                                                                           <  1:908174/148‑1 (MQ=255)
                      cccATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCt                                                                                                                          <  2:3127964/148‑1 (MQ=255)
                      cccATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCt                                                                                                                          <  2:3808824/148‑1 (MQ=255)
                        cATTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtga                                                                                                                        <  2:3153731/148‑1 (MQ=255)
                         aTTCGTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaa                                                                                                                       >  1:3013287/1‑148 (MQ=255)
                             gTCTGGCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCtgaatgaa                                                                                                                   <  1:1386635/148‑1 (MQ=255)
                                 ggCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTa                                                                                                               <  2:257510/148‑1 (MQ=255)
                                  gCTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAt                                                                                                              >  2:1495576/1‑148 (MQ=255)
                                   cTTTACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTAtt                                                                                                             <  1:284744/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:2896280/148‑1 (MQ=255)
                                      tACCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGaa                                                                                                          <  2:3753013/148‑1 (MQ=255)
                                       aCCGCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGACTGAAGTTATTGAAc                                                                                                         >  1:967796/1‑148 (MQ=255)
                                          gCCTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcc                                                                                                      >  1:2122967/1‑148 (MQ=255)
                                            cTGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAActcctc                                                                                                    <  2:117681/148‑1 (MQ=255)
                                             tGCCAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACtcctcc                                                                                                   <  2:2971787/148‑1 (MQ=255)
                                                cAGCACAACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagc                                                                                                >  2:3451375/1‑148 (MQ=255)
                                                   cacaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTTACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTccagccag                                                                                             >  1:1624401/1‑148 (MQ=255)
                                                    acaACGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCcagccagc                                                                                            <  2:134203/148‑1 (MQ=255)
                                                       aCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTg                                                                                         >  1:553270/1‑148 (MQ=255)
                                                       aCGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTg                                                                                         >  1:631063/1‑148 (MQ=255)
                                                        cGCCTGTTATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTgg                                                                                        <  1:2325903/148‑1 (MQ=255)
                                                               tATCGTCGTCATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAag                                                                                 >  1:2542245/1‑148 (MQ=255)
                                                                     cgtcATCAGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAg                                                                           <  1:109293/148‑1 (MQ=255)
                                                                        catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        >  2:3756415/1‑148 (MQ=255)
                                                                        catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        >  2:1047739/1‑148 (MQ=255)
                                                                        catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        >  2:120146/1‑148 (MQ=255)
                                                                        catcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGAcc                                                                        >  2:2053072/1‑148 (MQ=255)
                                                                         ataaGATTCTCCTATATAGTCTCTTTTATTGTGGGATTGACCCTGCTAACTTGTCTGTCTGAAGAGTTAGTATCATCACTATTAAATGTAAACTTCTGAATGAAGTTATTGATCTCCTCCAGCCTGCCTGGCAAAAAGAGCAAGACCt                                                                       <  2:2431299/145‑1 (MQ=255)
                                                                          tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                      >  2:2288826/1‑148 (MQ=255)
                                                                          tcaGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCtt                                                                      >  2:2288707/1‑148 (MQ=255)
                                                                           caGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCaaaa                                                                                   <  1:1713888/134‑1 (MQ=255)
                                                                            aGATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTaa                                                                    >  1:283169/1‑148 (MQ=255)
                                                                             gATTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAc                                                                   <  2:917052/148‑1 (MQ=255)
                                                                              aTTCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAAcc                                                                  >  1:3580928/1‑148 (MQ=255)
                                                                               ttCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACct                                                                 >  2:710818/1‑148 (MQ=255)
                                                                                tCACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctg                                                                <  2:1840395/148‑1 (MQ=255)
                                                                                 cACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgc                                                               <  2:2752929/148‑1 (MQ=255)
                                                                                 cACCGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACctgc                                                               <  1:698534/148‑1 (MQ=255)
                                                                                   ccGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCATGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGGTATTGAACTCCTACAGCCAGCCTGGCACAAAGAGCCAGACCTTAACctgctg                                                             >  2:2183748/1‑148 (MQ=255)
                                                                                    cGATATCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCtgctgc                                                            >  2:472548/1‑148 (MQ=255)
                                                                                      atatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCaa                                                          <  2:3202772/148‑1 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGACTGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         >  2:3367757/1‑148 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         >  2:1467189/1‑148 (MQ=255)
                                                                                       tatCGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAt                                                         <  1:2643487/148‑1 (MQ=255)
                                                                                          cGCCTCTTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAAtttt                                                      >  2:3525731/1‑148 (MQ=255)
                                                                                             ctctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGc                                                   >  1:738798/1‑148 (MQ=255)
                                                                                               ctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAg                                                 >  1:1176152/1‑148 (MQ=255)
                                                                                               ctTTTATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAg                                                 >  2:1669047/1‑148 (MQ=255)
                                                                                                 tttATTGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGaa                                               <  1:2470523/148‑1 (MQ=255)
                                                                                                     ttGTGGGATTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTg                                           >  1:3579616/1‑148 (MQ=255)
                                                                                                            aTTGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAag                                    <  1:3756345/148‑1 (MQ=255)
                                                                                                             ttGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGAGCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAaga                                   >  1:426947/1‑148 (MQ=255)
                                                                                                              tGACCGTGCTAAAGTGAATGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGATTGAAGTTATTGAACTCCTCCAGGCAGCATGGCAAACAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAagag                                  <  1:3654648/148‑1 (MQ=255)
                                                                                                              tGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAagag                                  <  1:32568/148‑1 (MQ=255)
                                                                                                              tGACCCTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAagag                                  >  2:1015789/1‑148 (MQ=255)
                                                                                                                acacTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTc                                <  2:2024565/145‑1 (MQ=255)
                                                                                                                 cccTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCa                               >  2:701281/1‑148 (MQ=255)
                                                                                                                   cTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAgg                             >  2:578418/1‑148 (MQ=255)
                                                                                                                   cTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAgg                             <  2:527640/148‑1 (MQ=255)
                                                                                                                   cTGCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAgg                             >  2:1536050/1‑148 (MQ=255)
                                                                                                                     gCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGtt                           <  2:1721890/148‑1 (MQ=255)
                                                                                                                     gCTAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGtt                           >  1:3142975/1‑148 (MQ=255)
                                                                                                                       tAACGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGtttt                         >  1:2802337/1‑148 (MQ=255)
                                                                                                                        aaCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTg                        >  2:2533173/1‑148 (MQ=255)
                                                                                                                        aaCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTg                        >  2:2691995/1‑148 (MQ=255)
                                                                                                                        aaCGTGTCTGTCTCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTg                        <  2:1125407/148‑1 (MQ=255)
                                                                                                                             gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGc                   <  1:3000493/148‑1 (MQ=255)
                                                                                                                             gtctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGc                   <  2:1562804/148‑1 (MQ=255)
                                                                                                                               ctgtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGa                 >  1:2521960/1‑148 (MQ=255)
                                                                                                                                 gtctCAAGAGTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAAc               <  1:1270593/148‑1 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        >  1:1931235/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        <  2:3520496/148‑1 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        >  1:2446944/1‑148 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        <  1:2638346/148‑1 (MQ=255)
                                                                                                                                        gagTAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGa        >  1:2447512/1‑148 (MQ=255)
                                                                                                                                           tAAGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGAttt     >  1:3117391/1‑148 (MQ=255)
                                                                                                                                            aaGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTg    <  2:3699013/148‑1 (MQ=255)
                                                                                                                                            aaGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTg    >  2:1815042/1‑148 (MQ=255)
                                                                                                                                             aGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGa   <  1:2047193/148‑1 (MQ=255)
                                                                                                                                             aGGATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGa   >  2:1641343/1‑148 (MQ=255)
                                                                                                                                              ggATCATCACTATGAAATGTAAACGTCTGAATGAAGTTATTGAACTCCTCCAGCCAGCCTGGCAAAAAGAGCCAGACCTTAACCTGCTGCAATTTTTGCAGAAACTGGCGAAAGAGTCAGGTTTTGACGGCGAACTGGCGGATTTGAc  >  2:2154658/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: