Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,655,149	C→A	5.5%	A404S (GCG→TCG)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,149	0	C	A	5.5%	228.0 / 2.1	72	A404S (GCG→TCG)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base C (31/37);  new base A (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAC  >  NC_000913/2655005‑2655291                                                                                                                                                 |                                                                                                                                               ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                                                                                                            <  1:1025743/149‑1 (MQ=255) ttCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTa                                                                                                                                            <  2:1128951/149‑1 (MQ=255)   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                                                                                                                         >  2:1275518/1‑149 (MQ=255)                    gCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGt                                                                                                                         <  1:867308/149‑1 (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                                                                                                                <  1:364948/149‑1 (MQ=255)                             gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGca                                                                                                                <  2:1380910/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                                                                                 >  2:652632/1‑149 (MQ=255)                                                            ggCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTg                                                                                 >  1:424929/1‑149 (MQ=255)                                                             gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                <  2:53971/149‑1 (MQ=255)                                                             gCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGc                                                                                <  2:495990/149‑1 (MQ=255)                                                                 atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                          <  1:1052692/87‑1 (MQ=255)                                                                 atTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACg                                                                                                                                          >  2:1052692/1‑87 (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:540268/149‑1 (MQ=255)                                                                                           cTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCt                                                  <  1:425245/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                                     >  1:41804/1‑149 (MQ=255)                                                                                                        acCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATAcgccg                                     >  2:1550084/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:1550084/149‑1 (MQ=255)                                                                                                                                    tGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCaaa         <  1:1577948/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)                                                                                                                                                 |                                                                                                                                               TTCGCAATTTGTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAC  >  NC_000913/2655005‑2655291

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