Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,655,161	T→G	5.2%	T400P (ACT→CCT)	pepB ←	aminopeptidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,655,161	0	T	G	5.2%	221.7 / 4.3	78	T400P (ACT→CCT)	pepB	aminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (34/40);  new base G (2/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCGGT  >  NC_000913/2655015‑2655307                                                                                                                                                   |                                                                                                                                                   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                                                                                                                                <  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                                                                                                 >  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)                                                        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: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:958803/148‑1 (MQ=255)                                                                                                                                 tcGGCGACTTACCGTAAAGCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGAc                  >  2:541115/2‑149 (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)                                                                                                                                              gTAAGTCCCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGc     <  1:800136/149‑1 (MQ=255)                                                                                                                                                aaGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGTTACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCgg   <  1:1449263/149‑1 (MQ=255)                                                                                                                                                aaGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCgg   <  2:561411/149‑1 (MQ=255)                                                                                                                                                aaGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCGCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCgg   <  1:1455737/149‑1 (MQ=255)                                                                                                                                                 aGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCGCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCGGt  <  2:298511/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTTGAATTTGCAGAATTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATATTAGTTTACGCCGTTAACAGATTAGCTATCGTGCGCACACCAAGTCCCGTAGCGCCCGCAGACCACTGTTCAACCGGCGCTTTACGGTAAGTCGCCGAGCAGTCGATATGCAGCCAGCCTTGCTGATAGTTCTCAACAAAGTGCGACAGGAAGCCCGCCGCCGTGCTCGCGCCTGCCGGATACGCCGCGCTTCCGGTATTGTTCAGTTCGGCAAAGTTAGACGGCAGCTGGCTGCGGT  >  NC_000913/2655015‑2655307

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