Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 300,770 C→A 86.3% V49V (GTG→GTT ecpE ← ECP production pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273300,7700CA86.3% 398.2 / 69.2 175V49V (GTG→GTTecpEECP production pilus chaperone
Reads supporting (aligned to +/- strand):  ref base C (16/8);  new base A (85/66);  total (101/74)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 175 total aligned reads displayed.

ACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACA  >  CP009273/300625‑300915
                                                                                                                                                 |                                                                                                                                                 
aCTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAAc                                                                                                                                                 <  2:1767171/148‑1 (MQ=255)
aCTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAAc                                                                                                                                                 <  2:592510/148‑1 (MQ=255)
 cTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACt                                                                                                                                                <  2:3256528/148‑1 (MQ=255)
 cTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACt                                                                                                                                                <  2:2962601/148‑1 (MQ=255)
  tCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTa                                                                                                                                               <  2:845119/148‑1 (MQ=255)
  tCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTGTTGTTGTTATTAACTa                                                                                                                                               >  2:2710600/1‑148 (MQ=255)
   cGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTac                                                                                                                                              <  2:2907651/148‑1 (MQ=255)
    gCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTaca                                                                                                                                             >  2:2869048/1‑148 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           >  2:2107198/1‑148 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           >  2:3035546/1‑148 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           >  2:2107301/1‑148 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           >  2:1851982/1‑148 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           <  2:1342/148‑1 (MQ=255)
      tctcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACg                                                                                                                                           >  2:1066011/1‑148 (MQ=255)
        tcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGtt                                                                                                                                         <  2:1816436/148‑1 (MQ=255)
        tcACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGtt                                                                                                                                         >  1:2219583/1‑148 (MQ=255)
         cACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGttt                                                                                                                                        >  1:334504/1‑148 (MQ=255)
          aCCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTg                                                                                                                                       <  1:111810/148‑1 (MQ=255)
           ccAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTAAACGTTTGc                                                                                                                                      >  2:2284668/1‑148 (MQ=255)
              gCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGa                                                                                                                                   >  1:346528/1‑148 (MQ=255)
               ccTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGAc                                                                                                                                  >  1:708033/1‑148 (MQ=255)
                  gCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGNTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACaaa                                                                                                                               <  1:1681549/148‑1 (MQ=255)
                   cAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAag                                                                                                                              <  1:371771/148‑1 (MQ=255)
                   cAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAag                                                                                                                              >  2:2564899/1‑148 (MQ=255)
                   cAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAag                                                                                                                              >  2:826759/1‑148 (MQ=255)
                     aCGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtc                                                                                                                            <  2:3093715/148‑1 (MQ=255)
                     aCGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtc                                                                                                                            <  2:1098533/148‑1 (MQ=255)
                     aCGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtc                                                                                                                            >  2:706924/1‑148 (MQ=255)
                        ccAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtcagt                                                                                                                         >  1:1502975/1‑148 (MQ=255)
                        ccAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtcagt                                                                                                                         >  1:1888876/1‑148 (MQ=255)
                        ccAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtcagt                                                                                                                         <  1:770204/148‑1 (MQ=255)
                         cAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAagtcagtc                                                                                                                        <  2:175562/148‑1 (MQ=255)
                          aGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGtct                                                                                                                       >  2:553633/1‑148 (MQ=255)
                             tGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCgg                                                                                                                    <  2:474861/148‑1 (MQ=255)
                                 ggggggCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACgg                                                                                                                <  1:1899487/148‑1 (MQ=255)
                                 ggggggCGAAAAGCAGTCCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACgg                                                                                                                <  2:522513/148‑1 (MQ=255)
                                    gggCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAg                                                                                                             >  1:1607937/1‑148 (MQ=255)
                                     ggCGAAACGCAGGTCACCATCCACCGGTCGGGTGCGCAAGTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGc                                                                                                            <  2:880513/148‑1 (MQ=255)
                                     ggCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGc                                                                                                            >  2:3026866/1‑148 (MQ=255)
                                         aaaaGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGaaa                                                                                                        <  1:2686357/148‑1 (MQ=255)
                                          aaaGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGaaaa                                                                                                       >  1:3103685/1‑148 (MQ=255)
                                           aaGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAAt                                                                                                      >  1:1032289/1‑148 (MQ=255)
                                            aGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATg                                                                                                     <  1:497159/148‑1 (MQ=255)
                                            aGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATg                                                                                                     >  2:430796/1‑148 (MQ=255)
                                              cAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTc                                                                                                   >  1:1134128/1‑148 (MQ=255)
                                              cAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTc                                                                                                   >  1:2786417/1‑148 (MQ=255)
                                                  tCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGAt                                                                                               >  1:2621343/1‑148 (MQ=255)
                                                    aCCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTg                                                                                             >  1:2403321/1‑148 (MQ=255)
                                                    aCCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTg                                                                                             <  1:800440/148‑1 (MQ=255)
                                                           aCCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGccgaccg                                                                                      >  2:1650615/1‑148 (MQ=255)
                                                             cGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaa                                                                                    <  1:2721565/148‑1 (MQ=255)
                                                              ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaaa                                                                                   <  2:182232/148‑1 (MQ=255)
                                                              ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaaa                                                                                   <  2:1197172/148‑1 (MQ=255)
                                                              ggTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGaaa                                                                                   <  2:1196778/148‑1 (MQ=255)
                                                                tCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                 >  2:202342/1‑148 (MQ=255)
                                                                tCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                 >  1:892853/1‑148 (MQ=255)
                                                                tCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATc                                                                                 >  1:2548140/1‑148 (MQ=255)
                                                                    gTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCt                                                                             <  1:1609974/148‑1 (MQ=255)
                                                                            ttCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCaaa                                                                     >  1:1728322/1‑148 (MQ=255)
                                                                            ttCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCaaa                                                                     >  2:1187219/1‑148 (MQ=255)
                                                                            ttCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCaaa                                                                     <  2:334504/148‑1 (MQ=255)
                                                                               aCTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGt                                                                  >  2:3149367/1‑148 (MQ=255)
                                                                                ctgctgCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTg                                                                 <  1:2979905/148‑1 (MQ=255)
                                                                                          ggCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCAt                                                       >  1:2672301/1‑148 (MQ=255)
                                                                                          ggCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCAt                                                       >  1:1839680/1‑148 (MQ=255)
                                                                                              aTCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAg                                                   >  1:1239455/1‑148 (MQ=255)
                                                                                               tCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGa                                                  >  1:2434658/1‑148 (MQ=255)
                                                                                               tCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGa                                                  <  2:1255335/148‑1 (MQ=255)
                                                                                                cAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGAc                                                 >  1:13575/1‑148 (MQ=255)
                                                                                                cAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGAc                                                 >  1:331072/1‑148 (MQ=255)
                                                                                                cAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGAc                                                 >  1:878242/1‑148 (MQ=255)
                                                                                                    aGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCt                                             <  1:2124808/148‑1 (MQ=255)
                                                                                                        cggATGGCAATACGGTAAATCCGAGCGCTTAGGTTGTTATTCACTCCACGTTTGCTGACAAAGTCAGTCTCGGACGGCCGCGAAAATGACAGATTGCCTACCGAAATCGCCAGTGTCAAAGCTGCAGCCAGCATACAGACCGATGTca                                         >  2:723302/3‑148 (MQ=255)
                                                                                                          gATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTcagc                                       >  2:797649/1‑148 (MQ=255)
                                                                                                           aTGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTcagca                                      <  2:1332915/148‑1 (MQ=255)
                                                                                                            tGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTcagcag                                     >  1:3229588/1‑148 (MQ=255)
                                                                                                               cAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGc                                  >  1:3002286/1‑148 (MQ=255)
                                                                                                               cAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGc                                  >  1:2400694/1‑148 (MQ=255)
                                                                                                               cAATACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGc                                  >  1:2400684/1‑148 (MQ=255)
                                                                                                                 aTACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCCGGGccg                                >  2:1710448/1‑147 (MQ=255)
                                                                                                                 aTACGGTATATCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCt                                <  1:611995/148‑1 (MQ=255)
                                                                                                                       tataTCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTa                          >  2:2487491/1‑148 (MQ=255)
                                                                                                                        atatCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAAGCGCCTGTGACAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTGAGTTaa                         >  2:2170586/1‑148 (MQ=255)
                                                                                                                        atatCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTaa                         >  2:3188175/1‑148 (MQ=255)
                                                                                                                          atCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATg                       >  2:3222708/1‑148 (MQ=255)
                                                                                                                           tCCGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGt                      >  2:2039766/1‑148 (MQ=255)
                                                                                                                            ccGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                     >  1:2032571/1‑148 (MQ=255)
                                                                                                                            ccGTGCGCTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGtt                     >  2:2733132/1‑148 (MQ=255)
                                                                                                                                gcgcTTTTGTTGTTATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGc                 >  2:403532/1‑148 (MQ=255)
                                                                                                                                       tgttgtTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCg          >  2:1844039/1‑148 (MQ=255)
                                                                                                                                         ttgttATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGcc        <  2:1518604/148‑1 (MQ=255)
                                                                                                                                         ttgttATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGcc        <  2:1518588/148‑1 (MQ=255)
                                                                                                                                          tgttATTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCt       >  1:487844/1‑148 (MQ=255)
                                                                                                                                            ttATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGa     >  1:2803864/1‑148 (MQ=255)
                                                                                                                                              aTTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAAc   >  2:3153358/1‑148 (MQ=255)
                                                                                                                                              aTTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAAc   <  1:31340/148‑1 (MQ=255)
                                                                                                                                              aTTAACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAAc   >  1:3172200/1‑148 (MQ=255)
                                                                                                                                               ttCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACa  <  1:2874678/148‑1 (MQ=255)
                                                                                                                                               ttCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACa  <  1:2874649/148‑1 (MQ=255)
                                                                                                                                               ttCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACa  >  1:2551733/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ACTCGCTCTCACCAGCCTGCAACGCCAGCTGGCGGGGGGCGAAAAGCAGTTCACCATCCACCGGTCGGGTGCGCAATTCACTGCTGCCCGGGCTATCAATAGCACTGATGGCAATACGGTATATCCGTGCGCTTTTGTTGTTATTCACTACACGTTTGCTGACAAAGTCAGTCTCGGACGGCAGCGAAAATGTCAGATTGCCGACCGAAATCGCCTGTGTCAAAGGTGCCGCCAGCATACAGACCGCTGTCAGCAGGGCCTTAGTTAATGTTACGCCACGTCGCCTGAACA  >  CP009273/300625‑300915

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: