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% 917.0 / NA 288T64A (ACG→GCG) yhjAputative cytochrome C peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (151/136);  total (151/137)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666539‑3666718
                                                                                          |                                                                                         
 aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGcgtg                                                                                        <  1:720043/93‑1 (MQ=255)
 aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGcgtg                                                                                        <  2:3003914/93‑1 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       >  2:4998538/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       >  1:5536200/1‑93 (MQ=255)
  aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       >  2:3575498/1‑93 (MQ=255)
   aaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGc                                                                                      >  2:5604330/1‑93 (MQ=255)
   aaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGc                                                                                      >  2:2054539/1‑93 (MQ=255)
     gATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCaa                                                                                    <  1:5227029/93‑1 (MQ=255)
     gATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCaa                                                                                    >  1:2309229/1‑93 (MQ=255)
       tttATATCCAAGCTTAATGTCGTAATGCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  >  1:5418260/1‑93 (MQ=255)
         tataTCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAATGCGCGTGGCAATAGt                                                                                <  2:5903717/93‑1 (MQ=255)
          atatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                               >  1:1102506/1‑93 (MQ=255)
          atatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                               <  2:2002680/93‑1 (MQ=255)
             tctaaGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            >  2:4375468/4‑93 (MQ=255)
             tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            <  2:1026764/93‑1 (MQ=255)
             tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            >  2:3799195/1‑93 (MQ=255)
             tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            >  1:312526/1‑93 (MQ=255)
              ccAAGCTTAATGTCGTAATCCATCAACTGTCTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCAt                                                                           >  2:1438400/1‑93 (MQ=255)
               cAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATc                                                                          >  1:1578490/1‑93 (MQ=255)
                aaGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATcc                                                                         >  2:2323558/1‑93 (MQ=255)
                  gCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtt                                                                       <  2:3257555/93‑1 (MQ=255)
                   cTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttt                                                                      >  1:1874796/1‑93 (MQ=255)
                   cTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttt                                                                      <  1:4744163/93‑1 (MQ=255)
                    ttAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtttt                                                                     >  1:5778299/1‑93 (MQ=255)
                    ttAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtttt                                                                     >  2:1258417/1‑93 (MQ=255)
                     tAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttttt                                                                    >  2:2380286/1‑93 (MQ=255)
                     tAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCttttt                                                                    <  1:322480/93‑1 (MQ=255)
                      aaTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTc                                                                   <  2:621578/93‑1 (MQ=255)
                       aTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcg                                                                  >  1:4874796/1‑93 (MQ=255)
                       aTGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcg                                                                  >  2:4118487/1‑93 (MQ=255)
                         gTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTcgcg                                                                >  2:5524273/1‑93 (MQ=255)
                          tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               >  2:4638714/1‑93 (MQ=255)
                          tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               <  1:2095558/93‑1 (MQ=255)
                          tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               >  1:4597258/1‑93 (MQ=255)
                             tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                            >  2:2556830/1‑93 (MQ=255)
                              aaTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaa                                                           >  2:3463516/1‑93 (MQ=255)
                               aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  1:2425281/93‑1 (MQ=255)
                               aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  2:1105793/93‑1 (MQ=255)
                                 ccATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGc                                                        <  2:5625195/93‑1 (MQ=255)
                                  cATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGcc                                                       >  1:6284474/1‑93 (MQ=255)
                                   aTCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCt                                                      >  2:3506253/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  2:4408828/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  1:5887233/1‑93 (MQ=255)
                                    tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  1:1488584/1‑93 (MQ=255)
                                      aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   >  1:3329649/1‑93 (MQ=255)
                                      aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   >  2:1517925/1‑93 (MQ=255)
                                      aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   <  1:427185/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  <  2:5091842/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  <  1:3614264/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  <  2:5797496/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  <  2:288576/93‑1 (MQ=255)
                                       aCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaaa                                                  <  1:4028019/93‑1 (MQ=255)
                                        cTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAc                                                 <  2:3533763/93‑1 (MQ=255)
                                        cTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAc                                                 >  1:4530207/1‑93 (MQ=255)
                                        cTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAc                                                 >  2:3154127/1‑93 (MQ=255)
                                         tGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAAcc                                                <  2:5235181/93‑1 (MQ=255)
                                          gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               >  1:6159822/1‑93 (MQ=255)
                                            ttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtta                                             <  1:534935/93‑1 (MQ=255)
                                             tCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttat                                            <  1:3747405/93‑1 (MQ=255)
                                               gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTTGCGGAGAAAGCCTAAAACCttatta                                          >  2:7146/1‑93 (MQ=255)
                                               gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttatta                                          <  1:3066182/93‑1 (MQ=255)
                                               gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttatta                                          >  2:2635431/1‑93 (MQ=255)
                                               gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttatta                                          <  2:1275211/93‑1 (MQ=255)
                                                cgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         <  1:5325346/93‑1 (MQ=255)
                                                cgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         >  1:5133504/1‑93 (MQ=255)
                                                 gcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTattatt                                        <  2:2154823/93‑1 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                       >  2:3068161/1‑93 (MQ=255)
                                                  ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                       >  2:4273009/1‑93 (MQ=255)
                                                   cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                      >  2:3214526/1‑93 (MQ=255)
                                                   cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                      >  2:5424386/1‑93 (MQ=255)
                                                    aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                     >  1:2548338/1‑93 (MQ=255)
                                                      gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   >  2:4302202/1‑93 (MQ=255)
                                                      gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   <  1:4681647/93‑1 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                >  1:1010855/1‑93 (MQ=255)
                                                         tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                <  1:1693675/93‑1 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              <  1:1264169/93‑1 (MQ=255)
                                                           taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              <  2:4718894/93‑1 (MQ=255)
                                                             ataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagc                            >  1:4436661/1‑93 (MQ=255)
                                                               aataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagcag                          <  2:3377761/93‑1 (MQ=255)
                                                                ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         >  2:1498104/1‑93 (MQ=255)
                                                                ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  1:977212/93‑1 (MQ=255)
                                                                  aGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATg                       <  1:3434484/93‑1 (MQ=255)
                                                                      gggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTg                   >  1:1456271/1‑93 (MQ=255)
                                                                       ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                  <  2:1703157/93‑1 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAAc                <  2:2938656/93‑1 (MQ=255)
                                                                         tAATTCTGCCGAAGGCGCGTGGCAATAGACGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAAc                >  1:2504091/1‑93 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              <  2:5353035/93‑1 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              >  1:3171653/1‑93 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              >  2:4936634/1‑93 (MQ=255)
                                                                           aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              <  1:1769721/93‑1 (MQ=255)
                                                                            ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             <  1:4107240/93‑1 (MQ=255)
                                                                            ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             >  1:2662219/1‑93 (MQ=255)
                                                                               tGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGc          >  2:4811148/1‑93 (MQ=255)
                                                                               tGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGc          <  1:993540/93‑1 (MQ=255)
                                                                               tGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGc          <  1:2938364/93‑1 (MQ=255)
                                                                               tGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGc          <  1:2381693/93‑1 (MQ=255)
                                                                                   gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGt      >  2:5244205/1‑93 (MQ=255)
                                                                                    aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt     >  2:5737642/1‑93 (MQ=255)
                                                                                     aGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGttt    >  2:4894365/1‑93 (MQ=255)
                                                                                      ggCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTa   >  1:1939633/1‑93 (MQ=255)
                                                                                          |                                                                                         
TAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666539‑3666718

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