Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9700GA36.0% 257.6 / 180.9 178intergenic (‑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 (42/72);  new base A (63/1);  total (105/73)
Fisher's exact test for biased strand distribution p-value = 1.73e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 186 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705873‑3706069
                                                                                                 |                                                                                                         
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        >  1:2038733/1‑100 (MQ=255)
 ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                        >  1:2212408/1‑100 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCAACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTan                                                                                                  <  2:565421/100‑2 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 >  1:1981007/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                 <  1:850059/100‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  1:759376/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              >  1:2227200/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              >  2:1822529/1‑100 (MQ=255)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         >  2:1347484/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:2447734/1‑100 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  2:2179961/100‑1 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     >  1:275696/1‑100 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               >  1:1464307/1‑100 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               <  2:2227178/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                              <  1:1403005/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                              >  2:728874/1‑100 (MQ=255)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                           >  2:2044380/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          >  1:2845674/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          <  2:226413/100‑1 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                        <  1:2173809/100‑1 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                        <  2:1717997/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACACCTTATCCGGCCTAC‑AAAATCGTg                                                                       >  2:1391728/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTAACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  2:431988/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  1:43043/100‑1 (MQ=255)
                                    agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                    <  2:2204977/100‑1 (MQ=255)
                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                  <  1:2939124/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                 >  2:898896/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                 <  2:3004931/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                <  1:2240940/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                               <  1:1652236/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                               <  1:832330/100‑1 (MQ=255)
                                           caGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                             <  1:1286753/100‑1 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                             >  1:1842321/1‑100 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                            <  1:979603/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                           <  2:2982537/100‑1 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGCTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                          >  1:433216/1‑99 (MQ=255)
                                                 atccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                       <  2:2403369/98‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                   <  2:70775/100‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                   <  1:1353440/100‑1 (MQ=255)
                                                    cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                   <  1:2558830/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                  >  1:658298/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                       <  2:2369302/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                       <  2:2478712/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                       <  2:2952264/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑Gt                                      <  2:947949/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑Gt                                      >  2:2335949/1‑100 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTa                                     <  1:1924279/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTGATGCCGGATGCGGCGTTAGCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCt                                <  2:2772865/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGGTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                              <  2:220848/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                              >  2:632887/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                              <  2:2683300/100‑1 (MQ=255)
                                                                     aaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                              >  1:928465/2‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGAt                             >  2:977422/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATa                            <  2:2844081/100‑1 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                       <  2:3061109/100‑1 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGAtgcgctcc                      >  1:460638/1‑93 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                      >  2:189679/1‑100 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                      <  2:2614562/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                    <  2:390685/100‑1 (MQ=255)
                                                                                gattttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   >  1:3066054/11‑100 (MQ=17)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  1:1788351/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  1:1879093/100‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                   <  2:211501/100‑1 (MQ=255)
                                                                                  ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCa                 >  2:1769391/9‑100 (MQ=21)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                >  1:1605329/1‑100 (MQ=255)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                <  2:2758137/100‑1 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATc               >  1:2104585/7‑100 (MQ=17)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc               >  1:1200845/7‑100 (MQ=18)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATc               >  1:1046052/7‑100 (MQ=17)
                                                                                     tgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa              >  1:2281188/6‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  2:1329482/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAg             >  2:2105427/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAg             >  1:502569/5‑100 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAACCGCGCAAATTCAATATATTGCG‑GA‑GA‑AAAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  2:1919829/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:1836426/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:2416884/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1027827/4‑100 (MQ=21)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg            <  2:2636986/100‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2436226/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:78956/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:1760422/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GG‑AT‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc           >  2:2741117/3‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:1984487/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:474452/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1835910/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:777102/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga         >  1:2884951/2‑98 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa          >  1:1372238/2‑100 (MQ=25)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:3055622/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:3055516/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:595827/2‑100 (MQ=21)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1452139/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:2873465/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:680887/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1181244/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa         >  1:96037/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  2:84063/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:1001982/1‑99 (MQ=21)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag         >  1:1479721/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa         >  1:2612654/1‑99 (MQ=17)
                                                                                                 |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705873‑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: ‑