Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,047 C→T 61.7% intergenic (‑306/+352) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,0470CT61.7% 38.6 / 171.8 173intergenic (‑306/+352)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base T (44/62);  total (80/93)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 174 total aligned reads displayed.

GTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA  >  NC_000913/4325990‑4326145
                                                         |                                                                                                  
gTGCGGGAATTTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          >  2:2724589/1‑100 (MQ=14)
gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          <  2:726888/100‑1 (MQ=14)
gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          <  2:2206870/100‑1 (MQ=31)
gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          <  2:2402412/100‑1 (MQ=31)
 tGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCAGATGGGGCGTAAGCGCCTTAGCCGACCTGCGGATg                                                          >  1:1088362/1‑99 (MQ=11)
 tGCGGGAATTTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATg                                                          <  1:1491438/99‑1 (MQ=14)
   caagAATTTGTAGGACTGCGAAGACGGCCACGCGCCGCATCAGTCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGATGGTg                                                       <  1:481526/97‑1 (MQ=11)
   cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTg                                                       >  1:1272794/1‑100 (MQ=17)
   cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACggtggacg                                                       >  1:2542306/1‑94 (MQ=0)
    gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGc                                                      >  2:582923/1‑100 (MQ=17)
     ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                     >  2:2964368/1‑100 (MQ=17)
     ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                     >  2:13276/1‑100 (MQ=18)
     ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTg                                                       >  1:1658491/1‑98 (MQ=17)
      gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCgg                                                    <  1:2644802/100‑1 (MQ=18)
      gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCgg                                                    <  1:363333/100‑1 (MQ=33)
       aaTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCggg                                                   <  2:3316373/100‑1 (MQ=18)
       aaTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCggg                                                   >  1:2590781/1‑100 (MQ=17)
        acttGTAGGCCTGATAAGGCGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg                                                     <  1:3203551/95‑1 (MQ=255)
         tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGaa                                                 <  2:3149419/100‑1 (MQ=18)
         tttGTACGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGaa                                                 <  1:2331128/100‑1 (MQ=16)
          ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAt                                                >  2:1744812/1‑100 (MQ=18)
          ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAt                                                >  2:876301/1‑100 (MQ=17)
           tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAtt                                               <  1:963470/100‑1 (MQ=18)
           tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAtn                                               <  2:3258787/100‑2 (MQ=18)
            gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                              >  2:1439718/1‑100 (MQ=17)
            gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAttt                                              >  2:3316222/1‑100 (MQ=18)
             tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  2:1311957/100‑1 (MQ=18)
             tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  1:975300/100‑1 (MQ=17)
             tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  2:475725/100‑1 (MQ=18)
             tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                             <  2:2376301/100‑1 (MQ=33)
              aGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                            <  2:2919611/100‑1 (MQ=18)
                gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                            <  1:2007631/98‑1 (MQ=17)
                gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                          <  2:3193726/100‑1 (MQ=18)
                gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                          <  2:766521/100‑1 (MQ=33)
                 tcTGAGAAGACTGGCAAGTGGCGCATCAGGCGGTCGGCACTGTTGCCGGATGCAGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                            <  1:2420132/96‑1 (MQ=11)
                 cgtGAGAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  2:700505/98‑1 (MQ=255)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  2:134118/100‑1 (MQ=14)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  2:2591428/100‑1 (MQ=18)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  1:2193615/100‑1 (MQ=18)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  1:3442484/100‑1 (MQ=17)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         <  2:637596/100‑1 (MQ=18)
                 ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAgg                                         >  2:2238169/1‑100 (MQ=32)
                  cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGc                                        <  2:1413557/100‑1 (MQ=18)
                   tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                       <  1:3080753/100‑1 (MQ=33)
                     aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTg                                     <  1:2621024/100‑1 (MQ=18)
                      tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                    <  1:1906709/100‑1 (MQ=18)
                      tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                    <  1:2910200/100‑1 (MQ=18)
                       aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                    >  1:919197/1‑99 (MQ=17)
                       aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                   <  2:3123831/100‑1 (MQ=17)
                       aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                   <  2:1205433/100‑1 (MQ=18)
                        aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTGAACGGCTGAGCCGCCCTTCGGGTGGTGCGGGAATATGGAGTCCTgct                                   >  1:1025887/1‑97 (MQ=11)
                        aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                   >  1:396857/1‑99 (MQ=32)
                         gACGCGCAAGCGTCGCATCATGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 <  1:3162348/100‑1 (MQ=16)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAtga                                 >  2:1736419/1‑98 (MQ=16)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 <  2:844824/100‑1 (MQ=18)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 <  1:1277745/100‑1 (MQ=17)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  1:818847/1‑100 (MQ=17)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  2:1170447/1‑100 (MQ=17)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATaa                                 >  1:3226282/1‑100 (MQ=12)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGGATTTGTAGGCCTGAt                                   >  2:1585864/1‑98 (MQ=17)
                         gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATaa                                 >  2:2901747/1‑100 (MQ=33)
                         gACGCGCAAGCGCCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCt                                      <  2:1091424/95‑1 (MQ=18)
                          aCGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                <  1:2263707/100‑1 (MQ=18)
                          aCGCGCAAGCGTCGCATCAGGCAGTCCGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                  <  1:2706362/98‑1 (MQ=17)
                           cgcgCAAGCGTCGTATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:2322164/1‑100 (MQ=18)
                           cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAGCGCCTTATCCGGCCTACGGCTGGTGCGGGTATTTGTAGGCCTGATaa                                 >  1:1620558/1‑98 (MQ=11)
                           cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAAGGCCTTATCCGGCCTAGGGATGGTGTGGGAATGTGTAGAGCTGATaa                                 >  1:536/1‑98 (MQ=11)
                           cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:400585/1‑100 (MQ=17)
                           cgcgCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTCCGGGTGGTGCGGGAATTTGTGGGccggggttgg                               >  1:2086242/1‑92 (MQ=11)
                           cgcgCAAGCGTCGCATCAGGCAGTCGACACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGa                               >  2:3386969/1‑100 (MQ=255)
                            gcgcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAg                                <  1:3314461/98‑1 (MQ=18)
                            gcgcAAGCGTCGCATCAGGCAGCCGGCACCGTTGCCGGATGCGGCGTAAGCGCCTTATCCGGCCTACGGATGGCGCGGCAATTTGCAGACCTgctcaaa                               >  1:977292/1‑93 (MQ=11)
                             cgcAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcg                             >  2:1261175/1‑100 (MQ=17)
                              gcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                            >  2:2046951/1‑100 (MQ=33)
                               cAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGATGGTGCGGGAAGTTGTAGGCCCGATCAGAcgag                           >  1:275057/1‑98 (MQ=11)
                               cAAGCGTCGCATCAGGCAGTCGGCACCGTTGTCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                           <  2:51987/100‑1 (MQ=17)
                               cAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg                           <  2:1858761/100‑1 (MQ=33)
                                aaGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                          <  2:1315481/100‑1 (MQ=14)
                                aaGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                          >  2:2792729/1‑100 (MQ=11)
                                aaGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                          >  1:2343706/1‑100 (MQ=1)
                                    gTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAg                       <  1:2760115/99‑1 (MQ=32)
                                    gTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                      >  2:1137187/1‑100 (MQ=33)
                                     tCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                     >  2:2676316/1‑100 (MQ=33)
                                     tCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                     <  2:1244136/100‑1 (MQ=33)
                                        cATCAGGCCGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                  <  1:2988758/100‑1 (MQ=17)
                                        cATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                  <  1:1216968/100‑1 (MQ=33)
                                         aTCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                 <  2:3167911/100‑1 (MQ=33)
                                           cAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCAt               >  2:1490305/1‑100 (MQ=32)
                                             ggCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa             >  1:732310/1‑100 (MQ=32)
                                              gCAGTCGGCACCGTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTCTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg            <  1:2083357/100‑1 (MQ=16)
                                               cAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg           <  2:1173933/100‑1 (MQ=18)
                                                aGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc          <  2:768851/100‑1 (MQ=33)
                                                    ggCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt      >  2:53063/1‑100 (MQ=11)
                                                    ggCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt      >  2:1635315/1‑100 (MQ=1)
                                                    ggCACTGTTGCCGGATGCGGCGTAAACGCCTGATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt      <  2:3091560/100‑1 (MQ=2)
                                                    ggCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGGATTTGTACGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGtt      >  2:2429486/1‑100 (MQ=255)
                                                     gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTg     <  2:1185621/100‑1 (MQ=14)
                                                     gCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTg     <  1:391930/100‑1 (MQ=38)
                                                      cACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGGTAAGACGCGCGTGCGTGGCATCAGGCAGTTgg    >  1:1590339/1‑100 (MQ=14)
                                                      cACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTgg    >  2:4857/1‑100 (MQ=40)
                                                         |                                                                                                  
GTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCA  >  NC_000913/4325990‑4326145

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