Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 403,178 T→C 9.6% intergenic (+16/‑103) 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,1780TC 9.6% 307.8 / 20.5 116intergenic (+16/‑103)phoA/psiFbacterial alkaline phosphatase/PsiF family protein
Reads supporting (aligned to +/- strand):  ref base T (63/40);  new base C (7/4);  total (71/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-02

Only 100 of 118 total aligned reads displayed.

GAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403031‑403327
                                                                                                                                                   |                                                                                                                                                     
gagGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGa                                                                                                                                                     >  2:861269/1‑150 (MQ=255)
  ggATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAt                                                                                                                                                   >  1:1591505/1‑150 (MQ=255)
     ttCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAtttt                                                                                                                                                <  2:748927/150‑1 (MQ=255)
      ccacaAGAACATCCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGATGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTc                                                                                                                                               <  2:493777/149‑1 (MQ=255)
       cacaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCg                                                                                                                                              >  1:728631/1‑150 (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)
               cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGCGAAAATTCACTGCCggg                                                                                                                                      <  2:1232915/150‑1 (MQ=255)
                 tACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGggtg                                                                                                                                    >  2:1017561/1‑150 (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)
                           cagtTGCGTCTTGCGGCGTATGGCCCGCATGCCGCCCCTGGTGTTGGACTGACCGAGCAGCCCGAGCGCTTCAACACCAAGACCGCCGCGCTGCGGCTGTATTACATCCGCGCCGGTCAGTGAATTTTCGCTGCCGGGTGGGTTATTTGc                                                                                                                          <  1:345853/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)
                                              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)
                                                           cgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACgg                                                                                          >  1:446219/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:162062/1‑150 (MQ=255)
                                                                   ttgttgGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAc                                                                                  >  2:596926/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:580007/150‑1 (MQ=255)
                                                                             gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGtt                                                                        >  1:1244774/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)
                                                                                  acaaGACAGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCACTAAATTTTCGATGCCGGGAGGTTTTTTTGCTGTTAGCAACCAGACTTAATTGCAGATCAAGGGCGCATACGCTCATGGTTAAAAc                                                                   >  2:63729/4‑150 (MQ=255)
                                                                                   ccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACa                                                                  >  1:162003/1‑150 (MQ=255)
                                                                                     aGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCAGTGTTAAAACATg                                                                >  2:1095955/1‑150 (MQ=255)
                                                                                      gACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGa                                                               >  1:1429187/1‑150 (MQ=255)
                                                                                        ccGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAag                                                             >  1:485350/1‑150 (MQ=255)
                                                                                               cttctACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATg                                                      >  2:494372/1‑150 (MQ=255)
                                                                                                  ctACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                   >  2:942175/1‑150 (MQ=255)
                                                                                                   tACACCATTAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATTACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGc                                                  <  2:633097/150‑1 (MQ=255)
                                                                                                   tACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGc                                                  <  2:1253662/150‑1 (MQ=255)
                                                                                                    acacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTTCTGATAGCAACCAGACTTAATAGCATATCACGGTCGCATACGCTCATGGTTAAAACATGAATAGTGATGGTgat                                                 >  2:31653/1‑148 (MQ=255)
                                                                                                      cccttGAACGCCGCGATGGGGCTGAAATAAAACCGCGCCCGGCAGGGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGAATTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               <  2:1540012/146‑1 (MQ=255)
                                                                                                      acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               >  2:4536/1‑150 (MQ=255)
                                                                                                      acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               >  2:1624113/1‑150 (MQ=255)
                                                                                                      acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               <  1:1069069/150‑1 (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                                     >  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)
                                                                                                                       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)
                                                                                                                                aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGc                     >  2:818401/1‑150 (MQ=255)
                                                                                                                                aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGc                     >  1:1705548/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)
                                                                                                                                       cgcgCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCgg              <  2:1429187/150‑1 (MQ=255)
                                                                                                                                        ccacccGGCAGCGATAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGAGGGATGGTGCTATGATAATAACATTAGTGGTTACCTTGCTTTTCGGt             >  1:82330/4‑150 (MQ=255)
                                                                                                                                        ccacccGGCAGCGAAAATTCACGGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTGAACACATGAAGAGGGATGGTGCTATGAAAATAACATTCCTGGTTACCTTGCTTTTCGGt             >  1:875763/4‑150 (MQ=255)
                                                                                                                                         cgcCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATAGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTc            <  2:433002/150‑1 (MQ=255)
                                                                                                                                         cacccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTACCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATCAAGAGGGATGGTCCTATGAAAATAACATTACTGGTTACATAGCTTTTCggac            >  2:596817/3‑148 (MQ=255)
                                                                                                                                          acccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt           >  1:56324/2‑150 (MQ=255)
                                                                                                                                           cccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTg          >  2:128880/1‑150 (MQ=255)
                                                                                                                                             cGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCATCCAGACTTAATGGCAGCTCACGGGCGCATACGCTGATGGTTAAAACATGAAGAGGGATGGTGCGATGAAAATAACATTACTGGTTACCTGGCTTTTCGTTCTGGt        >  2:79753/1‑150 (MQ=255)
                                                                                                                                               gCAGTGAATTTTCGCTGCCCGGTGTTGTTTTTGCTGATAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGAGGGATGGTGCTATGAAAATAACATTAATGGTTACCTTGATTTTCGGTCTGGttt      >  2:975497/1‑150 (MQ=255)
                                                                                                                                                   tGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt  <  1:1789167/150‑1 (MQ=255)
                                                                                                                                                   tGAATTTGCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACCCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt  <  2:73935/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403031‑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: