Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 7.5% intergenic (+28/+32) yhcG → / ← yhcH DUF1016 domain‑containing protein YhcG/DUF386 domain‑containing protein YhcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,9820GC 7.5% 447.9 / 20.8 133intergenic (+28/+32)yhcG/yhcHDUF1016 domain‑containing protein YhcG/DUF386 domain‑containing protein YhcH
Reads supporting (aligned to +/- strand):  ref base G (71/52);  new base C (6/4);  total (77/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 137 total aligned reads displayed.

TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGGT  >  NC_000913/3368837‑3369131
                                                                                                                                                    |                                                                                                                                                     
tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTc                                                                                                                                                   <  2:468907/150‑1 (MQ=255)
tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTc                                                                                                                                                   <  2:1242753/150‑1 (MQ=255)
   gCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGc                                                                                                                                                <  2:868000/150‑1 (MQ=255)
     cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACACTGTCTGCAg                                                                                                                                              <  2:412168/150‑1 (MQ=255)
      tAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACACTGTCTGCAGa                                                                                                                                             <  1:1245422/150‑1 (MQ=255)
        acacCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACACTGTCTGCAGAAt                                                                                                                                           <  1:372741/150‑1 (MQ=255)
             gTTCTGGCAAGCGAATATCAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTa                                                                                                                                      <  2:1577161/150‑1 (MQ=255)
             gTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTa                                                                                                                                      >  1:204797/1‑150 (MQ=255)
               tCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTaaa                                                                                                                                    <  2:1411738/150‑1 (MQ=255)
                 tGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAAc                                                                                                                                  <  1:16439/150‑1 (MQ=255)
                  ggCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACg                                                                                                                                 <  1:355262/150‑1 (MQ=255)
                    cAAGCGAATAAAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                               <  2:1639508/150‑1 (MQ=255)
                      aGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACt                                                                                                                             >  2:147153/1‑150 (MQ=255)
                      aGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCAGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAAAGAACt                                                                                                                             >  2:1247784/1‑150 (MQ=255)
                         gAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAg                                                                                                                          <  1:626670/150‑1 (MQ=255)
                          aaTATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGt                                                                                                                         >  2:1056044/1‑150 (MQ=255)
                           atatAAGATGCAACTACCTGATGAGAAACGACTCGCAGATTAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGtt                                                                                                                        <  2:1404542/150‑1 (MQ=255)
                            tataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGttt                                                                                                                       <  1:1422513/150‑1 (MQ=255)
                              taAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTa                                                                                                                     >  2:1610211/1‑150 (MQ=255)
                                aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                   >  2:874121/1‑150 (MQ=255)
                                aGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                   >  1:786512/1‑150 (MQ=255)
                                   tGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCa                                                                                                                >  1:99247/1‑150 (MQ=255)
                                   tGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCa                                                                                                                <  1:593219/150‑1 (MQ=255)
                                    gCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAt                                                                                                               <  2:65932/150‑1 (MQ=255)
                                      aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGTCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTa                                                                                                             <  1:346256/150‑1 (MQ=255)
                                          aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                         <  2:1136182/150‑1 (MQ=255)
                                          aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                         >  1:778582/1‑150 (MQ=255)
                                          aCCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATc                                                                                                         >  1:1439090/1‑150 (MQ=255)
                                           ccTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAAGCa                                                                                                        >  1:1399258/1‑150 (MQ=255)
                                              gatgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCAGTGCTAGAGGAAGGCTATAGACGCCATTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                                                                     >  2:1424074/1‑150 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                                  >  1:267793/1‑150 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                                  >  2:800845/1‑150 (MQ=255)
                                                 gagaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTa                                                                                                  <  1:153052/150‑1 (MQ=255)
                                                  agaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                 <  2:1461143/150‑1 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                                >  1:1564973/1‑150 (MQ=255)
                                                   gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                                <  2:150237/150‑1 (MQ=255)
                                                      aCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCtt                                                                                             >  1:879465/1‑150 (MQ=255)
                                                        gACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcc                                                                                           >  2:554370/1‑149 (MQ=255)
                                                         aCTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAAACTTAGTTCGTTTTAAAGATTCTGCAGACAGTGTCTGCAGAATATTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAACtt                                                                                             >  2:354791/1‑150 (MQ=25)
                                                            cGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcaccac                                                                                       >  1:365085/1‑150 (MQ=255)
                                                              cAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACaa                                                                                     >  2:1588617/1‑150 (MQ=255)
                                                                    aaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                               >  1:1034114/1‑150 (MQ=255)
                                                                    aaCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                               >  2:2900/1‑150 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                              >  2:1214445/1‑150 (MQ=255)
                                                                     aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                              <  2:267793/150‑1 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTACGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                             >  2:958140/1‑150 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                             <  1:1465269/150‑1 (MQ=255)
                                                                      cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTa                                                                             <  1:146396/150‑1 (MQ=255)
                                                                       tCGTGCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                            <  1:1431685/150‑1 (MQ=255)
                                                                         gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACCGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtc                                                                          <  2:323773/150‑1 (MQ=255)
                                                                         gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtc                                                                          >  2:222376/1‑150 (MQ=255)
                                                                           tCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCATCCTTAACCTTCACCAAAACCGTTTTATtcta                                                                        >  2:733958/1‑149 (MQ=255)
                                                                             gAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctctc                                                                      >  1:172825/1‑150 (MQ=255)
                                                                               acacAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCa                                                                    <  1:455131/150‑1 (MQ=255)
                                                                                   aGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                               <  2:970946/135‑1 (MQ=255)
                                                                                   aGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                               >  1:970946/1‑135 (MQ=255)
                                                                                    ggCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTc                                                               >  1:181766/1‑150 (MQ=255)
                                                                                    ggCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTc                                                               >  1:918931/1‑150 (MQ=255)
                                                                                            tAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCCCAACCTTTTTAATCTCTCCAGGCTCGCCGACaa                                                       >  2:441307/1‑150 (MQ=255)
                                                                                             agagGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                      >  1:890065/1‑150 (MQ=255)
                                                                                                ggAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCa                                                   >  2:842227/1‑150 (MQ=255)
                                                                                                ggAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCa                                                   >  1:1630767/1‑150 (MQ=255)
                                                                                                      cTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTgg                                             >  2:474284/1‑150 (MQ=255)
                                                                                                       tataGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCATAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGt                                            >  2:706874/1‑150 (MQ=255)
                                                                                                        ataGACTCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTTCAGACAGTGTCTTCAGAATCTTTAAAACGTATTAATTTTTAAGCGATTAAATCAGCCTTAACCTTGACCACAACCTTTTTAATCTCTCGAGGCTCGCCGACTACGAATGCTGGtt                                           <  2:204797/150‑1 (MQ=255)
                                                                                                        ataGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGtt                                           <  1:874121/150‑1 (MQ=255)
                                                                                                         taGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGttt                                          <  2:662611/150‑1 (MQ=255)
                                                                                                           gACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTAt                                        >  2:825921/1‑150 (MQ=255)
                                                                                                             cGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGc                                      <  1:1588617/150‑1 (MQ=255)
                                                                                                               cctttAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                                    <  2:1439090/147‑1 (MQ=255)
                                                                                                                   tAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCa                                >  1:873445/1‑150 (MQ=255)
                                                                                                                         cAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGc                          <  2:778582/150‑1 (MQ=255)
                                                                                                                              cGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATa                       <  1:1323163/150‑1 (MQ=255)
                                                                                                                               ggTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaa                      <  2:1034114/150‑1 (MQ=255)
                                                                                                                                gTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                     >  1:957024/1‑150 (MQ=255)
                                                                                                                                gTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATaaa                     >  1:1021025/1‑150 (MQ=255)
                                                                                                                                  ttCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACTTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                   >  2:117201/1‑150 (MQ=255)
                                                                                                                                  ttCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                   >  1:1118798/1‑150 (MQ=255)
                                                                                                                                  ttCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                   <  1:15085/150‑1 (MQ=255)
                                                                                                                                  ttCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAaca                   >  1:657981/1‑150 (MQ=255)
                                                                                                                                   tCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAacac                  >  1:1023178/1‑150 (MQ=255)
                                                                                                                                      a‑aTTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGc               >  1:759904/1‑150 (MQ=255)
                                                                                                                                       gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGc               >  2:437104/2‑150 (MQ=255)
                                                                                                                                        aTTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACATTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCg              >  2:93527/1‑150 (MQ=255)
                                                                                                                                         ttCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGa             <  2:195100/150‑1 (MQ=255)
                                                                                                                                         ttCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGa             >  1:235504/1‑150 (MQ=255)
                                                                                                                                          tCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGaa            >  1:600290/1‑150 (MQ=255)
                                                                                                                                          tCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGaa            >  1:963117/1‑150 (MQ=255)
                                                                                                                                           cTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc           >  2:106758/1‑150 (MQ=255)
                                                                                                                                           cTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACATTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAAc           >  2:278573/1‑150 (MQ=255)
                                                                                                                                            tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACa          >  1:1388297/1‑150 (MQ=255)
                                                                                                                                            tGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTAACCTGGCATAAACACGGCGAACa          >  2:1498581/1‑150 (MQ=255)
                                                                                                                                             gCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACAt         >  1:96109/1‑150 (MQ=255)
                                                                                                                                              cAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACAtt        >  1:1140621/1‑150 (MQ=255)
                                                                                                                                               aGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAACTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACAAGGCGAACATTc       >  2:107095/1‑150 (MQ=255)
                                                                                                                                                gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTcc      >  1:643686/1‑150 (MQ=255)
                                                                                                                                                 aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTccc     >  1:1534655/1‑150 (MQ=255)
                                                                                                                                                  cAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCg    >  1:1652176/1‑150 (MQ=255)
                                                                                                                                                    gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGGt  >  2:1656369/1‑150 (MQ=255)
                                                                                                                                                    gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGGt  >  1:19397/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAG‑‑CGGTTTCGA‑ATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGGT  >  NC_000913/3368837‑3369131

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: ‑