Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,161,966 C→A L66M (CTG→ATG)  yijD → DUF1422 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,161,9660CA98.6% 276.0 / NA 69L66M (CTG→ATG) yijDDUF1422 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (36/32);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTA  >  NC_000913/4161820‑4162104
                                                                                                                                                  |                                                                                                                                          
cTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                        >  2:210599/1‑151 (MQ=255)
 tGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGc                                                                                                                                        >  1:452789/1‑150 (MQ=255)
  ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCa                                                                                                                                       <  2:6017/150‑1 (MQ=255)
  ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCa                                                                                                                                       <  1:88620/150‑1 (MQ=255)
  ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCa                                                                                                                                       <  2:154152/150‑1 (MQ=255)
  ggCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCa                                                                                                                                       <  2:149460/150‑1 (MQ=255)
   gCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGc                                                                                                                                     <  1:678070/151‑1 (MQ=255)
       aTCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGcc                                                                                                                                 >  1:59573/1‑151 (MQ=255)
         cGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTg                                                                                                                               <  2:925662/151‑1 (MQ=255)
          gaatAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCGCTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGt                                                                                                                              <  2:45176/148‑1 (MQ=255)
           aTTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGtt                                                                                                                             >  2:619850/1‑151 (MQ=255)
           aTTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGtt                                                                                                                             <  2:103484/151‑1 (MQ=255)
           aTTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGtt                                                                                                                             >  1:453970/1‑151 (MQ=255)
            ttAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGttt                                                                                                                            >  2:597716/1‑151 (MQ=255)
              aaTGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGtttt                                                                                                                           <  1:138784/150‑1 (MQ=255)
               aTGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTAt                                                                                                                         >  1:491689/1‑151 (MQ=255)
                    aCTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCg                                                                                                                    >  1:576293/1‑151 (MQ=255)
                     cTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCg                                                                                                                    >  2:614988/1‑150 (MQ=255)
                      tttCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCgg                                                                                                                   >  1:947442/1‑150 (MQ=255)
                       ttCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCgg                                                                                                                   <  2:394055/149‑1 (MQ=255)
                           cAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTGGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTAct                                                                                                             >  2:702471/1‑151 (MQ=255)
                              cgcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                          <  1:108209/151‑1 (MQ=255)
                              cgcTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTActgct                                                                                                          <  2:654975/151‑1 (MQ=255)
                                 tGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGt                                                                                                        >  1:72186/1‑150 (MQ=255)
                                      aGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTCCAg                                                                                                    >  2:752602/1‑149 (MQ=255)
                                       gCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCAc                                                                                                 >  2:622262/1‑151 (MQ=255)
                                          ccATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACgg                                                                                               >  2:462478/1‑150 (MQ=255)
                                          ccATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACg                                                                                                <  1:739176/149‑1 (MQ=255)
                                           cATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATTGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                              >  1:333448/1‑150 (MQ=255)
                                            aTTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGt                                                                                              >  1:648527/1‑149 (MQ=255)
                                              tgtgCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGtt                                                                                          >  1:913757/1‑151 (MQ=255)
                                               gtgCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTc                                                                                         <  2:280462/151‑1 (MQ=255)
                                                 gCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGt                                                                                       <  1:951106/151‑1 (MQ=255)
                                                 gCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGt                                                                                       <  1:382459/151‑1 (MQ=255)
                                                           gTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATc                                                                             <  1:529366/151‑1 (MQ=255)
                                                             aTTCCCGATTATTTCCCTTGTGCTGACGGTTTACCGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGGCTGTTTTATTCTCGGCGTACTGCTGTACGGCACGGTAGTTCGTGCGGAATATcc                                                                            <  2:80953/150‑1 (MQ=255)
                                                             aTTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCg                                                                           <  1:421649/151‑1 (MQ=255)
                                                                 ccGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGata                                                                       <  2:629966/151‑1 (MQ=255)
                                                                    attattTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATc                                                                     >  2:566975/1‑150 (MQ=255)
                                                                    attattTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCg                                                                    >  1:37941/1‑151 (MQ=255)
                                                                         ttCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctct                                                               >  2:934682/1‑151 (MQ=255)
                                                                          tCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctct                                                               >  1:843611/1‑150 (MQ=255)
                                                                          tCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGctct                                                               <  2:607746/150‑1 (MQ=255)
                                                                              tGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAActt                                                          >  1:670919/1‑151 (MQ=255)
                                                                                   cTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTcc                                                     >  2:327357/1‑151 (MQ=255)
                                                                                         gTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGc                                                 >  1:206047/1‑149 (MQ=255)
                                                                                           ttACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGt                                              >  2:576917/1‑150 (MQ=255)
                                                                                            tACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTAc                                            >  1:566687/1‑151 (MQ=255)
                                                                                             aCTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTAc                                            <  2:143521/150‑1 (MQ=255)
                                                                                                gcctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACtctc                                        >  2:51730/1‑151 (MQ=255)
                                                                                                 cctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTcc                                       <  2:383805/151‑1 (MQ=255)
                                                                                                  ctgcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCg                                      >  1:866570/1‑151 (MQ=255)
                                                                                                    gcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTc                                    >  1:905122/1‑151 (MQ=255)
                                                                                                    gcATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTGCCCGGCAGTACTCTCCGt                                     <  2:648577/150‑1 (MQ=255)
                                                                                                      aTCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAt                                  >  1:928635/1‑151 (MQ=255)
                                                                                                       tCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCAtt                                 >  1:842987/1‑151 (MQ=255)
                                                                                                          aCGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATg                              <  2:448702/151‑1 (MQ=255)
                                                                                                            gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGgt                            >  1:537136/1‑151 (MQ=255)
                                                                                                            gTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGgt                            <  2:267752/151‑1 (MQ=255)
                                                                                                              tATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGgtgt                          >  1:339354/1‑151 (MQ=255)
                                                                                                                       cGCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTg                 >  2:39042/1‑151 (MQ=255)
                                                                                                                        gCACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTgg                >  1:645797/1‑151 (MQ=255)
                                                                                                                         cACTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGgc               <  2:892875/151‑1 (MQ=255)
                                                                                                                          aCTATGCCGGTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGgcg              <  1:560346/151‑1 (MQ=255)
                                                                                                                                   gTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGc     <  1:831698/151‑1 (MQ=255)
                                                                                                                                   gTAGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGc     <  1:101979/151‑1 (MQ=255)
                                                                                                                                     aGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCg    <  2:204968/150‑1 (MQ=255)
                                                                                                                                     aGGCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCg    <  1:125509/150‑1 (MQ=255)
                                                                                                                                      ggCTTGCCGGGTATGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTa  <  2:893334/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
CTGGCTTATCGATTAATGGTACTTTCGCAGCGCTGTTTAGCTCCATTGTGCCATTTTCTGTATTCCCGATTATTTCCCTGGTGCTGACGGTTTACTGCCTGCATCAACGTTATCTTAATCGCACTATGCCGGTAGGCTTGCCGGGTCTGGCAGCTGCCTGTTTTATTCTCGGCGTACTGCTGTACAGCACGGTAGTTCGTGCGGAATATCCGGATATCGGCTCTAACTTCTTCCCGGCAGTACTCTCCGTCATTATGGTGTTCTGGATTGGCGCGAAGATGCGTA  >  NC_000913/4161820‑4162104

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑