Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 403,182 T→A 9.6% intergenic (+20/‑99) phoA → / → psiF bacterial alkaline phosphatase/PsiF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913403,1820TA 9.6% 296.6 / 18.0 114intergenic (+20/‑99)phoA/psiFbacterial alkaline phosphatase/PsiF family protein
Reads supporting (aligned to +/- strand):  ref base T (59/42);  new base A (7/4);  total (68/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 117 total aligned reads displayed.

GGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403033‑403327
                                                                                                                                                     |                                                                                                                                                 
   ttCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAtttt                                                                                                                                                <  2:748927/150‑1 (MQ=255)
    ccacaAGAACATCCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGATGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTc                                                                                                                                               <  2:493777/149‑1 (MQ=255)
     cacaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCg                                                                                                                                              <  2:1317560/150‑1 (MQ=255)
      acaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAACCCGCGCCCGGCATTGAATTTTCgg                                                                                                                                             >  1:988681/1‑149 (MQ=255)
      acaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGc                                                                                                                                             <  1:1868/150‑1 (MQ=255)
            aCATACCGGCAGGCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCACTGCCgg                                                                                                                                       <  1:921013/150‑1 (MQ=255)
             cATCCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCATCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCCCTGCCggg                                                                                                                                      <  2:619545/150‑1 (MQ=255)
             cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCggg                                                                                                                                      >  2:956556/1‑150 (MQ=255)
             cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCACTGCCggg                                                                                                                                      <  1:1007320/150‑1 (MQ=255)
                 ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGggtggt                                                                                                                                  >  2:397896/1‑150 (MQ=255)
                 ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCAGggtggt                                                                                                                                  >  2:164661/1‑150 (MQ=255)
                   ggCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTCCACCATGCAAGCCGCTCTGGGGCTGAAATAAATCCGCGCCCGGCAGTGAAGTTTCGCTGCCGGGTGGttt                                                                                                                                >  1:73935/1‑150 (MQ=255)
                      agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGtttttt                                                                                                                             <  1:512390/150‑1 (MQ=255)
                      agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGtttttt                                                                                                                             <  1:998115/150‑1 (MQ=255)
                        tcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCGAATTTTTTTTTTATGTATTAGGAGACCTATGTGTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTg                                                                                                                           <  2:584538/150‑1 (MQ=255)
                          cgtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCt                                                                                                                         <  1:794757/149‑1 (MQ=255)
                              gCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTa                                                                                                                     <  1:868015/150‑1 (MQ=255)
                                 tATTGCGTCGTATGGCCCGCATTCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCa                                                                                                                  <  2:321821/150‑1 (MQ=255)
                                 tATTGCGGCGTATGGCCCACATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCa                                                                                                                  <  1:1184397/150‑1 (MQ=255)
                                  aTTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCTTTAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCaa                                                                                                                 <  2:231473/150‑1 (MQ=255)
                                     gcggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCa                                                                                                              >  2:1126743/1‑150 (MQ=255)
                                         cgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACt                                                                                                          <  1:1069117/150‑1 (MQ=255)
                                          gTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACtt                                                                                                         <  2:1447276/150‑1 (MQ=255)
                                          gTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACtt                                                                                                         <  1:1039962/150‑1 (MQ=255)
                                           tATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGTCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa                                                                                                        <  1:1747054/150‑1 (MQ=255)
                                           tATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTa                                                                                                        <  1:861269/150‑1 (MQ=255)
                                            aTGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTaa                                                                                                       >  1:201616/1‑150 (MQ=255)
                                               gCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATgg                                                                                                    <  1:1609336/150‑1 (MQ=255)
                                                 ccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCa                                                                                                  >  1:1253662/1‑150 (MQ=255)
                                                  cGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAg                                                                                                 <  1:1323420/150‑1 (MQ=255)
                                                     aTGCCGCCAATGGTGTTGGACGTACCGACAAGACCGATCTCTTCTACACCATTAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATc                                                                                              <  2:92539/150‑1 (MQ=255)
                                                      tGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCa                                                                                             >  2:123524/1‑150 (MQ=255)
                                                        ccgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACg                                                                                           >  2:1789167/1‑150 (MQ=255)
                                                         cgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCCGCAGTGAATTTTATCTGCCGGGTGGTTTTTTTTCTGTTAGCAAACAGACTAAATTGCAGATCACgg                                                                                          >  2:850739/1‑150 (MQ=255)
                                                          gccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACggg                                                                                         >  2:1501665/1‑150 (MQ=255)
                                                            cAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAACGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcg                                                                                       >  1:1707820/1‑150 (MQ=255)
                                                             aaTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcgc                                                                                      >  1:1066491/1‑150 (MQ=255)
                                                             aaTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACAACATGAAAGCCGCTCTGGGACTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcgc                                                                                      <  1:362097/150‑1 (MQ=255)
                                                                 ttgttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAc                                                                                  >  2:596926/1‑150 (MQ=255)
                                                                 ttgttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAc                                                                                  >  2:162062/1‑150 (MQ=255)
                                                                   gttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                >  2:1087774/1‑150 (MQ=255)
                                                                   gtaggAATGACCGACCAGACCGATATCTTATACACCATGAAAGCAGCTATGGTGCAGAAATAAAACCGCGCCCTGCAGAGAATTTTCGATGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACATAAATGCAGATAACCGCCGCATACGc                                                                                >  2:1813710/4‑150 (MQ=255)
                                                                    ttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCt                                                                               <  2:1601478/150‑1 (MQ=255)
                                                                        aCTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATg                                                                           <  2:931381/150‑1 (MQ=255)
                                                                           gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGtt                                                                        >  1:1244774/1‑150 (MQ=255)
                                                                           gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGtt                                                                        <  1:580007/150‑1 (MQ=255)
                                                                           gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAAGCTCATGGtt                                                                        >  1:1552721/1‑150 (MQ=255)
                                                                            accgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTa                                                                       <  1:1169506/150‑1 (MQ=255)
                                                                               gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:318084/1‑150 (MQ=255)
                                                                               gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:970368/1‑150 (MQ=255)
                                                                               gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGGTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:238941/1‑150 (MQ=255)
                                                                                 ccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACa                                                                  >  1:162003/1‑150 (MQ=255)
                                                                                   aGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCAGTGTTAAAACATg                                                                >  2:1095955/1‑150 (MQ=255)
                                                                                   aGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACGTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATg                                                                >  1:255002/1‑150 (MQ=255)
                                                                                    gACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGa                                                               >  1:1429187/1‑150 (MQ=255)
                                                                                      ccGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAag                                                             >  1:485350/1‑150 (MQ=255)
                                                                                         aTCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAgg                                                          <  2:592140/150‑1 (MQ=255)
                                                                                                ctACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                   >  2:942175/1‑150 (MQ=255)
                                                                                                 tACACCATTAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATTACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGc                                                  <  2:633097/150‑1 (MQ=255)
                                                                                                 tACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGc                                                  <  2:1253662/150‑1 (MQ=255)
                                                                                                   cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                >  1:577877/1‑150 (MQ=255)
                                                                                                    cccttGAACGCCGCGATGGGGCTGAAATAAAACCGCGCCCGGCAGGGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGAATTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               <  2:1540012/146‑1 (MQ=255)
                                                                                                    acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               <  1:1069069/150‑1 (MQ=255)
                                                                                                    acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               >  2:4536/1‑150 (MQ=255)
                                                                                                     ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATg                                              >  1:428859/1‑150 (MQ=255)
                                                                                                     ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATg                                              >  1:1519488/1‑150 (MQ=255)
                                                                                                     ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAGTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTACAACATGAAGAGGGATGGTGCTagg                                              >  1:421983/1‑148 (MQ=255)
                                                                                                           aaGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATa                                        >  2:1335510/1‑150 (MQ=255)
                                                                                                            aGCCGCTCTGGGGCAGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTCCCGGGTGGTTATTTTGCTCTTAGCAACCAGACTAAATGGCAGATCACGGGCGAATACGCTCAGGGTAAAAACATGAAGAGGGATGGTGCTATGAaaaaaa                                       >  1:362079/1‑150 (MQ=255)
                                                                                                             gCCGCTCTGGGGCTGAAATAAATCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGATTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGATGAGGGATGGTGATATGAAAATAAc                                      >  1:244321/1‑150 (MQ=255)
                                                                                                              ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                     >  2:681422/1‑150 (MQ=255)
                                                                                                              ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                     >  1:1213755/1‑150 (MQ=255)
                                                                                                              ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGATTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGAATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                     >  1:870255/1‑150 (MQ=255)
                                                                                                                gCTCTGGGGCTGAAATAAAACCGCTCCCTGCAGTGAATTTTCGCTGCCGTGTGGTTTTTTTGCTGTTAGCAACCAGTCTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACAtt                                   <  1:1612775/150‑1 (MQ=255)
                                                                                                                gCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACAtt                                   >  2:1225335/1‑150 (MQ=255)
                                                                                                                   ctGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATAACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACt                                >  2:1518130/1‑150 (MQ=255)
                                                                                                                     ggggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTgg                              >  1:1066865/1‑150 (MQ=255)
                                                                                                                      gggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGt                             <  1:426379/150‑1 (MQ=255)
                                                                                                                       ggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGtt                            >  2:1025/1‑150 (MQ=255)
                                                                                                                          tGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTAcc                         <  1:436300/150‑1 (MQ=255)
                                                                                                                            aaaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATAACGGGCGCATAAGCTCATGGTTAAAACATGAAGAGCGATGGTGATATGAAAATAAAATTACTGGTTACCtt                       >  2:1743573/1‑150 (MQ=255)
                                                                                                                             aaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTCCCTTg                      >  2:1391381/1‑150 (MQ=255)
                                                                                                                             aaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTg                      >  1:1332351/1‑150 (MQ=255)
                                                                                                                              aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGATGGATGGTGCTATGAAAATAAGATTACTGGTTACATTGc                     >  2:338314/1‑150 (MQ=255)
                                                                                                                              aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGc                     >  1:1705548/1‑150 (MQ=255)
                                                                                                                              aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGc                     >  2:818401/1‑150 (MQ=255)
                                                                                                                               tAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCt                    >  1:659070/1‑150 (MQ=255)
                                                                                                                               tAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGAGGGTGCTATGAAAATAACATTACTGGTTACCTTgat                    >  1:1378098/1‑148 (MQ=255)
                                                                                                                               tAAAACCGCGACAGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTTTTAGCAACAAGACTTAATGGCAGATCACTTGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCt                    <  1:1185664/150‑1 (MQ=255)
                                                                                                                                      ccacccGGCAGCGATAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGAGGGATGGTGCTATGATAATAACATTAGTGGTTACCTTGCTTTTCGGt             >  1:82330/4‑150 (MQ=255)
                                                                                                                                       cgcCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATAGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTc            <  2:433002/150‑1 (MQ=255)
                                                                                                                                       cacccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTACCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATCAAGAGGGATGGTCCTATGAAAATAACATTACTGGTTACATAGCTTTTCggac            >  2:596817/3‑148 (MQ=255)
                                                                                                                                        acccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt           >  2:1111/2‑150 (MQ=255)
                                                                                                                                        acccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt           >  1:56324/2‑150 (MQ=255)
                                                                                                                                         cccGGCAGTGAATTTTCGCTGCCGGGGGGGTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGAGCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTg          <  1:98250/150‑1 (MQ=255)
                                                                                                                                         cccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTg          >  2:128880/1‑150 (MQ=255)
                                                                                                                                           cGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCATCCAGACTTAATGGCAGCTCACGGGCGCATACGCTGATGGTTAAAACATGAAGAGGGATGGTGCGATGAAAATAACATTACTGGTTACCTGGCTTTTCGTTCTGGt        >  2:79753/1‑150 (MQ=255)
                                                                                                                                             gCAGTGAATTTTCGCTGCCCGGTGTTGTTTTTGCTGATAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGAGGGATGGTGCTATGAAAATAACATTAATGGTTACCTTGATTTTCGGTCTGGttt      >  2:975497/1‑150 (MQ=255)
                                                                                                                                                gTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttttt   <  2:177940/150‑1 (MQ=255)
                                                                                                                                                 tGAATTTGCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACCCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt  <  2:73935/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403033‑403327

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: