Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,710,293 T→A 5.6% S63C (AGC→TGC)  PP_4168 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,710,2930TA 5.6% 218.4 / 2.7 72S63C (AGC→TGC) PP_4168hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base A (2/2);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACCCGCAGGAT  >  NC_002947/4710145‑4710440
                                                                                                                                                    |                                                                                                                                                    
tCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGAGGCCGATCTGTGCCAGCTACTGGCt                                                                                                                                                      >  1:253627/1‑149 (MQ=255)
   tCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                                      >  1:227283/1‑146 (MQ=255)
   tCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                                      >  4:407992/1‑146 (MQ=255)
          gCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGAGCCAGTACCTGGCACAGatcg                                                                                                                                               >  3:66838/1‑143 (MQ=255)
          gCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGAGCCAGTACCTGGCACAGatcg                                                                                                                                               <  4:66838/146‑4 (MQ=255)
                gATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAg                                                                                                                                     <  5:260781/149‑1 (MQ=255)
                                 aCGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCt                                                                                                                    >  4:173786/1‑149 (MQ=255)
                                 aCGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCt                                                                                                                    <  3:173786/149‑1 (MQ=255)
                                       tGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                                                                                                                                                  >  5:121767/1‑114 (MQ=255)
                                       tGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                                                                                                                                                  <  6:121767/114‑1 (MQ=255)
                                         cgcTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGgcgcgc                                                                                                            <  3:407992/149‑1 (MQ=255)
                                          gctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTg                                                                                                                   <  2:120367/141‑1 (MQ=255)
                                          gctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTg                                                                                                                   >  1:120367/1‑141 (MQ=255)
                                           ctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                          >  3:247141/1‑149 (MQ=255)
                                              ctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTtc                                                                                                                                  >  7:101691/1‑122 (MQ=255)
                                              ctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTtc                                                                                                                                  <  8:101691/122‑1 (MQ=255)
                                              ctTCACATTCCGGCGCAGGCGGAGCCGGCCAGTCTGCGAACGTCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGTCTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATac                                                                                                       <  1:284662/149‑1 (MQ=255)
                                                tCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGgcgcgc                                                                                                            >  1:161207/1‑142 (MQ=255)
                                                  aCCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATacac                                                                                                     >  6:203803/1‑147 (MQ=255)
                                                    cTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGt                                                                                                 <  4:247141/149‑1 (MQ=255)
                                                    cTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGt                                                                                                 >  3:309400/1‑149 (MQ=255)
                                                     ttCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                                >  2:253624/1‑149 (MQ=255)
                                                        cGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCgg                                                                                              <  6:244058/148‑1 (MQ=255)
                                                         ggCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGGCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTa                                                                                            <  8:303670/149‑1 (MQ=255)
                                                          gcgcACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCgg                                                                                              >  3:347840/1‑146 (MQ=255)
                                                             cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGAt                                                                                         >  6:75229/1‑148 (MQ=255)
                                                             cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGAt                                                                                         >  8:69814/1‑148 (MQ=255)
                                                             cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        <  1:253624/149‑1 (MQ=255)
                                                              aCGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCa                                                                                       >  8:148793/1‑149 (MQ=255)
                                                                 cgGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCt                                                                                    <  5:75229/149‑1 (MQ=255)
                                                                    ggccggccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTa                                                                                 <  4:309400/149‑1 (MQ=255)
                                                                     gccggccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTAc                                                                                >  5:128608/1‑149 (MQ=255)
                                                                          ccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                                      <  1:326976/75‑1 (MQ=255)
                                                                          ccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                                      >  2:326976/1‑75 (MQ=255)
                                                                              tCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAg                                                                       <  5:203803/149‑1 (MQ=255)
                                                                                tGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        >  3:372602/1‑130 (MQ=255)
                                                                                tGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        <  4:372602/130‑1 (MQ=255)
                                                                                 gCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGc                                                                     >  8:128697/1‑148 (MQ=255)
                                                                                  cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                          >  5:371134/1‑110 (MQ=255)
                                                                                  cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                          <  6:371134/110‑1 (MQ=255)
                                                                                  cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        >  4:305314/1‑128 (MQ=255)
                                                                                  cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        <  3:305314/128‑1 (MQ=255)
                                                                                           aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCAccc                                                                <  6:22165/143‑1 (MQ=255)
                                                                                           aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCAccc                                                                >  5:22165/1‑143 (MQ=255)
                                                                                           aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAg                                                           >  3:77903/1‑148 (MQ=255)
                                                                                                ttCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                                >  8:359446/1‑106 (MQ=255)
                                                                                                ttCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                                <  7:359446/106‑1 (MQ=255)
                                                                                                     gTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGc                                                  >  8:106940/1‑147 (MQ=255)
                                                                                                      tCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAg                                               >  2:222918/1‑149 (MQ=255)
                                                                                                         cTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        >  8:336496/1‑105 (MQ=255)
                                                                                                         cTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                        <  7:336496/105‑1 (MQ=255)
                                                                                                            tgGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGTCTCGCGCTGGTCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCAACCCAACAGCACCTCTGCCAGATGc                                           >  5:113906/1‑147 (MQ=255)
                                                                                                            tgGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGc                                           >  5:358931/1‑147 (MQ=255)
                                                                                                              gAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTGGCCGAGTTGctgt                                        <  2:227283/148‑4 (MQ=255)
                                                                                                               aaaCTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGac                                      <  7:148793/149‑1 (MQ=255)
                                                                                                                aaCTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGacc                                     <  7:69814/149‑1 (MQ=255)
                                                                                                                     ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCAGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccac                                   >  1:290767/1‑147 (MQ=255)
                                                                                                                     ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCaa                                                              >  3:188458/1‑119 (MQ=255)
                                                                                                                     ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCaa                                                              <  4:188458/119‑1 (MQ=255)
                                                                                                                             gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccacc                                  >  1:167839/1‑139 (MQ=255)
                                                                                                                             gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccacc                                  <  2:167839/139‑1 (MQ=255)
                                                                                                                             gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCa                        <  4:347840/149‑1 (MQ=255)
                                                                                                                              ccGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAg                       <  7:106940/149‑1 (MQ=255)
                                                                                                                               cGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAAATCAGCTGTACCAGCACCAGGCCACCCAACAGCaa                                                        <  6:350710/115‑1 (MQ=255)
                                                                                                                               cGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAAATCAGCTGTACCAGCACCAGGCCACCCAACAGCaa                                                        >  5:350710/1‑115 (MQ=255)
                                                                                                                                  tCTGAGCCAGTACCTGGCACAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCg                                                    >  7:358767/1‑116 (MQ=255)
                                                                                                                                  tCTGAGCCAGTACCTGGCACAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCg                                                    <  8:358767/116‑1 (MQ=255)
                                                                                                                                     gtgCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGtt                                             >  8:293119/1‑120 (MQ=255)
                                                                                                                                     gtgCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGtt                                             <  7:293119/120‑1 (MQ=255)
                                                                                                                                           ggTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCAcc          <  6:358931/149‑1 (MQ=255)
                                                                                                                                             tACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGc                                                          <  1:379728/99‑1 (MQ=255)
                                                                                                                                             tACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGc                                                          >  2:379728/1‑99 (MQ=255)
                                                                                                                                                tGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACGAGTCGGTAGATCAGCTGTAACAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACCCGCAg     >  6:101929/1‑149 (MQ=255)
                                                                                                                                                  gCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGc                      >  7:37769/1‑130 (MQ=255)
                                                                                                                                                  gCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGc                      <  8:37769/130‑1 (MQ=255)
                                                                                                                                                   cTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACCCGCAGgat  <  2:253627/149‑1 (MQ=255)
                                                                                                                                                   cTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACCCGCAGgat  <  2:290767/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCC‑GAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACCCGCAGGAT  >  NC_002947/4710145‑4710440

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