Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,691,325	G→C	5.7%	Q91H (CAG→CAC)	PP_1487 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,691,325	0	G	C	5.7%	222.0 / 4.5	70	Q91H (CAG→CAC)	PP_1487	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/27);  new base C (2/2);  total (41/29)
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.40e-01

CGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGT  >  NC_002947/1691181‑1691472                                                                                                                                                   |                                                                                                                                                    cgtccTGCTGGGTGGGGTCGGTGTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCATGTTCGGTGGCAAGGCCCTGACTTACGTCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCggg                                                                                                                                                 <  4:131452/146‑1 (MQ=255)  tgcgtgctgGGTGGGTGCGGTCTGGGCCTGTTGCTGTGTTGCAAGATGGCCCTCATTTAGGGTTGCATGGCACTGACTTACGTCGGC‑‑CTCGCCGCGCTGTGGGTACTGGCCTATAAAGCTTACGGCAACTTGCAGGCAAACCTGCtggg                                                                                                                                                <  2:152872/145‑1 (MQ=255)      tgctgGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGcc                                                                                                                                            >  5:9669/1‑149 (MQ=255)      tgctgGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGcc                                                                                                                                            >  6:104976/1‑149 (MQ=255)                      tctGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCaaa                                                                                                                            >  3:104511/1‑149 (MQ=255)                       ctGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCGGATCGCCGCGCTGGGGGTACTGGCCTATAAAGATTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCa                                                                                                                              >  1:315052/1‑148 (MQ=255)                        tGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAAc                                                                                                                           >  3:140084/1‑148 (MQ=255)                              tGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATc                                                                                                                       <  8:159465/146‑1 (MQ=255)                              tGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATc                                                                                                                       <  8:159463/146‑1 (MQ=255)                              tGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATc                                                                                                                       >  7:159463/1‑146 (MQ=255)                              tGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAAACATc                                                                                                                       >  7:159465/1‑146 (MQ=255)                               gTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGAcc                                                                                                                   >  3:253981/1‑149 (MQ=255)                                         aGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTgccgcc                                                                                                         >  2:275024/1‑149 (MQ=255)                                          gCAAGAAGGCCGGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CCCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGccgccg                                                                                                        <  3:98533/149‑1 (MQ=255)                                            aagaagGCCCGCAATTACGTTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGc                                                                                                      <  6:188573/149‑1 (MQ=255)                                             agaagGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGcc                                                                                                     >  7:327562/1‑149 (MQ=255)                                                agGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAg                                                                                                  >  2:146073/1‑149 (MQ=255)                                                    ccGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCACACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCTTCGACCGCCTGCCGCCGGCCCAGGCCg                                                                                              >  5:54065/1‑149 (MQ=255)                                                    ccGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGcc                                                                                               >  2:19273/1‑148 (MQ=255)                                                     cGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGa                                                                                             <  8:327562/149‑1 (MQ=255)                                                                tGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTACCGCCGGCCCAGGCCGAGCAGCACAGcc                                                                                  <  6:9669/149‑1 (MQ=255)                                                                  gCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCACCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGCGCAGCACAGCCAg                                                                                >  4:31225/1‑149 (MQ=255)                                                                  gCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTAGAAAGCTTACGGCAACTGGCAGGCAAACCACCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAg                                                                                >  4:36038/1‑149 (MQ=255)                                                                    aaGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAACCCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGCCCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGc                                                                              >  4:215327/1‑149 (MQ=255)                                                                     aGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAgg                                                                               <  5:67112/147‑1 (MQ=255)                                                                     aGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAgg                                                                               >  6:67112/1‑147 (MQ=255)                                                                         ccTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAg                                                                                >  7:233892/1‑142 (MQ=255)                                                                         ccTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAg                                                                                <  8:233892/142‑1 (MQ=255)                                                                          cTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACt                                                                        >  5:31030/1‑149 (MQ=255)                                                                               ttACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTgcgcg                                                                   >  7:268699/1‑149 (MQ=255)                                                                                 aCGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGcc                                                                 >  6:166638/1‑149 (MQ=255)                                                                                       gc‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCtggtg                                                           >  2:256869/1‑149 (MQ=255)                                                                                       gc‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGG‑CCAGGCCGAGCATCACAGCCAGGCGGTACTGCGCGCCCtggtgg                                                          <  2:315052/149‑1 (MQ=255)                                                                                           cTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggc                                                         <  4:140084/149‑1 (MQ=255)                                                                                           cTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggc                                                         >  7:190991/1‑149 (MQ=255)                                                                                            tCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcg                                                        >  2:233766/1‑149 (MQ=255)                                                                                              tccgcGAGGGGGGTACGGGCCTATAANTCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcggc                                                      <  4:163227/148‑1 (MQ=255)                                                                                              gccgcGGGGGGGGTACGGGCCTATACAGCGTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTggcggc                                                      <  4:157203/149‑1 (MQ=255)                                                                                              gccgcGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGTGCGCCCTGGTggcg                                                        <  1:19273/147‑1 (MQ=255)                                                                                                  cgcTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGcc                                                  <  5:166638/149‑1 (MQ=255)                                                                                                   gcgcggggTCCTGGCCTTTAAAGGTTACGGCAACGGGCTGGCAACCCCGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCa                                                 <  7:57467/145‑1 (MQ=255)                                                                                                    cggggggTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCCGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCaa                                                <  5:104976/147‑1 (MQ=255)                                                                                                           tACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGAt                                         <  8:268699/149‑1 (MQ=255)                                                                                                                   tataAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCatat                                 >  1:268021/1‑149 (MQ=255)                                                                                                                    ataAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATc                                >  2:357340/1‑149 (MQ=255)                                                                                                                     taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATggg                                      >  1:109053/1‑142 (MQ=255)                                                                                                                     taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATggg                                      <  2:109053/142‑1 (MQ=255)                                                                                                                     taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATc                                >  3:257796/1‑148 (MQ=255)                                                                                                                     taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATc                                >  4:333937/1‑148 (MQ=255)                                                                                                                     taAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGCGCCGCAGCCGATCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGCCCGCGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATcg                               >  7:215632/1‑149 (MQ=255)                                                                                                                       caGCTTACGGCAAGTGGCAGGCAAACCAGGGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATcgac                             <  8:219204/148‑1 (MQ=255)                                                                                                                               ggCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGAcgagcgcga                     <  1:256869/149‑1 (MQ=255)                                                                                                                               ggCAACTGGCAGGAAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGAcgagcgc                       >  8:170559/1‑147 (MQ=255)                                                                                                                                gCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGAcgagcgcgag                    >  4:12874/1‑149 (MQ=255)                                                                                                                                   aCTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCgg                 <  4:253981/149‑1 (MQ=255)                                                                                                                                     tGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgc               >  4:136706/1‑149 (MQ=255)                                                                                                                                       gcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcgc             >  3:340240/1‑149 (MQ=255)                                                                                                                                       gcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcgc             <  6:31030/149‑1 (MQ=255)                                                                                                                                       gcaggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGgcgc             >  8:309412/1‑149 (MQ=255)                                                                                                                                         aggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCt            >  6:276943/1‑148 (MQ=255)                                                                                                                                         aggcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCt            >  3:16116/1‑148 (MQ=255)                                                                                                                                           gcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGAt         >  8:30711/1‑149 (MQ=255)                                                                                                                                           gcaAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGAt         >  3:360168/1‑149 (MQ=255)                                                                                                                                              aaCCTGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGa      <  1:277067/149‑1 (MQ=255)                                                                                                                                               aCCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGCCGGCCAACTCCGATCGGCAAATCGACGAGCGCGCGCGGGCGCTGATCGAg     >  7:359208/1‑149 (MQ=255)                                                                                                                                                ccAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAgg    <  4:104511/149‑1 (MQ=255)                                                                                                                                                ccACCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAgg    <  3:31225/149‑1 (MQ=255)                                                                                                                                                ccACCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAgg    <  3:36038/149‑1 (MQ=255)                                                                                                                                                 cAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAggg   >  4:320842/1‑149 (MQ=255)                                                                                                                                                  aGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGt  <  1:146073/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGC‑‑CTCGCCGCGCTGGGGGTACTGGCCTATAAAGCTTACGGCAACTGGCAGGCAAACCAGCGGGGGGCCGCAGCCGAGCCTCAAACCATCGACCGCCTGCCGCCGGCCCAGGCCGAGCAGCACAGCCAGGCGGTACTGCGCGCCCTGGTGGCGGCGGCCAAGTCCGATGGGCATATCGACGAGCGCGAGCGGGCGCTGATCGAGGGT  >  NC_002947/1691181‑1691472

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