Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC22.3% 292.1 / 77.3 130intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (51/50);  new base C (27/2);  total (78/52)
Fisher's exact test for biased strand distribution p-value = 2.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 130 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705887‑3706080
                                                                                                  |                                                                                                       
     cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                   <  1:215620/100‑1 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                  >  1:648711/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                  <  2:1132136/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  2:1885251/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  1:1557532/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               >  1:416082/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                        <  1:2256075/100‑1 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                        >  1:203695/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                        <  1:1632161/100‑1 (MQ=255)
                 cgatctgagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  2:1800613/94‑1 (MQ=255)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  2:573304/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAAGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    >  1:782115/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCgagatcg                                                                                >  1:1800623/1‑94 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  1:67463/1‑100 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              <  2:1497895/100‑1 (MQ=255)
                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  2:765012/1‑100 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                             <  1:768207/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                             >  1:2057318/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGCTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            >  1:1523869/1‑100 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                            <  2:2157374/100‑1 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                         <  1:1054144/100‑1 (MQ=255)
                               gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                        >  2:1042361/1‑100 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                      >  1:835494/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     >  1:759855/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  1:398468/100‑1 (MQ=255)
                                  ccACCGGGTACACAGATTGTTGACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                     <  2:1475847/100‑1 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                            >  1:269434/1‑100 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                            >  1:1125812/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                            >  2:1330941/1‑98 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑Att‑‑‑tt                                                      <  2:673778/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑Att‑‑‑tt                                                      >  2:1740457/1‑100 (MQ=255)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTc                                                     <  1:1034851/100‑1 (MQ=255)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTc                                                     >  1:434616/1‑100 (MQ=255)
                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTc                                                     >  1:1485433/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑g                                                  <  2:1234742/100‑1 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑Gt                                                 <  1:1740466/100‑1 (MQ=255)
                                                gATTTTTGACAAAGTGCGCTTTGTTCATGCCGGATGCGACGTGAACGTCTTATCTGGCCTAC‑AAGACCATGCAAATTCAATATATTGCA‑Gatcggaagag                                                      >  2:1208088/1‑91 (MQ=21)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTa                                                <  1:2201934/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTa                                                >  1:745954/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTa                                                >  2:707993/1‑100 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAg                                               <  1:1677423/100‑1 (MQ=255)
                                                 aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAg                                               <  2:2057303/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGcc                                            >  2:490954/1‑100 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCt                                           <  2:1903414/100‑1 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTg                                          >  2:1350749/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGa                                         >  1:388103/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGa                                         >  1:221989/1‑100 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGAt                                        <  1:1990273/100‑1 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGAt                                        <  1:966821/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATa                                       <  1:374849/100‑1 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATa                                       <  2:2043059/99‑1 (MQ=255)
                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                  <  1:1042368/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                 >  1:131870/1‑100 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                <  1:1330950/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGAtgcgctccgc                               >  1:2176556/1‑91 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                               <  2:147348/100‑1 (MQ=255)
                                                                  gattttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAgcg                              >  1:2155339/11‑100 (MQ=18)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                              <  2:848863/100‑1 (MQ=255)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                            <  2:1523860/100‑1 (MQ=255)
                                                                    atttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCAt                           >  1:686133/9‑100 (MQ=21)
                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                           >  1:864845/1‑100 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑ATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATc                          >  1:206809/7‑100 (MQ=18)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATc                          >  1:274615/7‑100 (MQ=18)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                          <  2:1895739/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                         <  2:835487/100‑1 (MQ=255)
                                                                        ttATGCCCGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                        >  2:1188753/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                        >  1:2173221/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                        >  1:134630/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  1:1990739/4‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                       >  2:436859/4‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg                       >  2:1561143/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑ACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  2:1448688/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                      >  2:2077417/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                      >  1:648648/3‑100 (MQ=25)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑ATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                      >  1:787673/3‑100 (MQ=17)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                     >  1:384546/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGATT‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                     >  1:1138604/2‑100 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑ATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                    >  1:1216802/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                    >  1:417113/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                    >  2:4891/1‑99 (MQ=21)
                                                                             ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                   >  1:608759/9‑98 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                  <  2:686127/92‑1 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑ACG‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                 >  2:1019628/7‑96 (MQ=25)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg               <  1:1482651/100‑1 (MQ=255)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc              <  2:221986/100‑1 (MQ=255)
                                                                                   gcggcgTTAGCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc             <  1:1728581/100‑1 (MQ=255)
                                                                                   gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc             <  2:648707/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct            >  2:856877/1‑100 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct            >  2:225385/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt           >  1:976823/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt           >  1:2088813/1‑100 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc        <  1:1072804/100‑1 (MQ=255)
                                                                                        gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc        >  1:818232/1‑100 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc       <  1:412365/100‑1 (MQ=255)
                                                                                         ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATA‑‑‑TTC‑‑GTAAGCATGATAAGCTGAGCGCATCAGGCATTtttccctttgcc       >  2:906042/1‑100 (MQ=18)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc    <  1:1658935/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa   <  1:490957/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa   >  1:1769699/1‑100 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  <  1:492892/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg  >  1:1814811/1‑100 (MQ=255)
                                                                                                  |                                                                                                       
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAG  >  NC_000913/3705887‑3706080

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