Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 403,187 G→A 10.6% intergenic (+25/‑94) phoA → / → psiF alkaline phosphatase/PsiF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913403,1870GA10.6% 267.9 / 15.4 94intergenic (+25/‑94)phoA/psiFalkaline phosphatase/PsiF family protein
Reads supporting (aligned to +/- strand):  ref base G (38/46);  new base A (5/5);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCAC  >  NC_000913/403038‑403333
                                                                                                                                                     |                                                                                                                                                    
cacaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCg                                                                                                                                                      >  1:504091/1‑150 (MQ=255)
  caAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCt                                                                                                                                                    <  1:53928/150‑1 (MQ=255)
    aGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGc                                                                                                                                                  >  1:724704/1‑150 (MQ=255)
     gAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGcc                                                                                                                                                 <  2:509402/150‑1 (MQ=255)
     gAACATACCAGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGcc                                                                                                                                                 >  1:1223053/1‑150 (MQ=255)
      aaCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCACTGCCg                                                                                                                                                <  1:319561/150‑1 (MQ=255)
       aCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCgg                                                                                                                                               >  2:404712/1‑150 (MQ=255)
       aCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCACTGCCgg                                                                                                                                               <  2:997357/150‑1 (MQ=255)
       aCATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAAC‑‑CACCCGGCAGCGAAAATTCACTGCCgg‑‑gc                                                                                                                                           <  1:442853/150‑2 (MQ=255)
        cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCgg‑‑g                                                                                                                                            >  2:311763/1‑150 (MQ=255)
        cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCCCCCGGCAGCGAAAATTCACTGCCgg‑‑g                                                                                                                                            <  1:681450/150‑1 (MQ=255)
           aCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGg‑‑gtgg                                                                                                                                         <  1:190426/150‑1 (MQ=255)
            ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGg‑‑gtggt                                                                                                                                        >  2:91735/1‑150 (MQ=255)
            ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGg‑‑gtggt                                                                                                                                        >  2:358577/1‑150 (MQ=255)
                 agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGtttttt                                                                                                                                   <  1:192069/150‑1 (MQ=255)
                 agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGtttttt                                                                                                                                   <  1:186521/150‑1 (MQ=255)
                  gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                                  >  1:844490/1‑150 (MQ=255)
                  gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTAAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                                  >  2:724005/1‑150 (MQ=255)
                  gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTAAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                                  >  2:718531/1‑150 (MQ=255)
                   tcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTg                                                                                                                                 <  1:349194/150‑1 (MQ=255)
                   tcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTg                                                                                                                                 <  1:124887/150‑1 (MQ=255)
                    cagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGc                                                                                                                                <  1:404712/150‑1 (MQ=255)
                    cagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGc                                                                                                                                >  1:1251253/1‑150 (MQ=255)
                         gCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGtag                                                                                                                           >  1:23470/1‑148 (MQ=255)
                         gCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTa                                                                                                                           <  2:108287/150‑1 (MQ=255)
                         gCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTa                                                                                                                           >  2:966076/1‑150 (MQ=255)
                          cGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAg                                                                                                                          >  2:77757/1‑150 (MQ=255)
                           gTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGc                                                                                                                         <  2:1108092/150‑1 (MQ=255)
                             aTTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCaa                                                                                                                       <  1:378134/150‑1 (MQ=255)
                               tGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAAcc                                                                                                                     <  2:372185/150‑1 (MQ=255)
                                gcggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCa                                                                                                                    >  2:1106519/1‑150 (MQ=255)
                                 cggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAg                                                                                                                   <  2:693334/150‑1 (MQ=255)
                                 cggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAg                                                                                                                   <  2:49907/150‑1 (MQ=255)
                                 cggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAg                                                                                                                   <  1:165811/150‑1 (MQ=255)
                                  ggcgTATGGCCCGCATGCCGCCATTTTTTTTTTTATTTATTAGGAGACCTATGTGTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGa                                                                                                                  <  2:901656/150‑1 (MQ=255)
                                        tGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                            >  2:972944/1‑150 (MQ=255)
                                        tGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                            <  1:1354305/150‑1 (MQ=255)
                                        tGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                            <  1:870834/150‑1 (MQ=255)
                                         ggCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATg                                                                                                           <  2:820834/150‑1 (MQ=255)
                                          gCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATgg                                                                                                          <  2:1132745/150‑1 (MQ=255)
                                           cccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                         <  2:1385768/150‑1 (MQ=255)
                                           cccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                         >  2:30010/1‑150 (MQ=255)
                                                  gccgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCAc                                                                                                  >  1:310829/1‑150 (MQ=255)
                                                        aaTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgcgc                                                                                            <  1:127063/150‑1 (MQ=255)
                                                         aTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCa                                                                                           >  2:1302042/1‑150 (MQ=255)
                                                         aTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCa                                                                                           <  2:1154610/150‑1 (MQ=255)
                                                         aTGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGGTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCa                                                                                           <  2:943629/150‑1 (MQ=255)
                                                            ttgttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAc                                                                                        <  2:377049/150‑1 (MQ=255)
                                                              gttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                      >  2:15003/1‑150 (MQ=255)
                                                              gttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                      <  2:127491/150‑1 (MQ=255)
                                                                     tGACCTACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGt                                                                               <  2:724704/150‑1 (MQ=255)
                                                                           accAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAAc                                                                         <  2:1223053/150‑1 (MQ=255)
                                                                            ccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACa                                                                        >  1:271922/1‑150 (MQ=255)
                                                                                   gATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAagag                                                                 >  1:258798/1‑150 (MQ=255)
                                                                                      ctctTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGa                                                              <  1:120817/150‑1 (MQ=255)
                                                                                              cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCTG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                      >  1:1026095/1‑150 (MQ=255)
                                                                                              cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                      <  1:1106519/150‑1 (MQ=255)
                                                                                              cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                      >  2:1150000/1‑150 (MQ=255)
                                                                                                ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCATATTTAATGGCATATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTTTg                                                    <  2:844490/150‑1 (MQ=255)
                                                                                                ccATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATg                                                    >  1:1202947/1‑150 (MQ=255)
                                                                                                 cATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGa                                                   >  1:1167302/1‑150 (MQ=255)
                                                                                                      aaGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATa                                              <  1:91735/150‑1 (MQ=255)
                                                                                                      aaGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATa                                              <  2:92589/150‑1 (MQ=255)
                                                                                                        gCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAAc                                            >  1:766380/1‑150 (MQ=255)
                                                                                                         ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                           >  1:1093368/1‑150 (MQ=255)
                                                                                                           gCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACAtt                                         <  2:242770/150‑1 (MQ=255)
                                                                                                                ggggCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTgg                                    <  2:310829/150‑1 (MQ=255)
                                                                                                                    cTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTAc                                >  2:440667/1‑150 (MQ=255)
                                                                                                                    cTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTAc                                <  1:966076/150‑1 (MQ=255)
                                                                                                                       aaaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCtt                             >  2:150174/1‑150 (MQ=255)
                                                                                                                       aaaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCtt                             >  2:1367296/1‑150 (MQ=255)
                                                                                                                        aaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTg                            >  1:149542/1‑150 (MQ=255)
                                                                                                                            aaaCCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCttt                        <  1:972944/150‑1 (MQ=255)
                                                                                                                              aCCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTc                      <  1:718531/150‑1 (MQ=255)
                                                                                                                              aCCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTc                      <  1:724005/150‑1 (MQ=255)
                                                                                                                                cgcgCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCgg                    <  2:1167302/150‑1 (MQ=255)
                                                                                                                                 ccacccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCgg                    >  1:1397586/4‑149 (MQ=255)
                                                                                                                                 ccacccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCgg                    <  2:1397586/146‑1 (MQ=255)
                                                                                                                                  cgcCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATTGTGCTATGAAAATAACATTACTGGTTAACTTGCTTTTCGGTc                  <  2:1251253/150‑1 (MQ=255)
                                                                                                                                   acccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt                 >  1:1109855/2‑150 (MQ=255)
                                                                                                                                   acccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt                 >  1:937233/2‑150 (MQ=255)
                                                                                                                                   acccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt                 >  2:568519/2‑150 (MQ=255)
                                                                                                                                     ccGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTgg               <  1:551403/150‑1 (MQ=255)
                                                                                                                                     ccGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGGGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTgg               >  1:748037/1‑150 (MQ=255)
                                                                                                                                     ccGGCAGCGAAAATTCACTGCCGGGCGCGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt                 >  2:88322/1‑150 (MQ=255)
                                                                                                                                      cGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGt              >  2:1473331/1‑150 (MQ=255)
                                                                                                                                          aGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt          >  1:388121/1‑150 (MQ=255)
                                                                                                                                           gTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttttt         <  1:358577/150‑1 (MQ=255)
                                                                                                                                           gTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttttt         <  1:1367296/150‑1 (MQ=255)
                                                                                                                                           gTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGtttttt         >  2:23004/1‑150 (MQ=255)
                                                                                                                                            tGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt        <  1:1231776/150‑1 (MQ=255)
                                                                                                                                              aaTTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTg                            >  1:983231/1‑128 (MQ=255)
                                                                                                                                              aaTTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTg                            <  2:983231/128‑1 (MQ=255)
                                                                                                                                               aTTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAac     >  1:835521/1‑150 (MQ=255)
                                                                                                                                                  ttCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAaccac  >  2:372154/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCAC  >  NC_000913/403038‑403333

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