Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,925 +A coding (35/1431 nt) yaaJ ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,9221.A98.4% 197.1 / ‑3.8 62coding (38/1431 nt)yaaJputative transporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (31/30);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7786‑8056
                                                                                                                                         |                                                                                                                                       
tAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaa                                                                                                                                      <  1:837623/139‑1 (MQ=255)
  aaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCa                                                                                                                                        <  2:1124814/136‑1 (MQ=255)
  aaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCa                                                                                                                                        >  1:1124814/1‑136 (MQ=255)
   aCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaag                                                                                                                                   >  1:176552/1‑137 (MQ=255)
      gCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaaggac                                                                                                                                <  2:567306/139‑6 (MQ=255)
       ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:1381530/137‑5 (MQ=255)
       ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  1:1381530/1‑133 (MQ=255)
       acaggCTGTGGATGAATGCTATTTTTAAGACTGTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:235408/136‑1 (MQ=255)
        cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               >  2:639160/1‑139 (MQ=255)
         tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:485114/137‑1 (MQ=255)
         tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:485114/1‑137 (MQ=255)
                      aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATGCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaaggac                                                                                                                                >  1:1651252/1‑118 (MQ=255)
                      aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGgag                                                                                                                 >  2:562530/1‑139 (MQ=255)
                      aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCAGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaaggac                                                                                                                                <  2:1651252/123‑6 (MQ=255)
                      aTGCTATTTTTAAGACTTTTGCAAAACTGGCAAATGTAGCGAAACTGGACAAGTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCCACGAGCAGGTAAATCATTGCCGATCCCCAAAAAGACGCTGATAATGaactgt                                                                                                                 >  1:235408/1‑135 (MQ=255)
                         cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGaaa                                                                                                              >  1:235249/1‑139 (MQ=255)
                         cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGaaa                                                                                                              >  2:1622998/1‑139 (MQ=255)
                         cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaaggac                                                                                                                                >  2:583707/1‑115 (MQ=255)
                         cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaaggac                                                                                                                                <  1:583707/120‑6 (MQ=255)
                            tttttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGaaaaaa                                                                                                           >  1:665393/1‑139 (MQ=255)
                               ttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:322014/113‑5 (MQ=255)
                               ttAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  1:322014/1‑109 (MQ=255)
                                tAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTg                                                                                                       <  2:176552/139‑1 (MQ=255)
                                    aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcat                                                                                                   <  2:377315/139‑1 (MQ=255)
                                    aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  2:304302/108‑5 (MQ=255)
                                    aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  1:304302/1‑104 (MQ=255)
                                      ttttGCCAAACTCGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgc                                                                                                 >  1:152559/1‑139 (MQ=255)
                                       tttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  2:235249/139‑1 (MQ=255)
                                       tttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  1:639160/139‑1 (MQ=255)
                                        ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               <  2:244069/107‑1 (MQ=255)
                                        ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               >  1:244069/1‑107 (MQ=255)
                                        ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                               >  1:1547884/1‑139 (MQ=255)
                                                             cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                          <  1:1622998/139‑1 (MQ=255)
                                                              gAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  2:860345/84‑1 (MQ=255)
                                                              gAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  1:860345/1‑84 (MQ=255)
                                                                 aCTGCACAAATCTGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTg                                                                                                       >  1:430311/1‑106 (MQ=255)
                                                                 aCTGCACAAATCTGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTg                                                                                                       <  2:430311/106‑1 (MQ=255)
                                                                  cTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGa                                                                     <  2:152559/139‑1 (MQ=255)
                                                                   tGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCAATCCCCAAAAGGACGCTGTTAATGAAGGAGATAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGAc                                                                    >  1:1085781/1‑139 (MQ=255)
                                                                     cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                           >  1:641157/1‑114 (MQ=255)
                                                                     cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                           <  2:641157/114‑1 (MQ=255)
                                                                     cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                  >  2:834729/1‑139 (MQ=255)
                                                                       naaaNNNGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  1:834729/138‑1 (MQ=255)
                                                                       caAATCCGGTGTGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                                <  2:1085781/139‑1 (MQ=255)
                                                                          aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGATGGAGAAAAAATCTgg                                                                                                      >  2:905540/1‑98 (MQ=255)
                                                                          aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGATGGAGAAAAAATCTgg                                                                                                      <  1:905540/98‑1 (MQ=255)
                                                                          aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgc                                                                                                     <  1:576818/99‑1 (MQ=255)
                                                                          aTCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgc                                                                                                     >  2:576818/1‑99 (MQ=255)
                                                                             cGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTgtgt                                                          <  1:834937/139‑1 (MQ=255)
                                                                               gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGATAAAAAATCTGGCATGCATATCCCTCTTATTGc                                                                                 <  1:993244/114‑1 (MQ=255)
                                                                               gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGc                                                                                 >  2:993244/1‑114 (MQ=255)
                                                                                  cGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACg                                                     >  2:893590/1‑139 (MQ=255)
                                                                                        gTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                               <  1:893590/139‑1 (MQ=255)
                                                                                                  aacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               <  2:1044764/49‑1 (MQ=255)
                                                                                                  aacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGc                                                                                                                               >  1:1044764/1‑49 (MQ=255)
                                                                                                  aacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                     <  1:562530/139‑1 (MQ=255)
                                                                                                       gcgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                <  1:809806/43‑1 (MQ=38)
                                                                                                       gcgcCGAAGAGCAGGTAAATCATTACCGATCCCCAAAAGGACg                                                                                                                                >  2:809806/1‑43 (MQ=38)
                                                                                                        cgcCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 >  2:1340891/1‑36 (MQ=255)
                                                                                                        cgcCGAAGAGCAGGTAAATCATTACCGATCCCCaaaagga                                                                                                                                 <  1:1340891/40‑5 (MQ=255)
                                                                                                              agagCAGATAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                     >  1:435897/1‑127 (MQ=255)
                                                                                                              agagCAGATAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                     <  2:435897/127‑1 (MQ=255)
                                                                                                                     gTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCaa                                              >  2:1370846/1‑111 (MQ=255)
                                                                                                                     gTAAATCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCaa                                              <  1:1370846/111‑1 (MQ=255)
                                                                                                                       aaaTCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     >  1:374081/1‑134 (MQ=255)
                                                                                                                       aaaTCATTACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                     <  2:374081/134‑1 (MQ=255)
                                                                                                                             ttACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcat                                                                                                   >  2:427480/1‑50 (MQ=255)
                                                                                                                             ttACCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcat                                                                                                   <  1:427480/50‑1 (MQ=255)
                                                                                                                             ttACCGATCCCCAAAAGGACGCTGTTAACGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGAt          <  2:1000320/139‑1 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                               >  2:1478270/1‑52 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                               <  1:1478270/52‑1 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAc                                                      <  2:54356/93‑1 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAAc                                                      >  1:54356/1‑93 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACg                          <  2:474382/121‑1 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACg                          >  1:474382/1‑121 (MQ=255)
                                                                                                                               aCCGATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATaa        >  1:293156/1‑139 (MQ=255)
                                                                                                                                  gATCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACtt     >  2:378840/1‑139 (MQ=255)
                                                                                                                                  atccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                          <  1:1632488/66‑1 (MQ=255)
                                                                                                                                   tccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTAATTGATAACt      >  2:985431/3‑136 (MQ=255)
                                                                                                                                   aTCCCCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                          >  2:1632488/1‑69 (MQ=255)
                                                                                                                                    tCCCCACAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACt      <  1:985431/136‑1 (MQ=255)
                                                                                                                                    ccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                                <  2:793227/44‑1 (MQ=255)
                                                                                                                                    ccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                           >  1:1063584/2‑50 (MQ=255)
                                                                                                                                    ccccaAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTg  <  2:1217541/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCTG  >  NC_000913/7786‑8056

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