Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 572,129 T→G 40.1% pseudogene (153/1040 nt) nmpC ← DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273572,1290TG40.1% 77.4 / inf 250pseudogene (153/1040 nt)nmpCDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage
Reads supporting (aligned to +/- strand):  ref base T (60/90);  new base G (45/55);  total (105/145)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 256 total aligned reads displayed.

GCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTT  >  CP009273/571980‑572278
                                                                                                                                                     |                                                                                                                                                     
gCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCt                                                                                                                                                       <  1:780580/150‑1 (MQ=255)
gCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCt                                                                                                                                                       <  2:4548296/150‑1 (MQ=255)
 cAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTg                                                                                                                                                      >  2:4061508/1‑150 (MQ=255)
  aaGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGc                                                                                                                                                     >  1:1405195/1‑150 (MQ=255)
  aaGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGc                                                                                                                                                     >  1:3165664/1‑150 (MQ=255)
         gTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTa                                                                                                                                              <  2:4013751/150‑1 (MQ=255)
           tttATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTaga                                                                                                                                            <  2:2388773/150‑1 (MQ=255)
             tATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAgagg                                                                                                                                          <  1:2869846/150‑1 (MQ=255)
             tATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAgagg                                                                                                                                          <  1:2657538/150‑1 (MQ=255)
                cTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGgaga                                                                                                                                       >  2:2627916/1‑150 (MQ=255)
                cTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGgaga                                                                                                                                       <  2:3385218/150‑1 (MQ=255)
                 tttGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGaa                                                                                                                                      <  1:1404338/150‑1 (MQ=255)
                    ggaggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAagta                                                                                                                                   <  1:2929083/150‑1 (MQ=255)
                     gaggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAagtag                                                                                                                                  <  1:68899/150‑1 (MQ=255)
                       ggaACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTg                                                                                                                                <  1:3807631/150‑1 (MQ=255)
                          aCCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCtt                                                                                                                             >  1:119291/1‑150 (MQ=255)
                             ttGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGc                                                                                                                          >  2:3705342/1‑150 (MQ=255)
                             ttGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGc                                                                                                                          >  1:3501284/1‑150 (MQ=255)
                              tGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCa                                                                                                                         <  2:2606743/150‑1 (MQ=255)
                                agaTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCAtt                                                                                                                       <  2:2318748/150‑1 (MQ=255)
                                 gaTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTa                                                                                                                      <  1:35838/150‑1 (MQ=255)
                                  aTTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTaa                                                                                                                     <  1:1020876/150‑1 (MQ=255)
                                    tCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACt                                                                                                                   <  2:4280611/150‑1 (MQ=255)
                                     cAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACtt                                                                                                                  <  1:3669311/150‑1 (MQ=255)
                                      aGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACttt                                                                                                                 <  2:3032860/150‑1 (MQ=255)
                                      aGCGCGGTTGCCTTTGAATTCATAGTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACttt                                                                                                                 <  2:930471/150‑1 (MQ=255)
                                         gcgGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTccc                                                                                                              <  1:485782/150‑1 (MQ=255)
                                          cgGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCg                                                                                                             <  1:1628608/150‑1 (MQ=255)
                                           ggTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGAATCATTAGAGGAGAAGTAGTGCTTAGCATTAACATTCCCGt                                                                                                            >  2:1758177/1‑150 (MQ=255)
                                            gatGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTTGGTTTCACCATAGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTATCCCGTa                                                                                                           >  2:4209690/3‑150 (MQ=255)
                                             ttGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTAc                                                                                                          >  2:976085/1‑150 (MQ=255)
                                                   ttGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATcc                                                                                                    <  2:68170/150‑1 (MQ=255)
                                                   ttGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATcc                                                                                                    >  1:537622/1‑150 (MQ=255)
                                                      aaTTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGc                                                                                                 >  1:910135/1‑150 (MQ=255)
                                                      aaTTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGc                                                                                                 >  1:910103/1‑150 (MQ=255)
                                                        ttCTTATTCCCTCTGTCCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGAATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGAATTAACTTTCCCGTACAGATCCAGCtt                                                                                               <  1:370123/150‑1 (MQ=255)
                                                         tCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTg                                                                                              <  2:3915446/150‑1 (MQ=255)
                                                          cATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGt                                                                                             >  1:4059991/1‑150 (MQ=255)
                                                           atatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGtt                                                                                            <  2:2394755/150‑1 (MQ=255)
                                                           atatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGtt                                                                                            <  1:4002729/150‑1 (MQ=255)
                                                           atatTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGAAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTAGta                                                                                            >  1:2106707/1‑149 (MQ=255)
                                                              ttCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACt                                                                                         <  2:517044/150‑1 (MQ=255)
                                                               tCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTg                                                                                        >  1:4398448/1‑150 (MQ=255)
                                                               acccACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTg                                                                                        >  2:2869049/2‑150 (MQ=255)
                                                                cccACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGt                                                                                       >  2:1513174/1‑150 (MQ=255)
                                                                   aCTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCtt                                                                                    <  2:1559894/150‑1 (MQ=255)
                                                                   aCTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCtt                                                                                    <  1:1146635/150‑1 (MQ=255)
                                                                    cTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCttt                                                                                   <  1:1585165/150‑1 (MQ=255)
                                                                       aCCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTttatt                                                                                >  2:3285622/1‑150 (MQ=255)
                                                                         cGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTATCTTtattat                                                                              >  1:1839521/1‑150 (MQ=255)
                                                                         cGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTtattat                                                                              >  2:1898307/1‑150 (MQ=255)
                                                                          gAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATtata                                                                             <  1:1566770/150‑1 (MQ=255)
                                                                          gAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATtata                                                                             <  1:4063621/150‑1 (MQ=255)
                                                                          gAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACAGTCTTTATtata                                                                             >  1:2269793/1‑150 (MQ=255)
                                                                          gAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATtata                                                                             >  2:125706/1‑150 (MQ=255)
                                                                           aaaCCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATaa                                                                            <  2:1672390/150‑1 (MQ=255)
                                                                            aaCCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATaaa                                                                           <  1:3160219/150‑1 (MQ=255)
                                                                            aaCCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATaaa                                                                           >  2:4363631/1‑150 (MQ=255)
                                                                              ccAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAAtt                                                                         <  1:4352319/150‑1 (MQ=255)
                                                                               cagtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAAttt                                                                        <  2:1303404/150‑1 (MQ=255)
                                                                               cagtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAAttt                                                                        <  2:1946585/150‑1 (MQ=255)
                                                                                agtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTc                                                                       >  2:4562483/1‑150 (MQ=255)
                                                                                 gtcagtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCa                                                                      >  2:2187517/1‑150 (MQ=255)
                                                                                    agtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCt                                                                   <  2:2503535/150‑1 (MQ=255)
                                                                                     gtTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTg                                                                  <  1:1061799/150‑1 (MQ=255)
                                                                                          tCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGa                                                             <  1:3977103/150‑1 (MQ=255)
                                                                                           cGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAg                                                            <  2:1043058/150‑1 (MQ=255)
                                                                                           cGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCAACATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAg                                                            <  2:4148028/150‑1 (MQ=255)
                                                                                            gTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGc                                                           <  2:2117015/150‑1 (MQ=255)
                                                                                             ttGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCa                                                          <  1:1041508/150‑1 (MQ=255)
                                                                                             ttGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCa                                                          <  2:476062/150‑1 (MQ=255)
                                                                                              tGATTTGGGTTTCACCTTTGAAGCCTAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAg                                                         <  1:2983811/150‑1 (MQ=255)
                                                                                                aTTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGAc                                                       >  1:2387038/1‑150 (MQ=255)
                                                                                                  ttGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACAt                                                     <  2:3830718/150‑1 (MQ=255)
                                                                                                  ttGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACAt                                                     >  2:893283/1‑150 (MQ=255)
                                                                                                      gTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGcc                                                 >  1:4516450/1‑150 (MQ=255)
                                                                                                            ccTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGt                                           <  2:3168299/150‑1 (MQ=255)
                                                                                                             cTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTa                                          >  2:364176/1‑150 (MQ=255)
                                                                                                             cTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTa                                          <  2:1955630/150‑1 (MQ=255)
                                                                                                                tGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTg                                       <  2:1591600/150‑1 (MQ=255)
                                                                                                                  aaGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGAt                                     >  1:1736468/1‑150 (MQ=255)
                                                                                                                    gCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGc                                   >  1:4058314/1‑150 (MQ=255)
                                                                                                                      cAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAg                                 <  2:2942752/150‑1 (MQ=255)
                                                                                                                       aaGACGGGCATAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCCTTAGTACTGATGCAGc                                >  1:1954936/1‑150 (MQ=255)
                                                                                                                        aGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCt                               <  1:2492165/150‑1 (MQ=255)
                                                                                                                            gggCATAAGTAGTATCACCTTCATCTGCATCATTAGATTATAATTAGTGCTTAGAATTAACTTTCACGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTTATAAGACATCGACATTAGTACTGATGCAGCTAcag                           <  2:931880/150‑1 (MQ=255)
                                                                                                                              gCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTAcagca                         <  2:730420/150‑1 (MQ=255)
                                                                                                                                 tAAGTAGTATCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGaa                      <  1:752446/150‑1 (MQ=255)
                                                                                                                                   agtagtATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAAt                    >  1:3312777/1‑150 (MQ=255)
                                                                                                                                    gtagtaTCACCATCATCTACATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAAtt                   <  2:234468/150‑1 (MQ=255)
                                                                                                                                         aTCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCAc              >  1:3373500/1‑150 (MQ=255)
                                                                                                                                         aTCACCATCATCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCAc              <  2:3025003/150‑1 (MQ=255)
                                                                                                                                              catcatCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTa         <  1:2024568/150‑1 (MQ=255)
                                                                                                                                              catcatCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTa         <  1:2024574/150‑1 (MQ=255)
                                                                                                                                              catcatCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTa         <  1:1712495/150‑1 (MQ=255)
                                                                                                                                              catcatCGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTa         <  2:1819667/150‑1 (MQ=255)
                                                                                                                                               atcatcTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTaa        >  1:1188151/1‑150 (MQ=255)
                                                                                                                                               atcatcGGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTaa        <  2:575843/150‑1 (MQ=255)
                                                                                                                                                tcatcTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAt       >  2:3765984/1‑150 (MQ=255)
                                                                                                                                                     tGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAAtttttt  <  2:810809/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAAGACGGGTTTTATCTTTGGAGGAACCTTGAGATTCAGCGCGGTTGCCTTTGAATTCATATTCCCACTGACCGAAACCAGTCAGTTGATCGTTGATTTGGGTTTCACCTTTGAAGCCAAGACGGGCATAAGTAGTATCACCATCATCTGCATCATTAGAGGAGAAGTAGTGCTTAGCATTAACTTTCCCGTACAGATCCAGCTTGTTACTGTCTTTATTATAAATTTCAGCTGCCTGAGCAGACATCGCCATTAGTACTGATGCAGCTACAGCAGAAATTGCCACTGTTAATTTTTT  >  CP009273/571980‑572278

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: