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.6% 861.3 / ‑6.3 236intergenic (‑23/+9)b3530/yhjMputative oxidoreductase subunit/putative endoglucanase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (92/143);  minor base A (1/0);  total (93/143)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 262 total aligned reads displayed.

TCGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATC  >  NC_000913/3686683‑3686861
                                                                                             |                                                                                      
tCGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCg                                                                                         <  1:5340582/93‑1 (MQ=255)
tCGCGTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCg                                                                                         <  1:1239296/93‑1 (MQ=255)
 cGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgt                                                                                        <  1:3970378/93‑2 (MQ=255)
  gCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGGCAGACCGAGGGAAAGCgtg                                                                                       >  1:257688/1‑91 (MQ=255)
   cTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtga                                                                                      >  1:1134912/1‑90 (MQ=255)
   cTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtga                                                                                      >  2:5900262/1‑90 (MQ=255)
   cTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtga                                                                                      <  2:551917/93‑4 (MQ=255)
    ttCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCgtgaa                                                                                     >  1:5494803/1‑89 (MQ=255)
      cGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAata                                                                                   <  2:1577582/93‑1 (MQ=255)
          tAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGtt                                                                               >  2:5186093/1‑93 (MQ=255)
          tAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGtt                                                                               >  2:2744169/1‑93 (MQ=255)
          tAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGtt                                                                               <  2:4971318/93‑1 (MQ=255)
          tAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGACGGAAAGCGTGAATATGtt                                                                               >  1:4496409/1‑93 (MQ=255)
           aaCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGttt                                                                              >  1:2480797/1‑93 (MQ=255)
           aaCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGttt                                                                              >  2:4723044/1‑93 (MQ=255)
           aaCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGttt                                                                              <  1:3326580/93‑1 (MQ=255)
            aCCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTa                                                                             >  2:2651748/1‑93 (MQ=255)
             ccGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAg                                                                            <  1:1030359/93‑1 (MQ=255)
             ccGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAg                                                                            >  1:792916/1‑93 (MQ=255)
              cGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAgt                                                                           <  1:2498316/93‑1 (MQ=255)
              cGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAgt                                                                           <  2:787086/93‑1 (MQ=255)
                aaCTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAgtgt                                                                         >  1:4027200/1‑93 (MQ=255)
                aaCTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAgtgt                                                                         >  1:1223356/1‑93 (MQ=255)
                   ttGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGaa                                                                      >  1:5884295/1‑93 (MQ=255)
                   ctGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGaa                                                                      <  1:5181957/92‑1 (MQ=255)
                    tGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAt                                                                     <  1:1403950/93‑1 (MQ=255)
                     gCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAtt                                                                    <  2:3998477/93‑1 (MQ=255)
                      cTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAttt                                                                   >  2:3428764/1‑93 (MQ=255)
                      cTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAttt                                                                   <  1:3606638/93‑1 (MQ=255)
                      cTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAttt                                                                   >  1:3785864/1‑93 (MQ=255)
                      cTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAAttt                                                                   <  1:1887806/93‑1 (MQ=255)
                       tCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTg                                                                  <  2:3112797/93‑1 (MQ=255)
                       tCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTg                                                                  <  2:4126308/93‑1 (MQ=255)
                       cccAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTg                                                                  <  1:242445/92‑1 (MQ=255)
                        ccAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTgc                                                                 <  1:5192651/93‑1 (MQ=255)
                        ccAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTgc                                                                 >  2:4458254/1‑93 (MQ=255)
                           gCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCa                                                              <  1:3353997/93‑1 (MQ=255)
                           gCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCa                                                              >  2:5607609/1‑93 (MQ=255)
                           gCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCa                                                              <  2:3565817/93‑1 (MQ=255)
                            cAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCAt                                                             <  2:6058910/93‑1 (MQ=255)
                             aaCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCAtt                                                            >  2:575506/1‑93 (MQ=255)
                             aaCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGAGTGAATTTGCGCttt                                                            >  1:4994718/1‑93 (MQ=255)
                              aCTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTc                                                           >  1:6056299/1‑93 (MQ=255)
                               cTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTcc                                                          <  2:5136976/93‑1 (MQ=255)
                               cTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTcc                                                          <  1:2734725/93‑1 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCt                                                         <  2:5245541/93‑1 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCt                                                         <  2:2421703/93‑1 (MQ=255)
                                tGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCt                                                         <  2:1335971/93‑1 (MQ=255)
                                  ttGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTgg                                                       >  1:3027924/1‑93 (MQ=255)
                                    gCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGCTTAGTGTGAATTTGCGCATTCCTGGcc                                                     >  1:4884555/1‑93 (MQ=255)
                                      tGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGcccc                                                   >  2:2550480/1‑93 (MQ=255)
                                      tGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGcccc                                                   >  2:2033654/1‑93 (MQ=255)
                                       gAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCa                                                  <  2:932427/93‑1 (MQ=255)
                                       gAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCa                                                  >  1:202774/1‑93 (MQ=255)
                                       gAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCa                                                  <  2:5367710/93‑1 (MQ=255)
                                        aGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAg                                                 >  1:3192801/1‑93 (MQ=255)
                                           ggTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCa                                              <  2:4314985/93‑1 (MQ=255)
                                            gTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTACTGGCCCCAGTCAg                                             >  2:1068074/1‑93 (MQ=255)
                                              tGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGt                                           <  1:5555944/93‑1 (MQ=255)
                                               ggTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGCTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGta                                          >  2:4014506/1‑93 (MQ=255)
                                                gTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaa                                         <  1:2962080/93‑1 (MQ=255)
                                                  gctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaata                                       >  1:753748/1‑93 (MQ=255)
                                                  gctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaata                                       <  1:6141396/93‑1 (MQ=255)
                                                  gctgcCGCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATGCCTGGCCCCAGTCAGGtaata                                       <  2:1419409/93‑1 (MQ=255)
                                                   ctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaataa                                      >  1:56808/1‑93 (MQ=255)
                                                   ctgcCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGtaataa                                      <  2:373624/93‑1 (MQ=255)
                                                      ccTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTc                                   <  1:1380172/93‑1 (MQ=255)
                                                      ccTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTc                                   <  1:4641011/93‑1 (MQ=255)
                                                       cTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCa                                  <  1:3263195/93‑1 (MQ=255)
                                                        tCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCAc                                 >  1:5091793/1‑93 (MQ=255)
                                                          gACCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCt                               <  1:407346/93‑1 (MQ=255)
                                                           aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtt                              <  1:1533400/93‑1 (MQ=255)
                                                           aCCATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCACTCCTGGCCCCAGTCAGGTTATAACTCACCtt                              <  2:5042254/93‑1 (MQ=255)
                                                             cATCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCtttt                            <  1:306117/93‑1 (MQ=255)
                                                               tCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGt                          <  1:6260819/93‑1 (MQ=255)
                                                               tCGGCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGt                          <  1:2320751/93‑1 (MQ=255)
                                                                 ggCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCg                        >  2:3783260/1‑93 (MQ=255)
                                                                 ggCATGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCg                        <  2:2142301/93‑1 (MQ=255)
                                                                   gaTGACGCCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCgag                      <  2:745771/92‑1 (MQ=255)
                                                                     tGACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGaa                    >  2:3349051/1‑93 (MQ=255)
                                                                      gACGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAg                   <  1:6085823/93‑1 (MQ=255)
                                                                       gcGGCCAGACCGAGGGATAGCGTAAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGc                  <  2:1404812/92‑1 (MQ=255)
                                                                       aCGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGc                  <  1:2374251/93‑1 (MQ=255)
                                                                       aCGGCCAGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGc                  <  1:3101165/93‑1 (MQ=255)
                                                                             aGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAAc            <  1:4153/93‑1 (MQ=255)
                                                                             aGACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAAc            <  1:4270495/93‑1 (MQ=255)
                                                                              gACCGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACg           >  2:435744/1‑93 (MQ=255)
                                                                                ccGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgt         >  2:165962/1‑93 (MQ=255)
                                                                                ccGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgt         >  2:2358361/1‑93 (MQ=255)
                                                                                ccGAGAGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgt         <  1:5723489/93‑1 (MQ=255)
                                                                                 cGAGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGgtg        <  1:3650750/93‑1 (MQ=255)
                                                                                   aGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGtt      <  2:4608377/93‑1 (MQ=255)
                                                                                   aGGGAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGtt      >  1:3416367/1‑93 (MQ=255)
                                                                                    gggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTg     <  2:2867932/93‑1 (MQ=255)
                                                                                     ggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGa    <  2:2012572/93‑1 (MQ=255)
                                                                                     ggAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGa    <  1:1628056/93‑1 (MQ=255)
                                                                                      gAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAt   <  2:3037573/93‑1 (MQ=255)
                                                                                      gAAAGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGAt   >  1:2171308/1‑93 (MQ=255)
                                                                                       aaaGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATc  >  1:1667884/1‑93 (MQ=255)
                                                                                       aaaGCGTGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATc  >  1:4507340/1‑93 (MQ=255)
                                                                                             |                                                                                      
TCGCTTCGCCTAACCGAACTTGCTCCAGCAACTGTTGCTGAGCGGTTGGTGCTGCCTCGACCATCGGCATGACGGCCAGACCGAGGGAAAGCGGAATATGTTTAGTGTGAATTTGCGCATTCCTGGCCCCAGTCAGGTAATAACTCACCTTTTGTCGAGAAGCGGAAACGGTGTTGATC  >  NC_000913/3686683‑3686861

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