Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 661,369 T→C 100% N4D (AAT→GAT)  ybeF ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913661,3690TC97.0% 316.4 / ‑0.8 100N4D (AAT→GAT) ybeFLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (55/41);  total (57/43)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 105 total aligned reads displayed.

ATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCT  >  NC_000913/661272‑661464
                                                                                                 |                                                                                               
aTAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCAc                                                                                               <  1:2009477/100‑1 (MQ=255)
  tctCAGAAGGTTAAGATCAATATTTCGTAAGGTTTTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTa                                                                                             <  1:168835/98‑1 (MQ=255)
   gTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTAt                                                                                            >  1:1563638/1‑100 (MQ=255)
     cAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTaaga                                                                                          >  2:168833/1‑97 (MQ=255)
     cAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATcc                                                                                          >  2:908622/1‑100 (MQ=255)
          ggTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCa                                                                                     <  1:1847380/100‑1 (MQ=255)
          ggTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCa                                                                                     >  1:2028051/1‑100 (MQ=255)
          ggTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCa                                                                                     <  2:2336457/100‑1 (MQ=255)
           gTTAAGATCAATATTTCGTAAGGTTGTAAATAGTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGCTCAATTTGATTATCACTATCTACGCa                                                                                     >  1:2890549/1‑99 (MQ=255)
            ttAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACATGGTTCAATTTGATTATCTCTATCCACGCatt                                                                                   <  1:2047812/100‑3 (MQ=255)
              aaGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTa                                                                                 >  2:2044485/1‑100 (MQ=255)
               aGATCAATTTTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTAc                                                                                <  1:1891553/100‑1 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  2:773300/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  1:2836338/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  2:1480371/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  2:1999794/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  1:106711/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  2:2380639/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGCTCAATTTGATTATCACTATCCACGCACTTACt                                                                               >  2:392254/1‑100 (MQ=255)
                gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCGTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                               <  2:2901034/100‑1 (MQ=255)
                   cAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCa                                                                            >  1:1618796/1‑100 (MQ=255)
                     atatTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAt                                                                          <  1:1311945/100‑1 (MQ=255)
                     atatTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAt                                                                          <  1:1204192/100‑1 (MQ=255)
                       atTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAttt                                                                        <  2:106710/100‑1 (MQ=255)
                        tttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAtttt                                                                       >  2:469623/1‑100 (MQ=255)
                        tttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAtttt                                                                       >  2:181904/1‑100 (MQ=255)
                        tttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAtttt                                                                       >  1:2361867/1‑100 (MQ=255)
                         ttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTa                                                                      <  1:2169518/100‑1 (MQ=255)
                          tCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTAt                                                                     <  1:1652292/100‑1 (MQ=255)
                           cGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTAtt                                                                    <  2:1270164/100‑1 (MQ=255)
                            gTAAGGTTGTAAATATTTGTGGTTTATCTTCTGACGACTTTCTGCTTAAACAGGTTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTc                                                                   <  2:53565/100‑1 (MQ=255)
                            gTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTc                                                                   <  2:2868585/100‑1 (MQ=255)
                            gTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTc                                                                   <  2:2868755/100‑1 (MQ=255)
                             tAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCa                                                                  >  1:847442/1‑100 (MQ=255)
                              aaGGTTGTAAATATTTGTTGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCAt                                                                 <  2:1742782/100‑1 (MQ=255)
                               aGGTTGTAAATATTTGTGGTTTACCTTCTGTCGTCTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATg                                                                <  1:1605579/100‑1 (MQ=255)
                                 gTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGa                                                              <  1:1098932/100‑1 (MQ=255)
                                    gTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGaaaa                                                           >  2:2227484/1‑100 (MQ=255)
                                    gTAAATATTTGTGGTTTACCTTCTGAAGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGaaaa                                                           <  2:552147/100‑1 (MQ=255)
                                     tAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGaaaaa                                                          >  2:2322867/1‑100 (MQ=255)
                                         tatTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaataat                                                      <  2:1388887/100‑1 (MQ=255)
                                           tttGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACGTACTCCAATTTTATTCATGGAAAAATAatat                                                    <  1:1156333/100‑1 (MQ=255)
                                             tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                  >  1:1060568/1‑100 (MQ=255)
                                             tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                  >  1:694954/1‑100 (MQ=255)
                                             tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                  <  1:2942124/100‑1 (MQ=255)
                                             tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                  <  1:2940364/100‑1 (MQ=255)
                                              gtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCACTTTTATTCATGGAAAAATAATAttt                                                  >  1:1499786/1‑99 (MQ=255)
                                               tgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaa                                                <  1:2184852/100‑1 (MQ=255)
                                                ggTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaa                                               >  2:198038/1‑100 (MQ=255)
                                                ggTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaa                                               >  2:2156569/1‑100 (MQ=255)
                                                ggTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaa                                               <  2:2973750/100‑1 (MQ=255)
                                                ggTTTACCTTCTGACGACATTCTGCTTAAACAGGGGTTAATTTGATTATCACTATCCACGCACTTACTCCAATTTTTTTCATGGTAAAATAATATTTaaa                                               >  2:952908/1‑100 (MQ=255)
                                                   ttACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaaaaa                                            <  2:1577217/100‑1 (MQ=255)
                                                   ntACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaaaaa                                            >  2:1506402/2‑100 (MQ=255)
                                                    tACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                           >  2:1511184/1‑100 (MQ=255)
                                                     aCCTTCTGACGACTTTCTGCTTAAACATTGTTCAATTTGATTATCACTATCTACTCACTTACTCCAATTTTATTCATTTAAAAATAATATTTAAAAAAtt                                          <  2:951176/100‑1 (MQ=255)
                                                      ccTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTa                                         >  2:1422579/1‑100 (MQ=255)
                                                      ccTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTa                                         <  1:2156592/100‑1 (MQ=255)
                                                       cTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAc                                        <  2:1626624/100‑1 (MQ=255)
                                                       cTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAc                                        <  1:2380666/100‑1 (MQ=255)
                                                       cTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAc                                        <  2:1170421/100‑1 (MQ=255)
                                                        ttCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACa                                       <  2:1050827/100‑1 (MQ=255)
                                                         tCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaa                                      >  2:197982/1‑100 (MQ=255)
                                                          cTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaat                                     >  2:237225/1‑100 (MQ=255)
                                                           tGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaata                                    >  2:2196554/1‑100 (MQ=255)
                                                            gacgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaataa                                   >  1:1062964/1‑100 (MQ=255)
                                                            gacgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaataa                                   >  1:120572/1‑100 (MQ=255)
                                                            gacgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaataa                                   >  1:2436582/1‑100 (MQ=255)
                                                             acgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaataat                                  >  1:2192239/1‑100 (MQ=255)
                                                             acgacTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACaataat                                  >  1:1375849/1‑100 (MQ=255)
                                                                acTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCtt                               >  2:1127635/1‑100 (MQ=255)
                                                                acTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCtt                               <  2:326791/100‑1 (MQ=255)
                                                                 cTTTCTGCTTTAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGAAACAATAATATTTAAAAAATTACAATAATCTTa                              <  2:1997484/100‑1 (MQ=255)
                                                                   ttCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATg                            >  2:949608/1‑100 (MQ=255)
                                                                     cTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTc                          >  2:851288/1‑100 (MQ=255)
                                                                     cTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTc                          <  2:1467160/100‑1 (MQ=255)
                                                                      tGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCt                         >  2:1599763/1‑100 (MQ=255)
                                                                      tGCTTAAACAGGGTTCAAATTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTAAAAAAATTACAATATTCTTATGTCt                         >  2:1457797/1‑100 (MQ=255)
                                                                        cTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTaa                       >  1:1969694/1‑100 (MQ=255)
                                                                        cTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                        >  2:1459015/1‑99 (MQ=255)
                                                                            aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                   >  2:981282/1‑100 (MQ=255)
                                                                            aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                   >  2:788090/1‑100 (MQ=255)
                                                                            aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                   >  2:1556740/1‑100 (MQ=255)
                                                                             aCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgga                  <  1:281762/100‑1 (MQ=255)
                                                                             aCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTaga                  <  2:1418470/100‑1 (MQ=255)
                                                                             aCAGAGTTCAATTTGATTATCACTATCCACGCACTTATTCCAATTAAAAAAAAAAAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgga                  <  2:2084739/100‑1 (MQ=255)
                                                                              cAGTGTTCAATTTGATTATCACTATCCACGCACTTAATCCAATTTTATTTATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgga                  <  2:2745285/99‑1 (MQ=255)
                                                                              cAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaa                 >  2:2834031/1‑100 (MQ=255)
                                                                                  gTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacgga             >  1:1846235/1‑100 (MQ=255)
                                                                                   tttaaTTTGATTATCACTATCCACGCATTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacggaa            <  1:2750932/97‑1 (MQ=255)
                                                                                     cAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGgaacggaacg          >  2:693710/1‑100 (MQ=255)
                                                                                      aaTTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGgaacggaacg          >  1:2164715/1‑99 (MQ=255)
                                                                                         ttcattatcACTATTTACGCACTTACTCCAATTGTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtt       <  2:1802575/96‑1 (MQ=255)
                                                                                          tGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtttt     >  1:2729310/1‑100 (MQ=255)
                                                                                          tGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtttt     <  2:2406000/100‑1 (MQ=255)
                                                                                           gATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtttt     >  2:1097932/1‑99 (MQ=255)
                                                                                            attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTg  <  2:1737446/100‑2 (MQ=255)
                                                                                            attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc   <  1:1423508/100‑1 (MQ=255)
                                                                                            attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc   <  1:1209969/100‑1 (MQ=255)
                                                                                             ttatcaCTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCt  >  1:2308209/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ATAGTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGCT  >  NC_000913/661272‑661464

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