Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,412,000 C→T *412* (TAG→TAA intR ← Rac prophage, integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,412,0000CT100.0% 392.0 / NA 121*412* (TAG→TAAintRRac prophage, integrase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (59/62);  total (59/62)

Only 100 of 121 total aligned reads displayed.

TCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGG  >  NC_000913/1411866‑1412123
                                                                                                                                      |                                                                                                                           
tCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACg                                                                                                                         >  2:102652/1‑139 (MQ=255)
      aCGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAAcc                                                                                                                   <  1:1545136/139‑1 (MQ=255)
             gacgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                            <  2:1794165/139‑1 (MQ=255)
               cgCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCa                                                                                                          <  1:1687050/139‑1 (MQ=255)
                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                  <  1:917601/128‑1 (MQ=255)
                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                  >  2:917601/1‑128 (MQ=255)
                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                              >  2:1785655/1‑132 (MQ=255)
                  aGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTGTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGAAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                              <  1:1785655/132‑1 (MQ=255)
                    aCGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                  <  2:1666657/126‑1 (MQ=255)
                    aCGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                  >  1:1666657/1‑126 (MQ=255)
                         tttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGTATTCAAc                                                                                                        <  1:1177461/131‑1 (MQ=255)
                         tttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        >  2:1177461/1‑131 (MQ=255)
                         tttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATTGCTAt                                                                                                <  1:799654/139‑1 (MQ=255)
                          ttGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATc                                                                                               <  2:1063730/139‑1 (MQ=255)
                               aattGTTCACGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCTTTCAACATGGCTATCTGTtc                                                                                          <  1:1057751/137‑1 (MQ=255)
                                  tGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcgtc                                                                                       >  1:221095/1‑139 (MQ=255)
                                  tGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcgtc                                                                                       >  1:1126777/1‑139 (MQ=255)
                                    ttCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTcgtcgt                                                                                     >  2:1461454/1‑139 (MQ=255)
                                     tCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGtt                                                                                    <  1:102652/139‑1 (MQ=255)
                                       aGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCa                                                                                  <  1:808347/139‑1 (MQ=255)
                                          ttGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGt                                                                               <  2:221095/139‑1 (MQ=255)
                                           tGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtc                                                                              >  1:273871/1‑139 (MQ=255)
                                              aTTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                           <  2:1126777/139‑1 (MQ=255)
                                                tGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCaa                                                                         <  1:1506532/139‑1 (MQ=255)
                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                 <  1:1021234/128‑1 (MQ=255)
                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             <  2:1345885/132‑1 (MQ=255)
                                                   tctttcttTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             >  1:1345885/1‑132 (MQ=255)
                                                   tctttcttTGTAATTTGTTCATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATcc                                                                      >  2:1196438/1‑139 (MQ=255)
                                                    ctttctttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAGATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCca                                                                     >  1:413121/1‑139 (MQ=255)
                                                          ttGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTtcg                                                                                         >  1:1153654/1‑113 (MQ=255)
                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                             >  2:714986/1‑139 (MQ=255)
                                                            gTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                             >  2:553007/1‑139 (MQ=255)
                                                              aaTTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           <  1:553007/139‑1 (MQ=255)
                                                                  tgttgATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACg                                                                                                                         <  2:1228794/73‑1 (MQ=255)
                                                                  tgttgATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACg                                                                                                                         >  1:1228794/1‑73 (MQ=255)
                                                                  tgttgATTTTCTTGCATTATTTCAGTTATCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACcg                                                                                                                  >  1:417182/1‑80 (MQ=255)
                                                                   gttgATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAt                                                                                                             <  2:567532/84‑1 (MQ=255)
                                                                   gttgATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAt                                                                                                             >  1:567532/1‑84 (MQ=255)
                                                                         tttCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           >  2:1152220/1‑128 (MQ=255)
                                                                         tttCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           <  1:1152220/128‑1 (MQ=255)
                                                                            cTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGTGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTGCATGtcat                                                                            <  1:1242394/108‑1 (MQ=255)
                                                                            cTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcat                                                                            >  2:1242394/1‑108 (MQ=255)
                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             >  1:1618224/1‑98 (MQ=255)
                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             <  2:1618224/98‑1 (MQ=255)
                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcaca                                                                   >  1:1672177/1‑108 (MQ=255)
                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCcaca                                                                   <  2:1672177/108‑1 (MQ=255)
                                                                                     tttCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGt                                    >  1:1492427/1‑139 (MQ=255)
                                                                                         aGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCa                                <  1:714986/139‑1 (MQ=255)
                                                                                           ttCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        <  2:679126/65‑1 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCATATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        >  1:1054456/1‑63 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                            <  2:1666019/59‑1 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCAtt                                                                                                            >  1:1666019/1‑59 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        <  2:1054456/63‑1 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        >  1:1400140/1‑63 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        <  2:1400140/63‑1 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                             >  1:1180660/1‑106 (MQ=255)
                                                                                             ctctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGt                                                             <  2:1180660/106‑1 (MQ=255)
                                                                                               ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATcc                                                                      >  1:1382230/1‑95 (MQ=255)
                                                                                               ctGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           >  2:1738289/1‑106 (MQ=255)
                                                                                                tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             >  1:652458/1‑87 (MQ=255)
                                                                                                tGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             <  2:652458/87‑1 (MQ=255)
                                                                                                     cTAAATGGGGCAGATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  2:413121/139‑1 (MQ=255)
                                                                                                     cTAAATGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  1:992258/139‑1 (MQ=255)
                                                                                                        aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGc                                                                                                   >  2:1362938/1‑57 (MQ=255)
                                                                                                        aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           >  1:1370839/1‑97 (MQ=255)
                                                                                                        aaTGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           <  2:1370839/97‑1 (MQ=255)
                                                                                                          tGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtca                                                                             >  1:1232471/1‑77 (MQ=255)
                                                                                                          tGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  1:902539/134‑1 (MQ=255)
                                                                                                          tGGGGCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    >  2:902539/1‑134 (MQ=255)
                                                                                                            gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                              <  1:577467/122‑1 (MQ=255)
                                                                                                            gggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                              >  2:577467/1‑122 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                              >  2:1579110/1‑41 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                              <  1:1579110/41‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                 >  1:775853/1‑70 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCAt                                                                                 <  2:775853/70‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                           <  1:34572/76‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                           >  1:1005836/1‑76 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                           <  2:1005836/76‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGtcatc                                                                           >  2:34572/1‑76 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                             >  2:274026/1‑106 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCAt                                             <  1:274026/106‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  2:1453502/131‑1 (MQ=255)
                                                                                                             ggCAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCtat            >  1:423033/1‑139 (MQ=255)
                                                                                                               cAAATTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATaa          >  2:547835/1‑139 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCa                                                                                                              >  2:462081/1‑36 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        >  2:413742/1‑42 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAAc                                                                                                        <  1:413742/42‑1 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                  <  2:317959/80‑1 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACAt                                                                  >  1:317959/1‑80 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  1:1278330/126‑1 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    >  2:1278330/1‑126 (MQ=255)
                                                                                                                  aTTGGGGGCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATg       >  1:1385318/1‑139 (MQ=255)
                                                                                                                       gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           <  1:76672/82‑1 (MQ=255)
                                                                                                                       gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTaa                                                           >  2:76672/1‑82 (MQ=255)
                                                                                                                       gggCAAACTTTGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCggg  <  2:1385318/139‑1 (MQ=255)
                                                                                                                                ttGCAATTACGATAACCGCGCATTCAACATGGCTAt                                                                                                >  2:888238/1‑36 (MQ=255)
                                                                                                                                 tGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                              <  1:1413000/101‑1 (MQ=255)
                                                                                                                                 tGCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATg                              >  2:1413000/1‑101 (MQ=255)
                                                                                                                                  gCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    >  2:96675/1‑110 (MQ=255)
                                                                                                                                  gCAATTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTtgg                    <  1:96675/110‑1 (MQ=255)
                                                                                                                                      |                                                                                                                           
TCGCCCACGTTACGACGCAGACGTTTTTGTAATTTGTTCAGGTTGTATTGTTCTTTCTTTGTAATTTGTTGATTTTCTTGCATTATTTCAGTTCTCTGGTACTAAATGGGGCAAATTGGGGGCAAACTTTGCAACTACGATAACCGCGCATTCAACATGGCTATCTGTTCGTCGTTCATGTCATCAATCCACATACCGTAAATTTCATACACCATCTGCGCAGTTTCATGCCCCATTTGGCTGGCTATAAATGCCGGG  >  NC_000913/1411866‑1412123

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: