Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 386,569 C→A 8.6% pseudogene (863/1458 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273386,5690CA 8.6% 868.8 / 40.8 233pseudogene (863/1458 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base C (113/100);  new base A (7/13);  total (120/113)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 235 total aligned reads displayed.

CAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAACGGGGGTT  >  CP009273/386420‑386718
                                                                                                                                                     |                                                                                                                                                     
cAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTc                                                                                                                                                       <  1:481186/150‑1 (MQ=255)
 aaTATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCa                                                                                                                                                      >  2:1902342/1‑150 (MQ=255)
 aaTATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCa                                                                                                                                                      >  2:1902449/1‑150 (MQ=255)
 aaTATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCa                                                                                                                                                      <  2:3089137/150‑1 (MQ=255)
   tatCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAgg                                                                                                                                                    <  1:330128/150‑1 (MQ=255)
   tatCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAgg                                                                                                                                                    <  1:2081734/150‑1 (MQ=255)
   tatCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAgg                                                                                                                                                    >  1:3260690/1‑150 (MQ=255)
    atCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGt                                                                                                                                                   <  1:3897695/150‑1 (MQ=255)
    atCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGt                                                                                                                                                   >  2:4411844/1‑150 (MQ=255)
      ctctGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTc                                                                                                                                                 >  1:2116833/1‑150 (MQ=255)
          gTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGAc                                                                                                                                             >  2:524916/1‑150 (MQ=255)
          gTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATACTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGAc                                                                                                                                             >  1:80495/1‑150 (MQ=255)
            aaaaGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTa                                                                                                                                           <  1:358556/150‑1 (MQ=255)
             aaaGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTAt                                                                                                                                          <  1:3064713/150‑1 (MQ=255)
             aaaGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTAt                                                                                                                                          >  2:3748449/1‑150 (MQ=255)
             aaaGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTAt                                                                                                                                          <  1:268363/150‑1 (MQ=255)
              aaGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATg                                                                                                                                         >  1:1631962/1‑150 (MQ=255)
               aGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATgt                                                                                                                                        >  2:301206/1‑149 (MQ=255)
               aGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGc                                                                                                                                        >  1:1324718/1‑150 (MQ=255)
               aGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGc                                                                                                                                        <  2:1362223/150‑1 (MQ=255)
                gATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCg                                                                                                                                       <  1:4783521/150‑1 (MQ=255)
                 aTTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGa                                                                                                                                      >  1:1600960/1‑150 (MQ=255)
                     cacaGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATgaa                                                                                                                                  <  1:5035/150‑1 (MQ=255)
                                cGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTAt                                                                                                                       >  2:4072123/1‑150 (MQ=255)
                                 gggTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTAtt                                                                                                                      <  2:1745612/150‑1 (MQ=255)
                                     tCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTACAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTc                                                                                                                  >  1:1592085/1‑150 (MQ=255)
                                       aaaTTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAg                                                                                                                >  2:1551238/1‑150 (MQ=255)
                                         aTTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCa                                                                                                              <  1:4534993/150‑1 (MQ=255)
                                           tataTCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAAt                                                                                                            <  1:4100575/150‑1 (MQ=255)
                                            atatCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAAtt                                                                                                           >  1:2315077/1‑150 (MQ=255)
                                             tatCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTc                                                                                                          >  2:3600200/1‑150 (MQ=255)
                                             tatCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTc                                                                                                          >  2:3600235/1‑150 (MQ=255)
                                             tatCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTc                                                                                                          >  2:4757441/1‑150 (MQ=255)
                                              atCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCa                                                                                                         <  2:2480520/150‑1 (MQ=255)
                                              atCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCa                                                                                                         <  2:2480454/150‑1 (MQ=255)
                                              atCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCa                                                                                                         >  2:2576803/1‑150 (MQ=255)
                                                     aTAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGa                                                                                                  <  1:4609234/150‑1 (MQ=255)
                                                     aTAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGa                                                                                                  >  1:638412/1‑150 (MQ=255)
                                                     aTAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGa                                                                                                  >  2:421253/1‑150 (MQ=255)
                                                      tAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGAt                                                                                                 <  1:2796837/150‑1 (MQ=255)
                                                         cAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATggg                                                                                              >  2:1752271/1‑150 (MQ=255)
                                                          aaTACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGt                                                                                             >  2:1828669/1‑150 (MQ=255)
                                                            tACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGa                                                                                           <  1:296414/150‑1 (MQ=255)
                                                            tACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGa                                                                                           >  1:1244498/1‑150 (MQ=255)
                                                            tACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGa                                                                                           >  2:2207226/1‑150 (MQ=255)
                                                             aCTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGAt                                                                                          >  1:2239315/1‑150 (MQ=255)
                                                              cTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATg                                                                                         >  2:2329629/1‑150 (MQ=255)
                                                               ttATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGt                                                                                        <  2:2055257/150‑1 (MQ=255)
                                                               ttATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGt                                                                                        <  1:4719640/150‑1 (MQ=255)
                                                                tATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGtt                                                                                       >  1:2211384/1‑150 (MQ=255)
                                                                  tGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGc                                                                                     >  2:1991670/1‑150 (MQ=255)
                                                                  tGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGc                                                                                     >  2:1308323/1‑150 (MQ=255)
                                                                     aaaaaCTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCtt                                                                                  <  2:946616/149‑1 (MQ=255)
                                                                      aaaaCTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCttt                                                                                 >  1:2734931/1‑150 (MQ=255)
                                                                       aaaCTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTa                                                                                >  1:750388/1‑150 (MQ=255)
                                                                        aaCTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTAc                                                                               >  2:149525/1‑150 (MQ=255)
                                                                          cTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCa                                                                             >  1:3240199/1‑150 (MQ=255)
                                                                            ggTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGc                                                                           >  1:1553361/1‑150 (MQ=255)
                                                                               cACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCAcc                                                                        >  1:276395/1‑150 (MQ=255)
                                                                               cACTTTGGCAAGTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTAGCACTGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCAcc                                                                        <  2:3085770/150‑1 (MQ=255)
                                                                                 cTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTg                                                                      <  2:3884793/150‑1 (MQ=255)
                                                                                    tGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGaa                                                                   <  1:1064436/150‑1 (MQ=255)
                                                                                      gCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATa                                                                 <  2:3260690/150‑1 (MQ=255)
                                                                                       cAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAg                                                                <  1:1718853/150‑1 (MQ=255)
                                                                                           tCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAAt                                                            <  2:3751162/150‑1 (MQ=255)
                                                                                            cTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAAtt                                                           >  1:3699222/1‑150 (MQ=255)
                                                                                            cTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAAtt                                                           <  2:1271257/150‑1 (MQ=255)
                                                                                             tgatgaACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATACGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAAttt                                                          <  1:2674542/150‑1 (MQ=255)
                                                                                              gatgaACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAAtttt                                                         >  2:1114793/1‑150 (MQ=255)
                                                                                                     cGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGa                                                  >  1:948503/1‑150 (MQ=255)
                                                                                                     cGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGa                                                  >  1:1286859/1‑150 (MQ=255)
                                                                                                       aGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAAt                                                <  2:3992553/150‑1 (MQ=255)
                                                                                                        gTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATg                                               <  2:580057/150‑1 (MQ=255)
                                                                                                        gTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATg                                               <  2:2970114/150‑1 (MQ=255)
                                                                                                          gATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGt                                             <  2:2039142/150‑1 (MQ=255)
                                                                                                           aTTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTc                                            >  1:1628185/1‑150 (MQ=255)
                                                                                                            ttATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCa                                           >  2:3335616/1‑150 (MQ=255)
                                                                                                                gCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGAt                                       <  2:3882130/150‑1 (MQ=255)
                                                                                                                    gCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCa                                   >  2:284364/1‑150 (MQ=255)
                                                                                                                        gggTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACgg                               >  2:2105093/1‑150 (MQ=255)
                                                                                                                           tGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGa                            <  2:3224279/150‑1 (MQ=255)
                                                                                                                            gATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGAc                           <  1:3131832/150‑1 (MQ=255)
                                                                                                                            gATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGAc                           <  2:52603/150‑1 (MQ=255)
                                                                                                                             aTGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACg                          <  1:1446158/150‑1 (MQ=255)
                                                                                                                               gTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACggg                        >  2:2003418/1‑150 (MQ=255)
                                                                                                                                 tGCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGt                      <  2:2690039/150‑1 (MQ=255)
                                                                                                                                  gCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTg                     <  2:4538035/150‑1 (MQ=255)
                                                                                                                                  gCAATGTCCTTTACTGCTTAAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTg                     <  1:421253/150‑1 (MQ=255)
                                                                                                                                   cAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGa                    >  2:2795971/1‑150 (MQ=255)
                                                                                                                                   cAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGa                    >  2:797086/1‑150 (MQ=255)
                                                                                                                                    aaTGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGaa                   <  2:2009141/150‑1 (MQ=255)
                                                                                                                                    aaTGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGaa                   <  2:2836686/150‑1 (MQ=255)
                                                                                                                                      tGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAg                 <  2:4648296/150‑1 (MQ=255)
                                                                                                                                           tttACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATAcc            >  2:197777/1‑150 (MQ=255)
                                                                                                                                           tttACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATAcc            >  1:3846921/1‑150 (MQ=255)
                                                                                                                                            ttACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCa           >  2:4344896/1‑150 (MQ=255)
                                                                                                                                              aCTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAAc         <  1:1902449/150‑1 (MQ=255)
                                                                                                                                               cTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAACg        <  1:4411844/150‑1 (MQ=255)
                                                                                                                                                 gCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAACggg      <  1:1124727/150‑1 (MQ=255)
                                                                                                                                                  cTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAACgggg     <  2:2670544/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATGCGATGAAGAACAATGTATTCCTCAGCAATTCAACGCTGATGGGTGATGTTGCCTTTACCAGCACCTGGAATAGTAATTTTGATCCGAATGGTCATGATTCCAACGGTGACGGGGTGAAAGATACCAACGGGGGTT  >  CP009273/386420‑386718

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: