Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,971	0	C	A	23.9%	232.3 / 65.5	88	intergenic (‑181/+127)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (38/29);  new base A (20/1);  total (58/30)
Fisher's exact test for biased strand distribution p-value = 1.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070                                                                                                  |                                                                                                        ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                       <  2:1521831/100‑1 (MQ=255) ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                       >  2:1115126/1‑100 (MQ=255)  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                      >  2:673744/1‑100 (MQ=255)    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                    >  2:836256/1‑100 (MQ=255)    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                    <  1:1355523/100‑1 (MQ=255)      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  >  2:528286/1‑100 (MQ=255)      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  >  1:469234/1‑100 (MQ=255)        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                <  2:1542076/100‑1 (MQ=255)         ggaTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                               <  1:973226/98‑1 (MQ=255)            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                            <  1:520856/100‑1 (MQ=255)              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACTGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                          >  1:986906/1‑100 (MQ=255)              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                          >  2:538285/1‑100 (MQ=255)                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                        >  2:339617/1‑100 (MQ=255)                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTACGCTTTGTTTATGCCGGATGCGGCGTTCACGCATTAtaagg                                                                                     >  2:931684/1‑96 (MQ=37)                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  1:1505984/100‑1 (MQ=255)                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    >  2:705120/1‑100 (MQ=255)                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:1064011/100‑1 (MQ=255)                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:888890/100‑1 (MQ=255)                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                 >  2:751751/1‑100 (MQ=255)                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  1:18705/1‑100 (MQ=255)                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  2:737076/1‑100 (MQ=255)                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  1:854914/1‑100 (MQ=255)                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATc                                                                         >  1:350972/1‑100 (MQ=255)                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                        >  1:634000/1‑100 (MQ=255)                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                        <  2:918989/100‑1 (MQ=255)                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                        <  2:302358/100‑1 (MQ=255)                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                       <  2:532615/99‑1 (MQ=255)                                   cAGAGGCACATCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                    <  2:1495518/100‑1 (MQ=255)                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                    >  1:100443/1‑100 (MQ=255)                                    agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                   <  2:907419/100‑1 (MQ=255)                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                 >  1:1144143/1‑100 (MQ=255)                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAG‑AAAATCGTGCAAAtt                                                                >  2:1575530/1‑100 (MQ=255)                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                >  1:170565/1‑100 (MQ=255)                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  2:1175553/100‑1 (MQ=255)                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                       >  2:1510720/1‑100 (MQ=255)                                                ccACCGGCTATGCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTGTCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                       >  1:842816/1‑100 (MQ=255)                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                      >  2:1059455/1‑100 (MQ=255)                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                     >  2:202933/1‑100 (MQ=255)                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGa                                                  <  2:496025/100‑1 (MQ=255)                                                          taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑Atttt                                            >  2:990914/1‑100 (MQ=255)                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑Attttn                                           <  1:646763/100‑2 (MQ=255)                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATCTTc                                           >  1:680840/1‑100 (MQ=255)                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTc                                           <  1:291255/99‑1 (MQ=255)                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑Gt                                       <  1:1219047/100‑1 (MQ=255)                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑Gt                                       <  2:585107/100‑1 (MQ=255)                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAGACATTGCAGA‑ATTTTC‑‑g                                        >  1:1177669/1‑99 (MQ=255)                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAgg                                    <  2:726402/100‑1 (MQ=255)                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAgg                                    <  1:1304307/100‑1 (MQ=255)                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAgg                                    >  2:1450226/1‑100 (MQ=255)                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGcc                                  >  1:1487946/1‑100 (MQ=255)                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGcc                                  <  2:764765/100‑1 (MQ=255)                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGa                               >  2:980208/1‑100 (MQ=255)                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGa                               <  1:474666/100‑1 (MQ=255)                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGAt                              >  2:937472/1‑100 (MQ=255)                                                                        aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATaa                            <  1:566399/100‑1 (MQ=255)                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGa                       <  2:891177/100‑1 (MQ=255)                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAg                      >  2:390365/1‑100 (MQ=255)                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAg                      >  2:640628/1‑100 (MQ=255)                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:758227/1‑100 (MQ=255)                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:489267/1‑100 (MQ=255)                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGAtgcgctccgtt                    <  2:1081477/100‑11 (MQ=255)                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:1505346/1‑100 (MQ=255)                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:700163/1‑100 (MQ=255)                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAgcg                    <  1:71052/100‑1 (MQ=255)                                                                                 tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAgcgc                   >  2:968875/10‑100 (MQ=18)                                                                                 tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAgcgc                   >  2:1238734/10‑100 (MQ=18)                                                                                 cttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AGCTGC‑‑GTAGGCCTGATAAGCGTAgcgc                   >  1:635269/10‑100 (MQ=18)                                                                                 actttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TTATTG‑‑GTAGGCCTGATAagcgcagcgc                   >  2:1321513/10‑100 (MQ=18)                                                                                  ttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCa                  >  2:1158298/9‑100 (MQ=17)                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑ACCAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:171610/8‑100 (MQ=18)                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  2:26409/1‑100 (MQ=255)                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCa               >  2:818608/6‑100 (MQ=21)                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1493620/5‑100 (MQ=17)                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑ATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg             <  2:150007/97‑1 (MQ=18)                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑ATCCAC‑‑ATAGGCCTGATAAGCGTAGCGCATCAgg             >  2:453122/4‑100 (MQ=21)                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:990134/4‑100 (MQ=21)                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑ATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgga            >  1:150005/3‑99 (MQ=18)                                                                                        aggCCGGATAAGGCGTTTACGCCCCATCCGGCCTAC‑AGAAGCATGCAAATTCAATATATTGCAGA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:856046/3‑100 (MQ=17)                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTGGG‑CTCCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:721670/3‑100 (MQ=21)                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAAG‑ATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:914534/3‑100 (MQ=17)                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGAGATTAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:899645/3‑100 (MQ=17)                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:312695/3‑100 (MQ=21)                                                                                        aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            <  2:591316/100‑1 (MQ=255)                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑ATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1579579/2‑100 (MQ=21)                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:872374/2‑100 (MQ=14)                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  2:121702/2‑100 (MQ=21)                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1109736/1‑99 (MQ=18)                                                                                           ccGGATGCGGCGTGAACGCCTTATCCGGCTTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TTATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc            <  1:888948/97‑1 (MQ=21)                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       >  2:611160/1‑100 (MQ=255)                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    >  2:791490/1‑100 (MQ=255)                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  2:732153/100‑1 (MQ=255)                                                                                                 gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:1157720/100‑1 (MQ=255)                                                                                                 gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:1198682/100‑1 (MQ=255)                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:1198339/1‑100 (MQ=255)                                                                                                  |                                                                                                        GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑ATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070

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: ‑