Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,655,155	G→C	5.4%	R402G (CGT→GGT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,155	0	G	C	5.4%	247.1 / 4.6	74	R402G (CGT→GGT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base C (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGG  >  NC_000913/2655007‑2655300                                                                                                                                                     |                                                                                                                                                  cGCAATTTGCTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTGCAACCGGCGCGTTACg                                                                                                                                                   <  1:821967/149‑1 (MQ=255)      tttGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAg                                                                                                                                              <  2:487520/149‑1 (MQ=255)         gttgAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCg                                                                                                                                           >  2:1577948/1‑149 (MQ=255)            gAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCg                                                                                                                                        >  1:371548/1‑149 (MQ=255)             aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCGTATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGa                                                                                                                                       >  1:189024/1‑149 (MQ=255)             aaTTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAGCCGGCGCTTTACGGTAAGTCGCCGa                                                                                                                                       >  1:1106325/1‑149 (MQ=255)                ttGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCa                                                                                                                                    <  2:1245601/149‑1 (MQ=255)                  gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                                                                                                                                  <  1:867308/149‑1 (MQ=255)                  gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                                                                                                                                  >  2:1275518/1‑149 (MQ=255)                  gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGACTTACCGTAAAGCGCCggttgaa                                                                                                                                  <  1:531311/149‑7 (MQ=38)                    aGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTATGTCGCCGAGCAGTCg                                                                                                                                >  2:307627/1‑149 (MQ=255)                     gAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGCTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGa                                                                                                                               >  1:550470/1‑149 (MQ=255)                      aaTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGat                                                                                                                              <  1:1470589/149‑1 (MQ=255)                           gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                                                                                                                         <  2:1380910/149‑1 (MQ=255)                           gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                                                                                                                         <  1:364948/149‑1 (MQ=255)                              tgCCGGATGCGGCGTAAACGCCTTATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGGCGATATGcagcc                                                                                                                      >  1:37033/2‑149 (MQ=255)                              ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcc                                                                                                                      >  1:741629/1‑149 (MQ=255)                              ggCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcc                                                                                                                      >  1:733762/1‑149 (MQ=255)                               gCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGcagcca                                                                                                                     >  2:800136/1‑149 (MQ=255)                                   gATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCt                                                                                                                 >  2:849021/1‑149 (MQ=255)                                      aaGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGc                                                                                                              <  1:839315/149‑1 (MQ=255)                                       aGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCTCACCAAGTCCCGTAGCGCCCGCAGACGACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCt                                                                                                             <  1:595325/149‑1 (MQ=255)                                         gCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGa                                                                                                           >  2:996287/1‑149 (MQ=255)                                           gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                         <  1:1546332/149‑1 (MQ=255)                                           gTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                         <  1:1376122/149‑1 (MQ=255)                                           gTTTACGACGCATCAGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATa                                                                                                         <  2:447303/149‑1 (MQ=255)                                               aCGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTtc                                                                                                     <  2:427972/149‑1 (MQ=255)                                                  ccgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTca                                                                                                  <  2:1243531/149‑1 (MQ=255)                                                   cgcATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGCCCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTcaa                                                                                                 >  1:1302744/1‑149 (MQ=255)                                                          ggCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTg                                                                                          >  1:424929/1‑149 (MQ=255)                                                          ggCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTg                                                                                          >  2:652632/1‑149 (MQ=255)                                                           gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                         <  2:495990/149‑1 (MQ=255)                                                           gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                         <  2:53971/149‑1 (MQ=255)                                                               atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                                   >  2:1052692/1‑87 (MQ=255)                                                               atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                                   <  1:1052692/87‑1 (MQ=255)                                                                ttAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAg                                                                                    <  1:1571432/149‑1 (MQ=255)                                                                ttAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTGCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAg                                                                                    <  1:267943/149‑1 (MQ=255)                                                                   gTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGaa                                                                                 <  1:643745/149‑1 (MQ=255)                                                                    tttACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAg                                                                                <  2:1210740/149‑1 (MQ=255)                                                                          ccGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCccgcc                                                                          >  1:826291/1‑149 (MQ=255)                                                                             ttAACAGCTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCccgccgcc                                                                       >  1:73887/1‑149 (MQ=255)                                                                             ttAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCGTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACCGGAAGCccgccgcc                                                                       >  2:1406353/1‑149 (MQ=255)                                                                                 cAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGc                                                                   >  1:847576/1‑149 (MQ=255)                                                                                    aTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAGGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTcg                                                                >  1:5657/1‑149 (MQ=255)                                                                                    aTTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTcg                                                                >  2:687771/1‑149 (MQ=255)                                                                                       aGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCgcgc                                                             <  2:937517/149‑1 (MQ=255)                                                                                         cTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCt                                                           <  1:425245/149‑1 (MQ=255)                                                                                         cTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCt                                                           <  1:540268/149‑1 (MQ=255)                                                                                          tATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTg                                                          >  2:425446/1‑149 (MQ=255)                                                                                            tCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTTAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGcc                                                        >  2:1455737/1‑149 (MQ=255)                                                                                            tCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGcc                                                        >  2:1449263/1‑149 (MQ=255)                                                                                                      acCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATAcgccg                                              >  2:1550084/1‑149 (MQ=255)                                                                                                      acCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATAcgccg                                              >  1:41804/1‑149 (MQ=255)                                                                                                            tCCCGTAGCGCCAGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTc                                        <  1:1070147/149‑1 (MQ=255)                                                                                                             cccGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTcc                                       <  1:651707/149‑1 (MQ=255)                                                                                                              ccGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTGGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCg                                      >  2:1288944/1‑149 (MQ=255)                                                                                                                 tAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTa                                   <  1:1566736/149‑1 (MQ=255)                                                                                                                 tAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTa                                   <  1:652632/149‑1 (MQ=255)                                                                                                                       ccGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTc                             <  1:474528/149‑1 (MQ=255)                                                                                                                       ccGCAGACCACTGTTCAACCCGCGCTTTACGGTAAGTCGCCGAGCAGTCGATCTGCAGCCAGCCTTGCTGATAGGTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTc                             <  1:1468420/149‑1 (MQ=255)                                                                                                                        cGCAGACCACTGTTCAGCCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCa                            <  2:1106325/149‑1 (MQ=255)                                                                                                                         gCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAg                           <  1:849021/149‑1 (MQ=255)                                                                                                                               cACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGc                     >  1:189127/1‑149 (MQ=255)                                                                                                                                  tGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCaaa                  <  1:1577948/149‑1 (MQ=255)                                                                                                                                  tGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCaaa                  <  1:1550084/149‑1 (MQ=255)                                                                                                                                       aaCCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAg             >  2:506064/1‑149 (MQ=255)                                                                                                                                        ctcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGa            >  1:958803/3‑149 (MQ=255)                                                                                                                                         tcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc           >  2:541115/2‑149 (MQ=255)                                                                                                                                         tcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc           <  2:958803/148‑1 (MQ=255)                                                                                                                                         ccGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc           <  2:733762/149‑1 (MQ=255)                                                                                                                                                tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCgg                      >  1:260148/1‑131 (MQ=255)                                                                                                                                                tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCgg                      <  2:260148/131‑1 (MQ=255)                                                                                                                                                tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGCCGGCAgct    >  1:54287/1‑149 (MQ=255)                                                                                                                                                tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgct    >  1:192302/1‑149 (MQ=255)                                                                                                                                                tttACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgct    >  2:658051/1‑149 (MQ=255)                                                                                                                                                  tACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgctgg  <  2:550470/149‑1 (MQ=255)                                                                                                                                                  tACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAgctgg  <  2:41804/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGG  >  NC_000913/2655007‑2655300

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