Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,892,363 A→G 6.2% Y291Y (TAT→TAC yoaA ← putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,3630AG 6.2% 243.7 / 5.1 82Y291Y (TAT→TACyoaAputative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base G (3/2);  total (44/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGA  >  NC_000913/1892214‑1892511
                                                                                                                                                     |                                                                                                                                                    
aaCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGa                                                                                                                                                      <  2:598621/150‑1 (MQ=255)
 aCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                                                                                                     <  1:312707/150‑1 (MQ=255)
 aCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                                                                                                     <  1:39275/150‑1 (MQ=255)
     ggCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAAcc                                                                                                                                                 >  1:933376/1‑150 (MQ=255)
       cgcgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                                                                                                               >  1:913778/1‑150 (MQ=255)
       cgcgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                                                                                                               <  2:188274/150‑1 (MQ=255)
         cgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcg                                                                                                                                             <  1:20/150‑3 (MQ=255)
           cTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTc                                                                                                                                           >  2:1148552/1‑150 (MQ=255)
           cTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTc                                                                                                                                           >  2:583941/1‑150 (MQ=255)
             cAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                                                                                                         >  2:1039118/1‑150 (MQ=255)
              aaaTGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgtgagc                                                                                                                                        <  2:831469/150‑8 (MQ=255)
              aaaTGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAcc                                                                                                                                        <  2:1353927/150‑1 (MQ=255)
                    cGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGtt                                                                                                                                   >  1:517809/1‑149 (MQ=255)
                    cGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGtt                                                                                                                                   <  2:517809/149‑1 (MQ=255)
                    cGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTg                                                                                                                                  >  2:199226/1‑150 (MQ=255)
                         ccAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                                                                                                             >  2:1291986/1‑150 (MQ=255)
                         ccAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                                                                                                             >  2:91727/1‑150 (MQ=255)
                          cAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACg                                                                                                                            >  2:1328564/1‑150 (MQ=255)
                           aGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGa                                                                                                                           >  2:1318035/1‑150 (MQ=255)
                             cAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaaa                                                                                                                         >  2:1268983/1‑150 (MQ=255)
                                ggCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATc                                                                                                                      >  1:793382/1‑150 (MQ=255)
                                  cGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCt                                                                                                                    <  1:1039118/150‑1 (MQ=255)
                                   ggAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                                                                                                                   <  2:230177/150‑1 (MQ=255)
                                    gAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCTCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                                                                                  <  2:1211148/150‑1 (MQ=255)
                                         ccccAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCAct                                                                                                             >  2:637238/1‑150 (MQ=255)
                                           ccAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATGGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGGAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActct                                                                                                           <  1:1138829/150‑1 (MQ=255)
                                               tGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGc                                                                                                       <  1:1328564/150‑1 (MQ=255)
                                                    ggttCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGAc                                                                                                  >  2:133827/3‑150 (MQ=255)
                                                    ggttCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGAc                                                                                                  >  2:127882/3‑150 (MQ=255)
                                                    gTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGAc                                                                                                  >  2:1299192/1‑150 (MQ=255)
                                                         gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                             >  1:1632533/1‑150 (MQ=255)
                                                         gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                             >  2:1275775/1‑150 (MQ=255)
                                                         gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                             <  1:1291986/150‑1 (MQ=255)
                                                          ccACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAg                                                                                            <  1:127174/150‑1 (MQ=255)
                                                           cACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgc                                                                                           <  1:1318035/150‑1 (MQ=255)
                                                            aCGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcg                                                                                          <  1:178917/150‑1 (MQ=255)
                                                             cGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAgcgc                                                                                         <  1:1299192/150‑1 (MQ=255)
                                                                cATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCGCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACt                                                                                      <  2:1618909/150‑1 (MQ=255)
                                                                 aTAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGATCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACtt                                                                                     >  1:294380/1‑150 (MQ=255)
                                                                  tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                                                                                    >  1:40629/1‑150 (MQ=255)
                                                                  tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                                                                                    >  1:1438696/1‑150 (MQ=255)
                                                                  tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                                                                                    <  1:199226/150‑1 (MQ=255)
                                                                  tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGAGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                                                                                    >  1:1186446/1‑150 (MQ=255)
                                                                   aaCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACtttt                                                                                   >  1:557796/1‑150 (MQ=255)
                                                                    aCAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTg                                                                                  <  1:1148552/150‑1 (MQ=255)
                                                                     cAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCCCCGAGTTGCAGACGAAAATCCTGCGCACTCTGCGCAAGACGATCAGCGCACTTTTGt                                                                                 >  2:1372449/1‑150 (MQ=255)
                                                                        aGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                                                                              <  1:1287071/150‑1 (MQ=255)
                                                                         gTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAct                                                                             >  1:1364112/1‑150 (MQ=255)
                                                                            ccAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActgct                                                                          >  1:1051242/1‑150 (MQ=255)
                                                                               gggTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTggg                                                                       <  2:64038/150‑1 (MQ=255)
                                                                               gggTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCGCACCGAGTGGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGCTCAGCGCACTTTTGTAACTGCTggg                                                                       <  1:156318/150‑1 (MQ=255)
                                                                                   gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTc                                                                   <  2:654539/150‑1 (MQ=255)
                                                                                     cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCtt                                                                 >  2:1283617/1‑150 (MQ=255)
                                                                                         gAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTaa                                                             <  1:1448019/150‑1 (MQ=255)
                                                                                           gCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAAtt                                                           <  2:1632533/150‑1 (MQ=255)
                                                                                            cAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTc                                                          >  2:834005/1‑150 (MQ=255)
                                                                                             aGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCg                                                         <  2:1364112/150‑1 (MQ=255)
                                                                                              gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGCCGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg                                                        >  1:1574456/1‑150 (MQ=255)
                                                                                              gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg                                                        <  1:1275775/150‑1 (MQ=255)
                                                                                              gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg                                                        >  2:329920/1‑150 (MQ=255)
                                                                                                     tGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTa                                                 >  2:425289/1‑150 (MQ=255)
                                                                                                        ccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGc                                              <  2:933376/150‑1 (MQ=255)
                                                                                                              aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGt                                        >  2:127730/1‑150 (MQ=255)
                                                                                                              aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGt                                        <  1:91727/150‑1 (MQ=255)
                                                                                                                tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCGCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCGGGGTGGCTTTTAATTCGGTGCGGTAGGCGATGGTGa                                      <  2:1365717/150‑1 (MQ=255)
                                                                                                                        ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCtttt                              <  1:1510116/150‑1 (MQ=255)
                                                                                                                         tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTg                             <  1:1117255/150‑1 (MQ=255)
                                                                                                                          aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGc                            >  1:1105727/1‑150 (MQ=255)
                                                                                                                               aCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGCCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTAGTAACTGCTGGGGGTCTTATAATTCCGAGCGGTAGGCGATCGTGATGACTTTTGCCAGGt                       >  1:1148266/1‑150 (MQ=255)
                                                                                                                               aCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGAAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGt                       >  2:308718/1‑150 (MQ=255)
                                                                                                                                    tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGc                  <  2:1051242/150‑1 (MQ=255)
                                                                                                                                      acccAGGTTATCGCGCTAACCTGGCTCAACGAGATGAACACGAAAATCCTACGCACTCTGGGAACGACGATCAGAGCAATTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATCGTGATGTCTTTTGCCAGGACGAGCAg                >  2:1168344/4‑150 (MQ=255)
                                                                                                                                       gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGt               >  1:257740/3‑150 (MQ=255)
                                                                                                                                       gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACAGCTGGGTGTCTTATAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGt               >  1:821060/3‑150 (MQ=255)
                                                                                                                                           ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTc           <  2:1105727/150‑1 (MQ=255)
                                                                                                                                            gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCg          >  1:1573700/1‑150 (MQ=255)
                                                                                                                                                ccGCGATAACCTGGCTCACCGAGTTGCAGCCGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTGTTGGAACTGCTGGGTGTCTTGTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGAAGTTGTCGACTg      >  1:222416/1‑150 (MQ=255)
                                                                                                                                                ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTg      >  1:74270/1‑150 (MQ=255)
                                                                                                                                                 cgcgGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTgg     >  1:1413010/1‑150 (MQ=255)
                                                                                                                                                 cgcgATAACGTGGCTCACCGAGTTGCAGACGATTATCCTGCGCACTCTGGGCAAGAGGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTGGGTGCGGTAGGCGATGGTGATGTCTTTTGCCTGGTCGAGCATTTGTCGACTgg     <  2:913778/150‑1 (MQ=255)
                                                                                                                                                  gcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGga    >  2:1016293/1‑150 (MQ=255)
                                                                                                                                                   cgATAACCTGGCTCACCGAGTTGCAGACGAAACTCCTGCGCACTCTGGGCAAGCCGAGCCGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGGCGAGCAGTTGTCGACTGgag   <  2:1539992/150‑1 (MQ=255)
                                                                                                                                                    gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga  <  1:1283617/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGA  >  NC_000913/1892214‑1892511

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