Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 403,183 T→A 10.6% intergenic (+21/‑98) 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,1830TA10.6% 270.4 / 19.0 94intergenic (+21/‑98)phoA/psiFalkaline phosphatase/PsiF family protein
Reads supporting (aligned to +/- strand):  ref base T (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 = 6.94e-01

ATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAAC  >  NC_000913/403035‑403330
                                                                                                                                                    |                                                                                                                                                     
aTTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAAttt                                                                                                                                                      <  2:953970/150‑1 (MQ=255)
   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:186521/150‑1 (MQ=255)
                    agtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGtttttt                                                                                                                                <  1:192069/150‑1 (MQ=255)
                     gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                               >  1:844490/1‑150 (MQ=255)
                     gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTAAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                               >  2:718531/1‑150 (MQ=255)
                     gtcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTAAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGttttttt                                                                                                                               >  2:724005/1‑150 (MQ=255)
                      tcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTg                                                                                                                              <  1:124887/150‑1 (MQ=255)
                      tcagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTg                                                                                                                              <  1:349194/150‑1 (MQ=255)
                       cagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGc                                                                                                                             >  1:1251253/1‑150 (MQ=255)
                       cagtTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGc                                                                                                                             <  1:404712/150‑1 (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                                                                                                                <  1:165811/150‑1 (MQ=255)
                                    cggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAg                                                                                                                <  2:693334/150‑1 (MQ=255)
                                    cggcgTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAg                                                                                                                <  2:49907/150‑1 (MQ=255)
                                     ggcgTATGGCCCGCATGCCGCCATTTTTTTTTTTATTTATTAGGAGACCTATGTGTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGa                                                                                                               <  2:901656/150‑1 (MQ=255)
                                           tGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                         <  1:1354305/150‑1 (MQ=255)
                                           tGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAAt                                                                                                         >  2:972944/1‑150 (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:30010/1‑150 (MQ=255)
                                              cccGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGc                                                                                                      <  2:1385768/150‑1 (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                                                   >  2:1150000/1‑150 (MQ=255)
                                                                                                 cacCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTa                                                   <  1:1106519/150‑1 (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                                           <  2:92589/150‑1 (MQ=255)
                                                                                                         aaGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATa                                           <  1:91735/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:937233/2‑150 (MQ=255)
                                                                                                                                      acccGGCAGCGAAAATTCACTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCt              >  1:1109855/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)
                                                                                                                                                    |                                                                                                                                                     
ATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGG‑‑GTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAAC  >  NC_000913/403035‑403330

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: