Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,686,775 +T intergenic (‑23/+9) b3530 ← / ← yhjM putative oxidoreductase subunit/putative endoglucanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,686,7751.T99.8% 1471.9 / NA 411intergenic (‑23/+9)b3530/yhjMputative oxidoreductase subunit/putative endoglucanase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base T (160/249);  minor base A (1/0);  total (162/249)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 449 total aligned reads displayed.

TCGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATC  >  NC_000913/3686683‑3686861
                                                                                             |                                                                                      
 cGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgt                                                                                        <  1:6599507/93‑2 (MQ=255)
 cGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgt                                                                                        <  1:8319722/93‑2 (MQ=255)
  gCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtg                                                                                       <  1:4669467/93‑3 (MQ=255)
   cTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtga                                                                                      <  1:6080283/93‑4 (MQ=255)
    ttCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtgaa                                                                                     >  1:3228771/1‑89 (MQ=255)
      cGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGGGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAata                                                                                   >  1:8409444/1‑93 (MQ=255)
      cGCCTAACCGAACGTGCTCCAGCAGCTGTTCGTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAata                                                                                   <  1:3700485/93‑1 (MQ=255)
       gCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAatat                                                                                  <  1:1425884/93‑1 (MQ=255)
       gCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAatat                                                                                  >  1:8185357/1‑93 (MQ=255)
       gCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAatat                                                                                  <  2:6640688/93‑1 (MQ=255)
         cTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGCCGGCCAGACCGAGGGAAAGCGTGAATATGt                                                                                >  1:5045976/1‑93 (MQ=255)
         cTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGt                                                                                >  2:875682/1‑93 (MQ=255)
         cTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGt                                                                                >  1:902983/1‑93 (MQ=255)
         cTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGt                                                                                >  2:5346147/1‑93 (MQ=255)
         cTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGt                                                                                >  2:8455070/1‑93 (MQ=255)
          tAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGtt                                                                               <  1:7344930/93‑1 (MQ=255)
           aaCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGttt                                                                              <  1:1714649/93‑1 (MQ=255)
            aCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTa                                                                             <  1:4496932/93‑1 (MQ=255)
             ccGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAg                                                                            >  1:4310969/1‑93 (MQ=255)
              cGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAgt                                                                           >  2:818653/1‑93 (MQ=255)
                  cTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGa                                                                       >  2:1856688/1‑93 (MQ=255)
                   ttGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGaa                                                                      <  2:5761045/93‑1 (MQ=255)
                   ttGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGaa                                                                      <  1:7093490/93‑1 (MQ=255)
                     gCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAtt                                                                    <  2:6555370/93‑1 (MQ=255)
                      cTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAttt                                                                   <  2:5934824/93‑1 (MQ=255)
                        ccAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACTGCCAGACCGAGGGAAAGCGTGAATATGTTTAGCGTGAATTTgc                                                                 >  2:7300563/1‑93 (MQ=255)
                           gCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCa                                                              >  2:3987672/1‑93 (MQ=255)
                           gCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCa                                                              >  1:4872934/1‑93 (MQ=255)
                            cAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCAt                                                             >  2:7069966/1‑93 (MQ=255)
                             aaCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCAtt                                                            >  1:8059815/1‑93 (MQ=255)
                              aCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTc                                                           <  2:2915377/93‑1 (MQ=255)
                              aCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTc                                                           >  1:984947/1‑93 (MQ=255)
                               cTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTcc                                                          >  1:5243429/1‑93 (MQ=255)
                               cTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTcc                                                          <  2:4710537/93‑1 (MQ=255)
                               cTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTcc                                                          <  2:4919054/93‑1 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTctg                                                         >  2:2713605/1‑91 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCt                                                         >  2:753357/1‑93 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAAAATGTTTAGTGTGAATTTGCGCATTCCt                                                         >  2:3654218/1‑93 (MQ=255)
                                tGTTGCTGAGCGGGTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCt                                                         >  1:199666/1‑93 (MQ=255)
                                  ttGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTgg                                                       <  2:4482843/93‑1 (MQ=255)
                                  gtGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTgg                                                       <  1:7604473/92‑1 (MQ=255)
                                    gCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGcc                                                     <  2:8075735/93‑1 (MQ=255)
                                     cTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGccc                                                    <  1:6779531/93‑1 (MQ=255)
                                       gAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCa                                                  >  2:8524669/1‑93 (MQ=255)
                                        ggCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAg                                                 <  1:1051744/92‑1 (MQ=255)
                                        aggggtgggtgCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAg                                                 <  2:6454347/90‑1 (MQ=255)
                                        aGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAg                                                 <  2:6725661/93‑1 (MQ=255)
                                         gCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGt                                                <  1:1868193/93‑1 (MQ=255)
                                           ggatGGTGATGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTAGCGCATTCCTGGCCCCAGTCa                                              <  2:7831381/90‑1 (MQ=255)
                                           ggTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCa                                              <  2:7235620/93‑1 (MQ=255)
                                            gTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAg                                             <  2:5483126/93‑1 (MQ=255)
                                             ttGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAgg                                            <  1:8127031/93‑1 (MQ=255)
                                               ggTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGta                                          <  1:3938292/93‑1 (MQ=255)
                                                 tgctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaat                                        <  2:4262823/93‑1 (MQ=255)
                                                  gctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaata                                       >  1:7780337/1‑93 (MQ=255)
                                                  gctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaata                                       <  2:8646232/93‑1 (MQ=255)
                                                   ctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaataa                                      >  2:7116807/1‑93 (MQ=255)
                                                   ctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaataa                                      <  2:314362/93‑1 (MQ=255)
                                                    tgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAAc                                     >  2:6163971/1‑93 (MQ=255)
                                                        tCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCAc                                 >  2:2526366/1‑93 (MQ=255)
                                                        tCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCAc                                 <  2:635334/93‑1 (MQ=255)
                                                           aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtt                              >  1:6322754/1‑93 (MQ=255)
                                                           aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtt                              >  1:2887695/1‑93 (MQ=255)
                                                           aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtt                              <  1:8620230/93‑1 (MQ=255)
                                                             cATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtttt                            <  2:2277811/93‑1 (MQ=255)
                                                              aTCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTg                           <  2:7175896/93‑1 (MQ=255)
                                                              aTCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTg                           >  2:7654121/1‑93 (MQ=255)
                                                               tCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGt                          >  1:6980288/1‑93 (MQ=255)
                                                               tCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGt                          <  2:2006567/93‑1 (MQ=255)
                                                                 ggCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCg                        >  2:7965315/1‑93 (MQ=255)
                                                                  gCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCga                       >  1:3288090/1‑93 (MQ=255)
                                                                  gCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCga                       >  2:1266946/1‑93 (MQ=255)
                                                                  gCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCga                       >  1:148561/1‑93 (MQ=255)
                                                                    aTGACGGCGAGACAGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCgaga                     <  1:5353784/93‑1 (MQ=255)
                                                                    aTGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCgaga                     <  1:236591/93‑1 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    >  1:8710089/1‑93 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    >  1:2183595/1‑93 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    <  1:131947/93‑1 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    <  2:3412933/93‑1 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    <  1:5893693/93‑1 (MQ=255)
                                                                      gACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAg                   <  1:3062980/93‑1 (MQ=255)
                                                                       aCGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGc                  <  1:4954468/93‑1 (MQ=255)
                                                                         ggCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCgg                <  2:6967789/93‑1 (MQ=255)
                                                                          gCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGa               <  1:3356022/93‑1 (MQ=255)
                                                                           ccAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGaa              >  1:6395207/1‑93 (MQ=255)
                                                                            cAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGaaa             <  2:4797042/93‑1 (MQ=255)
                                                                            cAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGaaa             >  2:1502528/1‑93 (MQ=255)
                                                                               aCCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACgg          <  2:3609477/93‑1 (MQ=255)
                                                                               aCCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACgg          <  1:3820075/93‑1 (MQ=255)
                                                                               aCCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACgg          <  1:4758711/93‑1 (MQ=255)
                                                                                ccGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgt         <  2:457323/93‑1 (MQ=255)
                                                                                 cGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgtg        >  1:4577423/1‑93 (MQ=255)
                                                                                  gAGGGAAGGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgtgt       <  2:3742017/93‑1 (MQ=255)
                                                                                  gAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgtgt       <  2:1155059/93‑1 (MQ=255)
                                                                                  gAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgtgt       <  1:2535029/93‑1 (MQ=255)
                                                                                   aGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGtt      <  1:5516980/93‑1 (MQ=255)
                                                                                    gggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTg     <  2:3252687/93‑1 (MQ=255)
                                                                                    gggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTg     >  2:5094630/1‑93 (MQ=255)
                                                                                     ggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGa    <  1:1192774/93‑1 (MQ=255)
                                                                                       aaaGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATc  >  1:269767/1‑93 (MQ=255)
                                                                                             |                                                                                      
TCGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATC  >  NC_000913/3686683‑3686861

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