Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,369,301 (T)6→7 coding (178/465 nt) yhcH ← DUF386 domain‑containing protein YhcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,369,2951.T98.6% 242.7 / ‑3.6 69coding (184/465 nt)yhcHDUF386 domain‑containing protein YhcH
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (34/34);  total (35/34)
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.75e-01

TGCTGCAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTGTAACG  >  NC_000913/3369163‑3369429
                                                                                                                                     |                                                                                                                                      
tgctgcAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCttttttt                                                                                                                                 >  1:404314/1‑139 (MQ=255)
 gctgcAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCtttttttc                                                                                                                                >  1:379696/1‑138 (MQ=255)
   tgcAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAGTTCCGCtttttttctc                                                                                                                              >  2:775736/1‑136 (MQ=255)
   tgcAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCtttttttctc                                                                                                                              >  1:783550/1‑136 (MQ=255)
   tgcAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCtttttttctc                                                                                                                              >  2:196864/1‑136 (MQ=255)
     cAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGa                                                                                                                             >  1:771186/1‑140 (MQ=255)
      aaaGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGAc                                                                                                                            <  2:472504/140‑1 (MQ=255)
         gCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACggg                                                                                                                         <  2:404314/140‑1 (MQ=255)
          cTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGc                                                                                                                        >  2:518868/1‑140 (MQ=255)
            gATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGa                                                                                                                      >  1:278423/1‑140 (MQ=255)
            gATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGa                                                                                                                      >  2:325518/1‑140 (MQ=255)
                 tcatcCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAg                                                                                                                 <  2:43231/140‑1 (MQ=255)
                  catcCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAgt                                                                                                                >  1:239581/1‑140 (MQ=255)
                    tcCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGAtt                                                                                                                    >  2:170496/1‑134 (MQ=255)
                    tcCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGAtt                                                                                                                    <  1:170496/134‑1 (MQ=255)
                    tcCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAgtgt                                                                                                              >  2:691619/1‑140 (MQ=255)
                       tCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTaa                                                                                                           >  2:326982/1‑140 (MQ=255)
                          tggtggAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGa                                                                                                                  <  1:337455/130‑1 (MQ=255)
                          tggtggAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGa                                                                                                                  >  2:337455/1‑130 (MQ=255)
                               gAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTa                                                                                                            >  1:88394/1‑131 (MQ=255)
                               gAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTa                                                                                                            <  2:88394/131‑1 (MQ=255)
                                aaCTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGa                                                                                                  <  1:627701/140‑1 (MQ=255)
                                 aCTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGAc                                                                                                 <  2:440632/140‑1 (MQ=255)
                                  ctctTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACa                                                                                                <  2:771186/140‑1 (MQ=255)
                                  ctctTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACa                                                                                                <  1:325518/140‑1 (MQ=255)
                                  ctctTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACa                                                                                                <  2:115727/140‑1 (MQ=255)
                                      tCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCa                                                                                            <  1:326982/140‑1 (MQ=255)
                                      tCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCa                                                                                            <  2:278423/140‑1 (MQ=255)
                                       cacaCTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCAt                                                                                           <  2:265093/140‑1 (MQ=255)
                                       cacaCTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCAt                                                                                           >  2:244181/1‑140 (MQ=255)
                                        acacTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATa                                                                                          <  1:196864/140‑1 (MQ=255)
                                          acTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATaaa                                                                                        >  2:266850/1‑140 (MQ=255)
                                              aCGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCAt                                                                                           <  1:149695/133‑1 (MQ=255)
                                              aCGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCAt                                                                                           >  2:149695/1‑133 (MQ=255)
                                               cgcgCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGa                                                                                                  >  2:41675/1‑125 (MQ=255)
                                               cgcgCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGa                                                                                                  <  1:41675/125‑1 (MQ=255)
                                               cgcgCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATa                                                                                          <  1:287815/133‑1 (MQ=255)
                                               cgcgCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATa                                                                                          >  2:287815/1‑133 (MQ=255)
                                                 cgcAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTg                                                                                 <  2:783550/140‑1 (MQ=255)
                                                    aGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCg                                                                              <  2:349379/140‑1 (MQ=255)
                                                    aGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCg                                                                              <  1:691619/140‑1 (MQ=255)
                                                                gccaAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCg                                                                  >  1:495793/1‑140 (MQ=255)
                                                                gccaAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCg                                                                  >  2:412142/1‑140 (MQ=255)
                                                                gccaAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCg                                                                  >  1:407256/1‑140 (MQ=255)
                                                                    aaCAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAAc                                                              >  1:803170/1‑140 (MQ=255)
                                                                        gAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGAcc                                                          <  1:266850/140‑1 (MQ=255)
                                                                               ttCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCt                                                   >  1:50914/1‑140 (MQ=255)
                                                                                tCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCtt                                                  >  2:518756/1‑140 (MQ=255)
                                                                                    cACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCt                                             >  1:830116/1‑140 (MQ=255)
                                                                                     aCCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCt                                             <  2:239581/140‑1 (MQ=255)
                                                                                       cGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTg                                           >  1:381184/1‑140 (MQ=255)
                                                                                       cGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTg                                           >  1:500463/1‑140 (MQ=255)
                                                                                               aaCAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGCCTgg                                   <  1:518868/140‑1 (MQ=255)
                                                                                                   gCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGc                               <  1:411952/140‑1 (MQ=255)
                                                                                                   gCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGc                               <  2:803170/140‑1 (MQ=255)
                                                                                                   gCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGc                               >  2:711268/1‑140 (MQ=255)
                                                                                                    cTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCt                              >  1:566899/1‑140 (MQ=255)
                                                                                                    cTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCt                              >  1:106512/1‑140 (MQ=255)
                                                                                                        atatCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATg                          >  1:683357/1‑140 (MQ=255)
                                                                                                         tatCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGc                         <  1:711268/140‑1 (MQ=255)
                                                                                                           tCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGt                                       >  2:461054/1‑124 (MQ=255)
                                                                                                           tCAATGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGt                                       <  1:461054/124‑1 (MQ=255)
                                                                                                              aTGTATTGCTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCg                    <  2:683357/140‑1 (MQ=255)
                                                                                                                      cTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAAcgcg            <  1:56681/140‑1 (MQ=255)
                                                                                                                      cTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAAcgcg            >  1:649425/1‑140 (MQ=255)
                                                                                                                      cTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAAcgcg            >  2:128316/1‑140 (MQ=255)
                                                                                                                      cTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAAcg              >  2:435722/1‑138 (MQ=255)
                                                                                                                      cTCGTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAAcg              <  1:435722/138‑1 (MQ=255)
                                                                                                                         gTGCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTcc         >  1:325551/1‑140 (MQ=255)
                                                                                                                           gCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTg       >  2:551850/1‑140 (MQ=255)
                                                                                                                           gCAATTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTg       >  1:814172/1‑140 (MQ=255)
                                                                                                                              aTTCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTGTaa    <  1:465716/140‑1 (MQ=255)
                                                                                                                                tCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTGTAACg  <  2:106512/140‑1 (MQ=255)
                                                                                                                                tCCGCTTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTGTAACg  <  2:814172/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TGCTGCAAAGCTGATAATCATCCTCATGGTGGAACTCTTCACACTGACGCGCAGTGCCTGCCATGCCAAACAGAATCCGTTCCTCACCGTTTAATAACAGCTGGATATCAATGTATTGCTCGTGCAATTCCGCTTTTTTCTCGACGGGCGATTGAGTGTTAAACGTCATGACATTCATAAAGATATTGTCGCCCTGTAATTCGTAACGACCCGGCGCTTTTTCTTGCGGTCTGGCAGCTAATGCCAGCGTTAACGCGTCCTGTAACG  >  NC_000913/3369163‑3369429

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