Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,748	0	C	A	44.4%	84.6 / 48.2	72	E55* (GAG→TAG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (35/5);  new base A (0/32);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.55e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGT  >  NC_000913/3177599‑3177897                                                                                                                                                      |                                                                                                                                                      acTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTc                                                                                                                                                       >  2:135949/1‑150 (MQ=255)       aCCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAAtt                                                                                                                                                <  2:1181159/151‑1 (MQ=255)          ctctTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTcc                                                                                                                                             >  1:165829/1‑151 (MQ=255)              tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                         >  2:1009170/1‑151 (MQ=255)              tCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                         >  1:719566/1‑151 (MQ=255)               cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgccg                                                                                                                                        >  1:921809/1‑151 (MQ=255)               cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                         <  1:704751/150‑1 (MQ=255)               cAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCccgcc                                                                                                                                         <  2:1153708/150‑1 (MQ=255)                aaTCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCcgccgc                                                                                                                                       >  1:638341/1‑151 (MQ=255)                aaTCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCcgccgc                                                                                                                                       >  2:312004/1‑151 (MQ=255)                aaTCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCcgccgc                                                                                                                                       >  2:143629/1‑151 (MQ=255)                      cccGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTc                                                                                                                                  >  2:542055/1‑150 (MQ=255)                        cGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATg                                                                                                                               >  1:120159/1‑151 (MQ=255)                                  ctcCAGTAGCCAAGGGTTTTCGCTGGTGCCGTCCGCGGCAATCCGAACCTGCTCAATCACCACATCTTCATCCCCCACTGGCTCAAAGGGTAGCAAGACGCCGGCGGGACCGCCCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTg                                                                                                                      <  1:310491/150‑1 (MQ=255)                                    ccAGTAGCAAGGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc                                                                                                                    <  1:438530/150‑1 (MQ=255)                                         aGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                               >  1:456444/1‑150 (MQ=255)                                          gCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                               >  2:36532/1‑149 (MQ=255)                                           ccAGGGGGTTTCGCTGTCGTCTTAGCCGGCAATCCAATTCTGCTAAATCAGCACAACTTAAACCCCCACTGGGAAAAAGGGTAGAAAGACGCGGGGGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg                                                                                                               <  1:557052/148‑1 (MQ=255)                                             aaGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAAt                                                                                                            >  2:1020169/1‑149 (MQ=255)                                                    ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                    >  1:1062779/1‑150 (MQ=255)                                                    ttCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATg                                                                                                    >  1:525217/1‑150 (MQ=255)                                                     tCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGAc                                                                                                  >  1:1017136/1‑151 (MQ=255)                                                      cGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                 >  1:1035707/1‑151 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACg                                                                                                 >  2:1079994/1‑150 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTAAAAGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                                <  1:795556/151‑1 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCACCCCGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                                <  1:429082/151‑1 (MQ=255)                                                       gCTGGTGTCGTACGCGGCAATCCAAATCTGCTCAATCAGCACAACTCCATCACGCACTGGGTCAAGGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGa                                                                                                   <  2:772951/148‑1 (MQ=255)                                                           gtgtCAAACGCGGCAATCCAAACCTCCTCAATCACCACACTTCCATCCCGCATTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCCCGTTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGa                                                                                                <  2:359116/147‑1 (MQ=255)                                                              tCGTACGCGGCAATCCGAATCTGCTCAATCAGCACACCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCt                                                                                           <  1:941480/149‑1 (MQ=255)                                                                ggaCGCGGCAATCCAAATCTGCTCAACCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACGGCGGCGTGCCCGCGTTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat                                                                                         <  2:1032654/147‑1 (MQ=255)                                                                    gcgGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATaa                                                                                       >  2:888609/1‑147 (MQ=255)                                                                      ggCAATCCGATTCTGCAAAATCAGCCCACCTTCACCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGa                                                                                    <  2:258730/148‑1 (MQ=255)                                                                       gCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                                 >  1:808897/1‑150 (MQ=255)                                                                        cAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGTTCAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTa                                                                                <  1:724247/150‑1 (MQ=255)                                                                         aaTCCGAATCTGCTCAATCAGCACAACTTCACCACGCACGGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt                                                                                 <  2:403664/148‑1 (MQ=255)                                                                            ccaaaTCTGCTCAACCACCACAACTTCACCCCGCCCGGGGTCAAAGGGTAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa                                                                            <  1:1108292/147‑1 (MQ=255)                                                                               attCTGCCCAATCCGCACACTTTCACCACGCCCGGGGTAAAGGGGTAGCAAGCCTGCGTCGTCCCCCCTCCAAAAAATTACCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaa                                                                           <  2:194951/146‑1 (MQ=255)                                                                               aaTCTTCTAAATCAGCACACCTGCTTCCCGCACGGGGTAAAAGGGTACAAAGCCGGCGGCGTAACCGCGCTAAAAAATTTGCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                         <  2:748944/150‑1 (MQ=255)                                                                               aaTCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaa                                                                         >  2:283466/1‑150 (MQ=255)                                                                               aaTCTGCTCAATCAGCAAAACTTCATCACGCACGGGGTAAAAGGGTAGCAAGACGGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                        <  1:412850/151‑1 (MQ=255)                                                                               aaTCTGCTCAACCAGCACACCTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaa                                                                        <  1:1041607/151‑1 (MQ=255)                                                                                   tGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGc                                                                    >  2:739947/1‑151 (MQ=255)                                                                                   tGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGc                                                                    >  1:426036/1‑151 (MQ=255)                                                                                       cAACCATCACAACTTCTCCACGCACGGGGTCAAAGGGTAGCAAGCCTGGGGCGTGACCCCTCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcg                                                                  <  2:241590/149‑1 (MQ=255)                                                                                       cAACCAGCACAACTCCACCCCGCACTGGGTCAAAGGGTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCcgcg                                                                <  2:1144514/151‑1 (MQ=255)                                                                                         aTCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt                                                              >  2:258755/1‑151 (MQ=255)                                                                                           cAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATa                                                             >  2:119534/1‑150 (MQ=255)                                                                                            aGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                            >  1:1107885/1‑150 (MQ=255)                                                                                             gCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaa                                                            >  1:201966/1‑149 (MQ=255)                                                                                              cacaCCTTCATCCCGCACTGGGTAAAAGGGTACCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATaaa                                                           <  2:144175/149‑1 (MQ=255)                                                                                                     tcatcaCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct                                                   >  1:644552/1‑150 (MQ=255)                                                                                                       atcaCGCACTGGGTAAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTtct                                                   <  1:275927/148‑1 (MQ=255)                                                                                                         cccgaccggggTAAAGGGGTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                              <  2:85472/143‑1 (MQ=255)                                                                                                         cacgcacgGGGTCAAAGGGTAGCAAGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTg                                               <  1:545446/150‑1 (MQ=255)                                                                                                         caCGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGc                                              >  1:565392/1‑151 (MQ=255)                                                                                                           cGCACTGGGTCAAAGGGTAGAAAGACTGCGGCGTGACCGCGCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGCGTTTCTCGt                                                <  1:231314/147‑1 (MQ=255)                                                                                                             cACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaat                                           >  2:370917/1‑150 (MQ=255)                                                                                                               cTGGGCCAAGGGGTGAAAAGACGGCGGCGTGCCCGCGCTAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata                                          <  2:154235/149‑1 (MQ=255)                                                                                                                 gggTCAAAGGGTAGCAAGCCTGGGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAGGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa                                         <  2:412356/148‑1 (MQ=255)                                                                                                                 gggTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAtt                                       >  2:1159687/1‑150 (MQ=255)                                                                                                                  gttCAAAGGGTAACAACCCTGGGGAGTGACCCCCCTAAAAATTTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt                                      <  2:1054589/148‑1 (MQ=255)                                                                                                                      aaggggTAGCAAGACTGGGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTa                                   <  1:906384/146‑1 (MQ=255)                                                                                                                      aaaGGGTAGCAAGCCGGCGGCGTGCCCGCGCTAAAAATTTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACa                                 <  2:712195/151‑1 (MQ=255)                                                                                                                      aaaGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTAc                                  >  2:735589/1‑150 (MQ=255)                                                                                                                         gggTAGCAAGCCTGCGGCGTGACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATc                               <  1:722670/150‑1 (MQ=255)                                                                                                                             aGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGtttt                          >  2:691225/1‑151 (MQ=255)                                                                                                                                aaGACTGCGGCGTGCCCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTCTCTCGTGCAATAATTTCTACATCGttt                           <  2:235578/147‑1 (MQ=255)                                                                                                                                 aGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGcc                      >  1:318705/1‑151 (MQ=255)                                                                                                                                    cTGGGGCTTGACCGCGCAAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGttttt                         <  2:273241/145‑1 (MQ=255)                                                                                                                                    cTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaa                    >  2:732933/1‑150 (MQ=255)                                                                                                                                           gggACCGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg            <  2:811969/149‑1 (MQ=255)                                                                                                                                            tgcccGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACg            <  1:656336/147‑1 (MQ=255)                                                                                                                                                     tCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGt  >  2:125321/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAAGGGTAGCAAGACTGCGGCGTGACCGCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGGGCAGGTTGT  >  NC_000913/3177599‑3177897

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑