Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,667,169 T→C intergenic (‑351/‑53) yhjA ← / → treF putative cytochrome C peroxidase/cytoplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,667,1690TC84.6% 591.8 / 98.2 260intergenic (‑351/‑53)yhjA/treFputative cytochrome C peroxidase/cytoplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base C (105/115);  total (124/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 263 total aligned reads displayed.

AATTCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGAC  >  NC_000913/3667077‑3667257
                                                                                            |                                                                                        
aaTTCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACt                                                                                          <  1:3099345/93‑1 (MQ=255)
 aTTCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTa                                                                                         <  1:2684901/93‑1 (MQ=255)
  ttCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTAc                                                                                        >  2:5937672/1‑93 (MQ=255)
   tCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGaccacc                                                                                       >  1:4211820/1‑93 (MQ=255)
    ccATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCAccc                                                                                      >  2:6393573/1‑93 (MQ=255)
    ccATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCAccc                                                                                      <  2:2546306/93‑1 (MQ=255)
     cATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCt                                                                                     >  2:4708990/1‑93 (MQ=255)
      aTTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCtt                                                                                    >  1:446390/1‑93 (MQ=255)
       ttGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTa                                                                                   >  2:3715630/1‑93 (MQ=255)
        tGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTaa                                                                                  <  1:2138215/93‑1 (MQ=255)
        tGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTaa                                                                                  <  1:3639793/93‑1 (MQ=255)
         gATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTaaa                                                                                 <  1:6423107/93‑1 (MQ=255)
         gATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTaaa                                                                                 <  1:633692/93‑1 (MQ=255)
          atatTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAg                                                                                <  1:4134126/93‑1 (MQ=255)
          atatTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAg                                                                                <  1:1386709/93‑1 (MQ=255)
           tatTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGa                                                                               >  1:5408021/1‑93 (MQ=255)
              tAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGaaaa                                                                            <  1:277974/93‑1 (MQ=255)
              tAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGaaaa                                                                            <  2:6468662/93‑1 (MQ=255)
                 tGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGcc                                                                         >  1:6643497/1‑93 (MQ=255)
                 tGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGcc                                                                         >  2:5368539/1‑93 (MQ=255)
                  gCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCg                                                                        >  2:664491/1‑93 (MQ=255)
                   cGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGa                                                                       >  1:1608570/1‑93 (MQ=255)
                     ttATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATa                                                                     >  1:3018343/1‑93 (MQ=255)
                     ttATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATa                                                                     >  2:1087541/1‑93 (MQ=255)
                      tataTTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATaa                                                                    >  2:5801935/1‑93 (MQ=255)
                      tataTTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATaa                                                                    >  2:4820281/1‑93 (MQ=255)
                          ttttttGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGGGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAg                                                                >  1:4934121/1‑93 (MQ=255)
                          ttttttGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAg                                                                >  2:3950324/1‑93 (MQ=255)
                          ttttttGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAg                                                                >  2:2657351/1‑93 (MQ=255)
                           tttttGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGc                                                               <  2:4146406/93‑1 (MQ=255)
                             tttGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGcga                                                             <  1:3776107/93‑1 (MQ=255)
                                gACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGcgacga                                                          <  1:7209243/93‑1 (MQ=255)
                                gACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGcgacga                                                          >  2:736441/1‑93 (MQ=255)
                                 aCAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGaa                                                         <  1:32867/93‑1 (MQ=255)
                                   aGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAAtt                                                       >  2:5259524/1‑93 (MQ=255)
                                     gTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTc                                                     >  1:6945381/1‑93 (MQ=255)
                                      ttAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCg                                                    <  1:5865454/93‑1 (MQ=255)
                                        aataaACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGa                                                  >  2:708821/1‑93 (MQ=255)
                                        aataaACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGa                                                  >  2:7149515/1‑93 (MQ=255)
                                         ataaACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAg                                                 <  2:4059021/93‑1 (MQ=255)
                                         ataaACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAg                                                 <  2:4229361/93‑1 (MQ=255)
                                         ataaACGTGATCCACCGCACGCTTTGTCGCCCACCAGCCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAg                                                 <  2:1504936/93‑1 (MQ=255)
                                            aaCGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGtt                                              >  2:6094229/1‑93 (MQ=255)
                                             aCGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTg                                             >  2:6278230/1‑93 (MQ=255)
                                              cGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTgg                                            <  1:4710632/93‑1 (MQ=255)
                                              cGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTgg                                            <  1:1044913/93‑1 (MQ=255)
                                               gTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGa                                           <  2:3431394/93‑1 (MQ=255)
                                                tGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGAt                                          <  2:1332004/93‑1 (MQ=255)
                                                  aTCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATcc                                        >  1:857397/1‑93 (MQ=255)
                                                  aTCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATcc                                        <  1:6127663/93‑1 (MQ=255)
                                                   tCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCt                                       >  1:7203232/1‑93 (MQ=255)
                                                   tCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCt                                       >  1:3357849/1‑93 (MQ=255)
                                                   tCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCt                                       <  2:1823874/93‑1 (MQ=255)
                                                    ccACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCtt                                      <  1:2738583/93‑1 (MQ=255)
                                                     cACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTa                                     <  1:1003755/93‑1 (MQ=255)
                                                     cACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTa                                     <  1:3591051/93‑1 (MQ=255)
                                                     cACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTa                                     >  1:5663922/1‑93 (MQ=255)
                                                      aCCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTAt                                    >  1:3962660/1‑93 (MQ=255)
                                                      aCCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTAt                                    <  2:1911178/93‑1 (MQ=255)
                                                        cGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGc                                  >  2:1747321/1‑93 (MQ=255)
                                                        cGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGc                                  >  1:17720/1‑93 (MQ=255)
                                                        cGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGc                                  <  2:7475387/93‑1 (MQ=255)
                                                        cGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGc                                  >  1:3849703/1‑93 (MQ=255)
                                                        cGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGc                                  <  2:4013424/93‑1 (MQ=255)
                                                          cACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTc                                >  1:6484026/1‑93 (MQ=255)
                                                          cACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTc                                >  2:2320373/1‑93 (MQ=255)
                                                          cACGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTc                                <  1:658503/93‑1 (MQ=255)
                                                           aCGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCa                               >  2:3213570/1‑93 (MQ=255)
                                                           aCGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCa                               >  1:531460/1‑93 (MQ=255)
                                                            cGCTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCaa                              <  2:4532942/93‑1 (MQ=255)
                                                              cTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATc                            >  1:3370749/1‑93 (MQ=255)
                                                              cTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATc                            >  2:4656347/1‑93 (MQ=255)
                                                              cTTTGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATc                            >  1:4616520/1‑93 (MQ=255)
                                                               tttGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCa                           >  1:5498877/1‑93 (MQ=255)
                                                               tttGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCa                           <  1:6689535/93‑1 (MQ=255)
                                                                 tGTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGa                         <  1:1433363/93‑1 (MQ=255)
                                                                  gTCGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGaa                        >  2:1562868/1‑93 (MQ=255)
                                                                    cGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGaaaa                      <  2:3047798/93‑1 (MQ=255)
                                                                    cGCCCACCAGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGaaaa                      <  2:6606538/93‑1 (MQ=255)
                                                                     gCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAAt                     <  2:768083/93‑1 (MQ=255)
                                                                       ccaccaGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTc                   <  2:5137111/93‑1 (MQ=255)
                                                                        caccaGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCa                  <  1:5614986/93‑1 (MQ=255)
                                                                        caccaGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCa                  >  2:6676584/1‑93 (MQ=255)
                                                                         accaGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCaa                 >  2:7310300/1‑93 (MQ=255)
                                                                           caGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCaaaa               <  2:5620048/93‑1 (MQ=255)
                                                                           caGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCaaaa               <  1:2488415/93‑1 (MQ=255)
                                                                           caGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCaaaa               <  2:4883798/93‑1 (MQ=255)
                                                                           caGGCGGAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCaaaa               <  2:5036885/93‑1 (MQ=255)
                                                                              gCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAAccc            <  2:5303067/93‑1 (MQ=255)
                                                                                ggAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTa          <  2:5438748/93‑1 (MQ=255)
                                                                                 gAGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTaa         >  2:2944041/1‑93 (MQ=255)
                                                                                  aGCGAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAAt        >  2:1044469/1‑93 (MQ=255)
                                                                                     gAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCa     >  1:5991089/1‑93 (MQ=255)
                                                                                     gAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCa     >  1:3491823/1‑93 (MQ=255)
                                                                                     gAATGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCa     <  1:4353198/93‑1 (MQ=255)
                                                                                      aaTGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAg    >  1:3253112/1‑93 (MQ=255)
                                                                                      aaTGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAg    >  1:5774285/1‑93 (MQ=255)
                                                                                        tGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGAc  <  1:2502350/93‑1 (MQ=255)
                                                                                        tGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGAc  >  1:1850168/1‑93 (MQ=255)
                                                                                        tGACCACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGAc  >  2:4419227/1‑93 (MQ=255)
                                                                                            |                                                                                        
AATTCCATTGATATTAGTGCGTTATATTTTTTGACAGGTTAATAAACGTGATCCACCGCACGCTTTGTCGCCCACCAGGCGGAGCGAATGACTACCCTTAAAGAAAAGCCGATAATTAGCGACGAATTTCGGAGGTTGGATCCTTATGCTCAATCAGAAAATTCAAAACCCTAATCCAGAC  >  NC_000913/3667077‑3667257

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