Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,666,629 T→C T64A (ACG→GCG)  yhjA ← putative cytochrome C peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,666,6290TC99.7% 1987.6 / ‑4.6 619T64A (ACG→GCG) yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (311/306);  total (313/306)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.27e-01

Only 100 of 651 total aligned reads displayed.

TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666539‑3666718
                                                                                          |                                                                                         
tAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCgcgt                                                                                         >  1:903812/1‑93 (MQ=255)
 aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGcgtg                                                                                        >  2:1520702/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       >  1:1581184/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       <  1:6353442/93‑1 (MQ=255)
   aaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGc                                                                                      >  1:5918917/1‑93 (MQ=255)
      aTTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAAt                                                                                   >  1:9552851/1‑93 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  <  1:5224360/93‑1 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  >  2:2784616/1‑93 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  <  2:1800244/93‑1 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  <  1:6252432/93‑1 (MQ=255)
        gtataTCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCGGCCGAAGGCGCGTTGCAATAg                                                                                 >  2:6164516/2‑93 (MQ=255)
         tataTCCAAGCTTAATATCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGt                                                                                >  2:7657218/1‑93 (MQ=255)
          atatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                               >  1:5764768/1‑93 (MQ=255)
           tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCg                                                                              >  2:3619949/1‑93 (MQ=255)
           tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCg                                                                              >  2:2879639/1‑93 (MQ=255)
           tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCAAGGCGCGTGGCAATAGTCGc                                                                             <  2:1325879/93‑1 (MQ=255)
               cAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATc                                                                          >  2:6211105/1‑93 (MQ=255)
                 aGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCt                                                                        >  1:3932978/1‑93 (MQ=255)
                   cTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttt                                                                      <  1:1184681/93‑1 (MQ=255)
                     tAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttttt                                                                    >  2:5913074/1‑93 (MQ=255)
                     tAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttttt                                                                    <  2:3102460/93‑1 (MQ=255)
                     tAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGCGGCGTGGCAATAGTCGCATCCttttt                                                                    >  2:706965/1‑93 (MQ=255)
                       aTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcg                                                                  >  1:7277360/1‑93 (MQ=255)
                        tGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcgc                                                                 <  2:939858/93‑1 (MQ=255)
                          tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               >  2:3530421/1‑93 (MQ=255)
                             tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGGAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                            >  2:1226458/1‑93 (MQ=255)
                               aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  1:5513716/93‑1 (MQ=255)
                               aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  1:1049206/93‑1 (MQ=255)
                                tCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAg                                                         <  2:5202882/93‑1 (MQ=255)
                                tCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAg                                                         >  1:185377/1‑93 (MQ=255)
                                   aTCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCt                                                      >  2:7942688/1‑93 (MQ=255)
                                   aTCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCt                                                      >  2:8356041/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  2:9202119/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  1:2072871/1‑93 (MQ=255)
                                      aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   <  2:8758762/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  >  2:5740086/1‑93 (MQ=255)
                                        cTGTTGCGCGCCAGGAATATAATAATAGGCGGGTAATGCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAc                                                 <  2:4497306/93‑1 (MQ=255)
                                          gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               >  1:2714004/1‑93 (MQ=255)
                                          gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               >  1:6466087/1‑93 (MQ=255)
                                          gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               >  1:6954664/1‑93 (MQ=255)
                                           tttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtt                                              >  2:246100/1‑93 (MQ=255)
                                           tttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtt                                              <  1:9971192/93‑1 (MQ=255)
                                            ttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtta                                             <  1:7120798/93‑1 (MQ=255)
                                             tCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttat                                            <  1:9713437/93‑1 (MQ=255)
                                             tCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttat                                            >  2:9714171/1‑93 (MQ=255)
                                             tCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttat                                            <  1:2305303/93‑1 (MQ=255)
                                              cgcgCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCATAAAACCttatt                                           >  1:8248812/1‑93 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                        <  1:4718718/93‑1 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                        <  2:3376726/93‑1 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                        <  2:5122040/93‑1 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                        >  2:6405431/1‑93 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                       >  1:3489141/1‑93 (MQ=255)
                                                   cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                      >  1:5634622/1‑93 (MQ=255)
                                                   cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                      >  2:4351838/1‑93 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                     >  1:7565689/1‑93 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                     >  1:8569459/1‑93 (MQ=255)
                                                     ggAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCa                                    >  2:6634529/1‑93 (MQ=255)
                                                      gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   <  1:7888117/93‑1 (MQ=255)
                                                      gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   <  1:4090073/93‑1 (MQ=255)
                                                       aaTATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACt                                  >  2:959022/1‑93 (MQ=255)
                                                        atatAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACta                                 >  1:6006831/1‑92 (MQ=255)
                                                        atatAATAATAGGCGGGTAATTCTGCCGAAGGCGCGCGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTg                                 >  2:2666443/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                >  2:679207/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                >  1:4860135/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                >  1:2517646/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                <  1:2476588/93‑1 (MQ=255)
                                                          ataataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTCATTATTTTCACTGAc                               >  2:8969760/1‑93 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              >  1:2820465/1‑93 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              <  2:204207/93‑1 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              >  2:9085583/1‑93 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              <  1:9362134/93‑1 (MQ=255)
                                                            aataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcag                             <  1:2939428/93‑1 (MQ=255)
                                                             ataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAGAGCCTAAAACCTTATTATTTTCACTGAcagc                            <  2:1916389/93‑1 (MQ=255)
                                                              taataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagca                           >  2:3780121/1‑93 (MQ=255)
                                                              taataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagca                           >  1:3774518/1‑93 (MQ=255)
                                                                ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  2:8544642/93‑1 (MQ=255)
                                                                 taGGCGGGCAATTCTGCCGAAGGAGCGAGGCAAAAGACGACAGCTTTTTCGCAGAGGAAGCCTACAACCATAGTAATATCACTGCCAACACAt                        >  2:7333079/1‑93 (MQ=255)
                                                                  aGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATg                       <  1:10008314/93‑1 (MQ=255)
                                                                  aGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATg                       >  1:8736639/1‑93 (MQ=255)
                                                                  aGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGACTAAAACCTTATTATTTGCACTGACGGCCGacg                       >  1:3029042/1‑91 (MQ=255)
                                                                   ggCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGc                      <  2:5638119/93‑1 (MQ=255)
                                                                    gCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGcc                     <  1:3614621/93‑1 (MQ=255)
                                                                       ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                  >  1:8091985/1‑93 (MQ=255)
                                                                        gTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGaa                 <  1:3215602/93‑1 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAAc                >  1:1370220/1‑93 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTAACTGACAGCAGATGCCTGAAc                >  1:6825823/1‑93 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              <  1:5051972/93‑1 (MQ=255)
                                                                            ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             >  2:4765079/1‑93 (MQ=255)
                                                                            ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             <  2:3937013/93‑1 (MQ=255)
                                                                             tCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTCTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCg            <  2:4611296/93‑1 (MQ=255)
                                                                              cTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCgg           <  2:1115804/93‑1 (MQ=255)
                                                                                gCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGACTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGcc         >  2:4916912/1‑93 (MQ=255)
                                                                                  cGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTg       <  2:7516715/93‑1 (MQ=255)
                                                                                   gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGt      <  1:4396019/93‑1 (MQ=255)
                                                                                    aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt     >  1:9236576/1‑93 (MQ=255)
                                                                                    aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt     >  2:8260399/1‑93 (MQ=255)
                                                                                    aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt     >  2:8061877/1‑93 (MQ=255)
                                                                                     aGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGttt    <  2:7797246/93‑1 (MQ=255)
                                                                                      ggCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTa   <  1:3339388/93‑1 (MQ=255)
                                                                                       gcgcGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAc  <  2:4174306/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666539‑3666718

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