Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,892,358 C→T 6.2% G293D (GGT→GAT)  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,3580CT 6.2% 245.4 / 4.1 80G293D (GGT→GAT) yoaAputative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base T (3/2);  total (41/39)
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.59e-01

ACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTG  >  NC_000913/1892212‑1892507
                                                                                                                                                  |                                                                                                                                                     
aCAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                                                                    <  2:227165/150‑1 (MQ=255)
aCAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                                                                    <  2:1377240/150‑1 (MQ=255)
aCAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCACTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                                                                    <  1:1579621/150‑1 (MQ=255)
 cAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                                                                                                   >  2:1510116/1‑150 (MQ=255)
  aaCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGa                                                                                                                                                  <  2:598621/150‑1 (MQ=255)
   aCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                                                                                                 <  1:39275/150‑1 (MQ=255)
   aCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                                                                                                 <  1:312707/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:1291986/150‑1 (MQ=255)
                                                           gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                         >  2:1275775/1‑150 (MQ=255)
                                                           gCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                         >  1:1632533/1‑150 (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:1438696/1‑150 (MQ=255)
                                                                    tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACttt                                                                                >  1:40629/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)
                                                                                                                                                  |                                                                                                                                                     
ACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTG  >  NC_000913/1892212‑1892507

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: