Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,621,330	G→C	6.1%	D49H (GAT→CAT) ‡	PP_2294 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,621,330	0	G	C	6.1%	204.5 / 3.4	66	D49H (GAT→CAT) ‡	PP_2294	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base C (2/2);  total (32/34)
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.56e-01

ATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGG  >  NC_002947/2621186‑2621473                                                                                                                                                 |                                                                                                                                                atgatgTCGAGACGGGCGATTGGGATCGTGTTAGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGc                                                                                                                                             >  1:245919/1‑149 (MQ=255)  tgctgTCGAGACGGGCGCTTGCGATCGTGCTGGCCATTGGAGCTGCTTACGTTGTCCTCAAAGGAGTTGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGGTCAGCTTGATGCt                                                                                                                                            <  4:109588/146‑1 (MQ=255)       tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGCGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGt                                                                                                                                       <  1:217315/149‑1 (MQ=255)       tCGAGACGGGCGAGTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgca                                                                                                                                       <  6:79030/149‑3 (MQ=255)         gagaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGtt                                                                                                                                      <  1:68038/148‑1 (MQ=255)          agaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTgg                                                                                                                                    <  8:267688/149‑1 (MQ=255)                       ggATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCaaa                                                                                                                       <  4:80146/149‑1 (MQ=255)                                 agCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGAt                                                                                                              <  2:245919/147‑1 (MQ=255)                                  gCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGAt                                                                                                              >  6:311434/1‑147 (MQ=255)                                  gCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCa                                                                                                            >  7:352049/1‑149 (MQ=255)                                   ccATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcat                                                                                                                                      >  2:109347/1‑119 (MQ=255)                                   ccATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcat                                                                                                                                      <  1:109347/122‑4 (MQ=255)                                    cATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCag                                                                                                           >  7:182331/1‑148 (MQ=255)                                         gAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAgaagaa                                                                                                      <  8:352049/148‑1 (MQ=255)                                            ctgctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                  <  5:271665/149‑1 (MQ=255)                                                ttATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACa                                                                                              >  6:26675/1‑149 (MQ=255)                                                  aTGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACAtt                                                                                            >  2:307952/1‑149 (MQ=255)                                                     ttgtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTc                                                                                         >  6:285603/1‑149 (MQ=255)                                                      tgtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTc                                                                                         >  3:239169/1‑148 (MQ=255)                                                        tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTg                                                                                                                <  6:357462/123‑1 (MQ=255)                                                        tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTg                                                                                                                >  5:357462/1‑123 (MQ=255)                                                        tCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGa                                                                                       <  7:109862/148‑1 (MQ=255)                                                          cTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAATAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTc                                                                                    <  8:182331/149‑1 (MQ=255)                                                             aaaGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCctc                                                                                 >  4:318962/1‑149 (MQ=255)                                                             aaaGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCctc                                                                                 >  8:322407/1‑149 (MQ=255)                                                                ggAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGAGATTTTCGATTCCTCtt                                                                               >  6:358622/1‑148 (MQ=255)                                                                      gCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACaaaaa                                                                        >  1:145277/1‑149 (MQ=255)                                                                       cTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAAt                                                                       <  3:318962/149‑1 (MQ=255)                                                                           cgcgCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACt                                                                   <  4:137034/149‑1 (MQ=255)                                                                              gcCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGa                                                                <  2:203768/149‑1 (MQ=255)                                                                                cTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGAt                                                               >  6:300025/1‑148 (MQ=255)                                                                                      cttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                         >  3:283688/1‑148 (MQ=255)                                                                                       ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCTCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                         <  7:322407/147‑1 (MQ=255)                                                                                       ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATg                                                        >  2:132824/1‑148 (MQ=255)                                                                                           tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAgg                                                   >  1:194876/1‑149 (MQ=255)                                                                                           tcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAgg                                                   >  3:302776/1‑149 (MQ=255)                                                                                                 aGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAAt                                             >  2:69834/1‑149 (MQ=255)                                                                                                     cTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAc                                         <  5:300025/149‑1 (MQ=255)                                                                                                      tCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAcc                                        >  7:261214/1‑149 (MQ=255)                                                                                                      tCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAcc                                        >  4:369902/1‑149 (MQ=255)                                                                                                       ccAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCt                                       >  1:205936/1‑149 (MQ=255)                                                                                                        cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAAt                                             >  4:369888/1‑142 (MQ=255)                                                                                                        cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAAt                                             <  3:369888/142‑1 (MQ=255)                                                                                                        cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAcc                                        <  8:261214/147‑1 (MQ=255)                                                                                                        cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                      >  3:98686/1‑149 (MQ=255)                                                                                                        cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                      >  5:365010/1‑149 (MQ=255)                                                                                                          gagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCGTGAAAGGCTCAATTGACCTtgt                                    <  4:239169/149‑1 (MQ=255)                                                                                                           agGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgtg                                   <  5:26675/149‑1 (MQ=255)                                                                                                               ttATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACt                                                                   >  3:62502/1‑113 (MQ=255)                                                                                                               ttATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACt                                                                   <  4:62502/113‑1 (MQ=255)                                                                                                                tATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAtt                               <  5:311434/148‑1 (MQ=255)                                                                                                                 aTCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTc                              <  1:132824/148‑1 (MQ=255)                                                                                                                  tCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAg                            <  5:285603/149‑1 (MQ=255)                                                                                                                       aaaGCGCATTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAAc                       <  4:302776/149‑1 (MQ=255)                                                                                                                       aaaGCGCAGTGACCTATGCTCAGCTTGATGCTTAGGTTGGCTGCAAGATCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAAc                       <  5:358622/149‑1 (MQ=255)                                                                                                                            gcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTTGATGACCATGAAAGCCTCAATTGACCTTGTGGATTCAGACAACGTATc                  >  2:245018/1‑149 (MQ=255)                                                                                                                            gcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATc                  <  2:194876/149‑1 (MQ=255)                                                                                                                                  aCCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCaa            <  3:369902/149‑1 (MQ=255)                                                                                                                                     tATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                    >  1:67586/1‑122 (MQ=255)                                                                                                                                     tATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                    <  2:67586/122‑1 (MQ=255)                                                                                                                                        catCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGa      >  5:265514/3‑149 (MQ=255)                                                                                                                                         cTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGAc     <  1:69834/149‑1 (MQ=255)                                                                                                                                            aGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACGAAAATCACTGGATGACCATGAAAg                                                    <  4:341682/99‑1 (MQ=255)                                                                                                                                            aGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACGAAAATCACTGGATGACCATGAAAg                                                    >  3:341682/1‑99 (MQ=255)                                                                                                                                            aGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACg    <  4:283688/147‑1 (MQ=255)                                                                                                                                             gCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACggg  >  8:156916/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                ATGATGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGG  >  NC_002947/2621186‑2621473

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