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.5% 206.9 / ‑3.9 65coding (38/1431 nt)yaaJputative transporter
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (29/35);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7788‑8054
                                                                                                                                        |                                                                                                                                    
aaCCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaa                                                                                                                                  >  1:172243/1‑138 (MQ=255)
    gCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaaggac                                                                                                                              <  1:1729355/139‑6 (MQ=255)
     ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGTGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACg                                                                                                                              <  1:1523972/139‑1 (MQ=255)
     ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACg                                                                                                                              <  2:1419840/139‑1 (MQ=255)
      cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               >  2:1493134/1‑132 (MQ=255)
      cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               <  1:1493134/136‑5 (MQ=255)
       tGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGTGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCt                                                                                                                            <  1:870698/139‑1 (MQ=255)
               gATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGa                                                                                                                    >  1:67902/1‑139 (MQ=255)
               gATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGa                                                                                                                    <  1:1623719/139‑1 (MQ=255)
                                  aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               >  2:1733745/1‑104 (MQ=255)
                                  aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               <  1:1733745/108‑5 (MQ=255)
                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATg                                                                                                                     >  1:1316083/1‑117 (MQ=255)
                                    ttttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATg                                                                                                                     <  2:1316083/117‑1 (MQ=255)
                                         ccAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          <  2:67902/139‑1 (MQ=255)
                                          cAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         <  2:1397314/139‑1 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCatat                                                                                           >  2:42081/1‑130 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCatat                                                                                           <  1:42081/130‑1 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTAt                                                                                  >  1:313081/1‑139 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTAt                                                                                  >  1:746567/1‑139 (MQ=255)
                                                 gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTAt                                                                                  >  1:1543886/1‑139 (MQ=255)
                                                        tAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACg                                                                                                                              <  1:865637/88‑1 (MQ=255)
                                                        tAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACg                                                                                                                              >  2:865637/1‑88 (MQ=255)
                                                           cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                        <  2:1543886/139‑1 (MQ=255)
                                                           cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                        >  1:408754/1‑139 (MQ=255)
                                                           cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                        <  2:746567/139‑1 (MQ=255)
                                                               aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGAGGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATg                                                                    >  1:1186650/1‑139 (MQ=255)
                                                               aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATg                                                                    >  2:267461/1‑139 (MQ=255)
                                                               aCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATg                                                                    >  2:1204483/1‑139 (MQ=255)
                                                                 tGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGAc                                                                  >  1:125813/1‑139 (MQ=255)
                                                                  gCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACt                                                                 <  1:267461/139‑1 (MQ=255)
                                                                   cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGa                                                                                                                    <  2:328795/87‑1 (MQ=255)
                                                                   cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGa                                                                                                                    >  1:328795/1‑87 (MQ=255)
                                                                   cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACtt                                                                >  1:282131/1‑139 (MQ=255)
                                                                     caAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGAGGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:1186650/139‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTTAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          <  1:1028278/103‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGaa                                                                                                             >  1:1703840/1‑84 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGaa                                                                                                             <  2:1703840/84‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                          >  2:1028278/1‑103 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  1:1561189/131‑1 (MQ=255)
                                                                             gTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  2:1561189/1‑131 (MQ=255)
                                                                                  aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACg                                                                                                                              <  1:1067284/62‑1 (MQ=255)
                                                                                  aaaaGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCACTACCGATCCTCA‑AAAGGACg                                                                                                                              >  2:1067284/1‑62 (MQ=255)
                                                                                         aaCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                      >  1:1267664/1‑127 (MQ=255)
                                                                                         aaCCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTaa                                                      <  2:1267664/127‑1 (MQ=255)
                                                                                              acaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              >  1:401610/1‑114 (MQ=255)
                                                                                              acaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTc                                                              <  2:401610/114‑1 (MQ=255)
                                                                                              acaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               <  2:1000247/48‑5 (MQ=255)
                                                                                              acaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               >  1:1000247/1‑44 (MQ=255)
                                                                                               caacCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                    <  1:1071536/139‑1 (MQ=255)
                                                                                                 acCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  <  2:313081/139‑1 (MQ=255)
                                                                                                    tGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTGTTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             <  2:374868/125‑1 (MQ=255)
                                                                                                    tGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGATAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             >  1:374868/1‑125 (MQ=255)
                                                                                                     gcgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  <  2:931483/135‑1 (MQ=255)
                                                                                                     gcgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                  >  1:931483/1‑135 (MQ=255)
                                                                                                     gcgcCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               >  1:814526/1‑37 (MQ=255)
                                                                                                     gcgcCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑aaaagga                                                                                                                               <  2:814526/41‑5 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         <  1:1706511/79‑1 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                         >  2:1706511/1‑79 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                           <  2:1471834/125‑1 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAAt                                           >  1:1471834/1‑125 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                    >  2:1322271/1‑132 (MQ=255)
                                                                                                      cgcCGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaa                                    <  1:1322271/132‑1 (MQ=255)
                                                                                                        ccGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             >  2:928039/1‑121 (MQ=255)
                                                                                                        ccGAAGAGCAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             <  1:928039/121‑1 (MQ=255)
                                                                                                                cAGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACg                   <  1:749728/139‑1 (MQ=255)
                                                                                                                 aGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                      <  2:258093/119‑1 (MQ=255)
                                                                                                                 aGGTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                      >  1:258093/1‑119 (MQ=255)
                                                                                                                   gTAAATCATTACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAgg                <  2:125813/139‑1 (MQ=255)
                                                                                                                            tACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGa                                                                   <  2:413082/79‑1 (MQ=255)
                                                                                                                            tACCGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGa                                                                   >  1:413082/1‑79 (MQ=255)
                                                                                                                              ccGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                  >  2:267196/1‑46 (MQ=255)
                                                                                                                              ccGATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                  <  1:267196/46‑1 (MQ=255)
                                                                                                                                gATCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACtt   <  1:1204483/139‑1 (MQ=255)
                                                                                                                                atccc‑ca‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCGGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc  >  2:1134529/4‑139 (MQ=255)
                                                                                                                                  tCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             >  2:515815/1‑95 (MQ=255)
                                                                                                                                  tCC‑CCAAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCa                                             <  1:515815/95‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
AACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCC‑CC‑AAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTC  >  NC_000913/7788‑8054

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