Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 403,183 T→A 9.4% intergenic (+21/‑98) 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,1830TA 9.4% 316.0 / 17.2 114intergenic (+21/‑98)phoA/psiFbacterial alkaline phosphatase/PsiF family protein
Reads supporting (aligned to +/- strand):  ref base T (61/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.88e-01

Only 100 of 115 total aligned reads displayed.

TTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403036‑403327
                                                                                                                                                   |                                                                                                                                                
ttCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAtttt                                                                                                                                                <  2:748927/150‑1 (MQ=255)
 ccacaAGAACATCCCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGATGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTc                                                                                                                                               <  2:493777/149‑1 (MQ=255)
  cacaAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCg                                                                                                                                              >  1:728631/1‑150 (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)
          cATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCggg                                                                                                                                      >  2:924574/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)
              ccGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCAGggtggt                                                                                                                                  >  2:164661/1‑150 (MQ=255)
                ggCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTCCACCATGCAAGCCGCTCTGGGGCTGAAATAAATCCGCGCCCGGCAGTGAAGTTTCGCTGCCGGGTGGttt                                                                                                                                >  1:73935/1‑150 (MQ=255)
                   agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGtttttt                                                                                                                             <  1:998115/150‑1 (MQ=255)
                   agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGtttttt                                                                                                                             <  1:512390/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)
                                        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)
                                                gCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAATCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGa                                                                                                <  1:1280636/150‑1 (MQ=255)
                                                  aTGCCGCCAATGGTGTTGGACGTACCGACAAGACCGATCTCTTCTACACCATTAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATc                                                                                              <  2:92539/150‑1 (MQ=255)
                                                     ccgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACg                                                                                           >  2:1789167/1‑150 (MQ=255)
                                                      cgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACgg                                                                                          >  1:446219/1‑150 (MQ=255)
                                                      cgccAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCCGCAGTGAATTTTATCTGCCGGGTGGTTTTTTTTCTGTTAGCAAACAGACTAAATTGCAGATCACgg                                                                                          >  2:850739/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)
                                                                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)
                                                                        gaccgaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATAAGCTCATGGtt                                                                        >  1:1552721/1‑150 (MQ=255)
                                                                            gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:970368/1‑150 (MQ=255)
                                                                            gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:318084/1‑150 (MQ=255)
                                                                            gaccAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGGTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTaaaa                                                                    >  1:238941/1‑150 (MQ=255)
                                                                             acaaGACAGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCACTAAATTTTCGATGCCGGGAGGTTTTTTTGCTGTTAGCAACCAGACTTAATTGCAGATCAAGGGCGCATACGCTCATGGTTAAAAc                                                                   >  2:63729/4‑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)
                                                                                          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)
                                                                                                cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                >  1:577877/1‑150 (MQ=255)
                                                                                                 cccttGAACGCCGCGATGGGGCTGAAATAAAACCGCGCCCGGCAGGGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGAATTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               <  2:1540012/146‑1 (MQ=255)
                                                                                                 acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               >  2:1624113/1‑150 (MQ=255)
                                                                                                 acCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTAt                                               >  2:4536/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)
                                                                                                           ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                     >  2:681422/1‑150 (MQ=255)
                                                                                                           ccGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACa                                     >  1:1213755/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)
                                                                                                                       tGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTAcc                         <  1:436300/150‑1 (MQ=255)
                                                                                                                         aaaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATAACGGGCGCATAAGCTCATGGTTAAAACATGAAGAGCGATGGTGATATGAAAATAAAATTACTGGTTACCtt                       >  2:1743573/1‑150 (MQ=255)
                                                                                                                          aaTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTCCCTTg                      >  2:1391381/1‑150 (MQ=255)
                                                                                                                           aTAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTTAAACATGAAGATGGATGGTGCTATGAAAATAAGATTACTGGTTACATTGc                     >  2:338314/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)
                                                                                                                              aaaCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCttt                  <  1:1126743/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)
                                                                                                                                     acccGGCAGCGAAAATTCACTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt           >  2:1111/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)
                                                                                                                                              tGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt  <  1:1789167/150‑1 (MQ=255)
                                                                                                                                              tGAATTTGCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACCCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttttt  <  2:73935/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTT  >  NC_000913/403036‑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: