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,6290TC98.1% 819.3 / 1.1 263T64A (ACG→GCG) yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base C (140/117);  total (141/122)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 280 total aligned reads displayed.

TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTA  >  NC_000913/3666539‑3666717
                                                                                          |                                                                                        
tAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCgcgt                                                                                        <  2:1674917/93‑1 (MQ=255)
tAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCgcgt                                                                                        >  1:3972366/1‑93 (MQ=255)
 aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGcgtg                                                                                       >  1:4495497/1‑93 (MQ=255)
 aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGACGAAGGCGcgtg                                                                                       >  2:573207/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                      >  2:3450417/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                      >  1:1608164/1‑93 (MQ=255)
   aaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGc                                                                                     <  2:1907871/93‑1 (MQ=255)
   aaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGc                                                                                     <  2:2237207/93‑1 (MQ=255)
     gATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCaa                                                                                   >  2:3908551/1‑93 (MQ=255)
      aTTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAAt                                                                                  <  2:1535819/93‑1 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                 <  1:1003295/93‑1 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                 <  1:2461486/93‑1 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                 <  1:3685196/93‑1 (MQ=255)
         tataTCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGt                                                                               <  1:3488939/93‑1 (MQ=255)
          atatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                              >  1:4749871/1‑93 (MQ=255)
           tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCg                                                                             >  1:1841580/1‑93 (MQ=255)
           tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCg                                                                             >  1:735844/1‑93 (MQ=255)
            atCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGc                                                                            <  1:2963576/93‑1 (MQ=255)
              ccAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCAt                                                                          >  1:1420821/1‑93 (MQ=255)
              ccAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCAt                                                                          <  2:3067789/93‑1 (MQ=255)
                aaGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATcc                                                                        <  2:2009042/93‑1 (MQ=255)
                aaGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATcc                                                                        <  2:1602023/93‑1 (MQ=255)
                 aGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCt                                                                       >  1:1151398/1‑93 (MQ=255)
                  gCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtt                                                                      >  2:464662/1‑93 (MQ=255)
                  gCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtt                                                                      <  1:2231575/93‑1 (MQ=255)
                    ttAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtttt                                                                    >  2:435516/1‑93 (MQ=255)
                      aaTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTc                                                                  <  1:2850558/93‑1 (MQ=255)
                      aaTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTc                                                                  <  2:2978314/93‑1 (MQ=255)
                        tGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcgc                                                                >  2:314065/1‑93 (MQ=255)
                          tCGTAATCCATCAACTGTTTCCCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                              <  2:1400338/93‑1 (MQ=255)
                            gTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgag                                                            >  1:4353548/1‑93 (MQ=255)
                            gTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgag                                                            <  1:996944/93‑1 (MQ=255)
                             tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                           >  1:1659750/1‑93 (MQ=255)
                             tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                           >  1:3807441/1‑93 (MQ=255)
                             tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                           <  2:4037958/93‑1 (MQ=255)
                                 ccATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAGTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGc                                                       >  2:3837868/1‑93 (MQ=255)
                                 ccATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGc                                                       >  1:2048640/1‑93 (MQ=255)
                                  cATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGcc                                                      >  2:45179/1‑93 (MQ=255)
                                  cATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGcc                                                      <  2:4153640/93‑1 (MQ=255)
                                   aTCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCt                                                     >  1:1944322/1‑93 (MQ=255)
                                   aTCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCt                                                     >  2:319743/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                    >  1:3914251/1‑93 (MQ=255)
                                     cAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaa                                                   >  2:3547890/1‑93 (MQ=255)
                                      aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                  <  2:3009568/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                 >  2:186175/1‑93 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                 <  2:932195/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                 >  1:3734258/1‑93 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                 >  1:1432400/1‑93 (MQ=255)
                                        cTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAc                                                >  1:2597744/1‑93 (MQ=255)
                                         tGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAcc                                               <  2:1509361/93‑1 (MQ=255)
                                         tGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAcc                                               <  2:653922/93‑1 (MQ=255)
                                          gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                              >  1:4279249/1‑93 (MQ=255)
                                           tttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtt                                             <  2:1450523/93‑1 (MQ=255)
                                             tCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttat                                           <  1:481736/93‑1 (MQ=255)
                                               gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttatta                                         <  2:1648335/93‑1 (MQ=255)
                                                cgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                        >  2:208520/1‑93 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                       >  1:4066245/1‑93 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGTTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                      >  1:815337/1‑93 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                      >  2:1076718/1‑93 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                      <  1:3554908/93‑1 (MQ=255)
                                                   cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                     >  1:690171/1‑93 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                    <  2:724827/93‑1 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                    >  2:473538/1‑93 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCCTATTATTTTc                                    >  1:3618047/1‑93 (MQ=255)
                                                     ggAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCa                                   >  2:2330206/1‑93 (MQ=255)
                                                      gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                  >  1:3148253/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                               <  1:4250736/93‑1 (MQ=255)
                                                          ataataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAc                              >  1:2006687/1‑93 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                             <  1:90446/93‑1 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                             >  2:3314674/1‑93 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                             <  1:1789538/93‑1 (MQ=255)
                                                             ataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagc                           >  1:3856178/1‑93 (MQ=255)
                                                              taataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagca                          >  1:4645713/1‑93 (MQ=255)
                                                                ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                        <  2:2710035/93‑1 (MQ=255)
                                                                 taGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGAt                       <  2:2760768/93‑1 (MQ=255)
                                                                   ggCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGc                     >  1:32666/1‑93 (MQ=255)
                                                                    gCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGcc                    <  2:420857/93‑1 (MQ=255)
                                                                    gCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGcc                    >  1:420042/1‑93 (MQ=255)
                                                                       ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                 >  1:595349/1‑93 (MQ=255)
                                                                       ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                 >  1:1196175/1‑93 (MQ=255)
                                                                       ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                 >  1:3447391/1‑93 (MQ=255)
                                                                        gTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGaa                <  1:191700/93‑1 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCGATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGACGCCTGAAc               >  1:1134436/1‑93 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAAc               >  1:2574256/1‑93 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCAGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAAc               <  2:2908664/93‑1 (MQ=255)
                                                                          aaTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACa              >  2:1330120/1‑93 (MQ=255)
                                                                          aaTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACa              >  2:776556/1‑93 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt             <  2:3776908/93‑1 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt             >  2:246590/1‑93 (MQ=255)
                                                                            ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc            <  1:3378065/93‑1 (MQ=255)
                                                                             tCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCg           <  2:4474524/93‑1 (MQ=255)
                                                                                tccGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGcc        <  1:2513224/92‑1 (MQ=255)
                                                                                gCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGcc        <  2:258090/93‑1 (MQ=255)
                                                                                  cGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTg      <  1:4413785/93‑1 (MQ=255)
                                                                                   gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTGACTGACAGCAGATGCCTGAACATCGGCCTGt     >  2:3949997/1‑93 (MQ=255)
                                                                                   gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGt     <  2:2554014/93‑1 (MQ=255)
                                                                                    aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt    <  2:1808874/93‑1 (MQ=255)
                                                                                     aGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGttt   <  2:1112179/93‑1 (MQ=255)
                                                                                     aGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGttt   >  2:3377216/1‑93 (MQ=255)
                                                                                      ggCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTa  <  2:2281275/93‑1 (MQ=255)
                                                                                          |                                                                                        
TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTA  >  NC_000913/3666539‑3666717

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