Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9700GA35.1% 302.2 / 211.1 211intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (74/63);  new base A (63/11);  total (137/74)
Fisher's exact test for biased strand distribution p-value = 4.44e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

Only 100 of 225 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                          
  gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  1:2654886/1‑100 (MQ=255)
   ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                         <  1:826429/100‑1 (MQ=255)
    gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  2:932983/100‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:1580793/100‑1 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   >  1:203148/1‑100 (MQ=255)
          aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                  >  1:224525/1‑100 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                  >  1:1824749/1‑99 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  1:1253389/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  2:10896/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  2:1651575/100‑1 (MQ=255)
            tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  2:3115030/100‑1 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                               <  2:1659285/100‑1 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              >  1:1087709/1‑100 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  2:2466725/100‑1 (MQ=255)
                   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         >  1:15593/1‑100 (MQ=255)
                     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                       >  1:101277/1‑100 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      >  2:2557413/1‑100 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  1:377218/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:1139664/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:1002249/100‑1 (MQ=255)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                     <  1:1807319/100‑2 (MQ=37)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                    <  2:520025/100‑3 (MQ=37)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                  >  2:588585/1‑100 (MQ=255)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                  <  1:1845536/100‑2 (MQ=25)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                <  2:210816/100‑1 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                 <  1:1461271/100‑3 (MQ=25)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                               >  1:1930134/1‑100 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                <  2:2290683/100‑7 (MQ=25)
                               tctaCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                            <  1:1624342/97‑1 (MQ=255)
                               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                            >  1:729796/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAtaga                                                                         >  2:1624338/1‑97 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                         >  2:1577149/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        <  2:1521408/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                  >  1:1948675/1‑100 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                  <  1:480460/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  2:1427971/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  2:2838098/100‑1 (MQ=255)
                                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  1:216976/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                >  2:1798117/1‑100 (MQ=255)
                                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                               >  2:2759536/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              <  2:529687/100‑1 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              >  1:1917666/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              <  2:1059568/100‑1 (MQ=255)
                                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                          <  2:1494721/100‑1 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                         >  2:181037/1‑100 (MQ=255)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                         <  2:118940/100‑1 (MQ=255)
                                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                        >  2:1264582/1‑100 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                       <  1:1762730/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                       >  1:1597698/1‑100 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                      <  1:2274060/100‑1 (MQ=255)
                                                      cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                    <  1:2392299/100‑1 (MQ=255)
                                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t‑‑‑t                                          >  1:3352256/1‑100 (MQ=255)
                                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t‑‑‑t                                          <  1:2149524/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTa                                     >  1:1043269/1‑100 (MQ=255)
                                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAg                                    >  2:369423/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGcc                                 <  2:718562/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGcc                                 <  2:2773431/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGcc                                 <  2:1650893/100‑1 (MQ=255)
                                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCt                                >  1:2843728/1‑100 (MQ=255)
                                                                     taaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCt                                >  1:1291537/4‑100 (MQ=18)
                                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGa                              >  1:58313/1‑100 (MQ=255)
                                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGAt                             >  2:3317488/1‑100 (MQ=255)
                                                                                   attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAgcgc                  >  1:1539805/10‑100 (MQ=18)
                                                                                     tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATc               >  1:669997/8‑100 (MQ=18)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCAt                >  1:1736916/1‑100 (MQ=255)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc               <  2:349710/100‑1 (MQ=255)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:1326802/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCa              >  2:2406714/6‑100 (MQ=18)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg             >  1:1336780/1‑100 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:738645/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2763613/5‑100 (MQ=17)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:2228552/5‑100 (MQ=18)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:2292466/5‑100 (MQ=18)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:1840102/4‑100 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAgg            >  1:3230938/4‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:3119162/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:602221/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg            >  2:2588132/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg            >  1:2498707/3‑100 (MQ=21)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:50031/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:1592886/3‑100 (MQ=18)
                                                                                           tGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGGGCAAAGTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  2:2675484/1‑100 (MQ=255)
                                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa          >  1:1370537/1‑100 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑AT‑TC‑T‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGct          >  1:986839/2‑99 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAATCCTTGCGAATTCAATATATTGTATGG‑AA‑A‑‑‑TTC‑‑AGGCCTGATAAGCGTAGCGCATCAGGctt         >  2:1444908/2‑98 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:3014478/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:2143376/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑‑‑CACGTAGGCCTGATAAGCGAAGCGCATCAGGCa          >  2:2893934/2‑100 (MQ=18)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:2568400/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:3344299/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgagat         >  2:1892466/1‑96 (MQ=21)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GG‑CT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1417800/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa         >  1:846655/1‑99 (MQ=25)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑C‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1213618/1‑99 (MQ=21)
                                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCagt        >  2:310193/1‑98 (MQ=17)
                                                                                                 aTGCGGCGTTAACGCCTTATCTGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:2056176/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑CT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  >  2:783277/1‑100 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  1:3067158/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   >  1:1156189/1‑99 (MQ=255)
                                                                                                   gcggcgTTAACAGGGTATCCGGCTTAT‑AATATCGTGCAAATTCAAGACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc   <  2:3032259/99‑1 (MQ=18)
                                                                                                   |                                                                                                          
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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