Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,710,287 A→C 5.3% Y65D (TAC→GAC) ‡ PP_4168 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,710,2870AC 5.3% 242.5 / 4.4 76Y65D (TAC→GAC) ‡PP_4168hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base C (2/2);  total (39/37)
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.64e-01

GGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACC  >  NC_002947/4710142‑4710432
                                                                                                                                                 |                                                                                                                                                  
ggTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACt                                                                                                                                                  <  1:168881/148‑1 (MQ=255)
ggTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACt                                                                                                                                                  <  6:35548/148‑1 (MQ=255)
   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)
                   gATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAg                                                                                                                             <  5:260781/149‑1 (MQ=255)
                                    aCGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCt                                                                                                            <  3:173786/149‑1 (MQ=255)
                                    aCGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCt                                                                                                            >  4:173786/1‑149 (MQ=255)
                                          tGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                                                                                                                                          >  5:121767/1‑114 (MQ=255)
                                          tGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGa                                                                                                                                          <  6:121767/114‑1 (MQ=255)
                                            cgcTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGgcgcgc                                                                                                    <  3:407992/149‑1 (MQ=255)
                                             gctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTg                                                                                                           <  2:120367/141‑1 (MQ=255)
                                             gctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTg                                                                                                           >  1:120367/1‑141 (MQ=255)
                                              ctgctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                  >  3:247141/1‑149 (MQ=255)
                                                 ctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtc                                                                                                                          >  7:101691/1‑122 (MQ=255)
                                                 ctTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtc                                                                                                                          <  8:101691/122‑1 (MQ=255)
                                                 ctTCACATTCCGGCGCAGGCGGAGCCGGCCAGTCTGCGAACGTCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGTCTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATac                                                                                               <  1:284662/149‑1 (MQ=255)
                                                   tCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTAc                                                                                                                                                   <  7:256605/96‑1 (MQ=255)
                                                   tCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTAc                                                                                                                                                   >  8:256605/1‑96 (MQ=255)
                                                   tCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGgcgcgc                                                                                                    >  1:161207/1‑142 (MQ=255)
                                                     aCCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATacac                                                                                             >  6:203803/1‑147 (MQ=255)
                                                       cTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGt                                                                                         >  3:309400/1‑149 (MQ=255)
                                                       cTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGt                                                                                         <  4:247141/149‑1 (MQ=255)
                                                        ttCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                        >  2:253624/1‑149 (MQ=255)
                                                           cGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCgg                                                                                      <  6:244058/148‑1 (MQ=255)
                                                            ggCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGGCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTa                                                                                    <  8:303670/149‑1 (MQ=255)
                                                             gcgcACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCgg                                                                                      >  3:347840/1‑146 (MQ=255)
                                                                cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGAt                                                                                 >  8:69814/1‑148 (MQ=255)
                                                                cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGAt                                                                                 >  6:75229/1‑148 (MQ=255)
                                                                cACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                <  1:253624/149‑1 (MQ=255)
                                                                 aCGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCa                                                                               >  8:148793/1‑149 (MQ=255)
                                                                    cgGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCt                                                                            <  5:75229/149‑1 (MQ=255)
                                                                       ggccggccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTa                                                                         <  4:309400/149‑1 (MQ=255)
                                                                        gccggccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTAc                                                                        >  5:128608/1‑149 (MQ=255)
                                                                             ccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                              >  2:326976/1‑75 (MQ=255)
                                                                             ccAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCt                                                                                                                                              <  1:326976/75‑1 (MQ=255)
                                                                                 tCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAg                                                               <  5:203803/149‑1 (MQ=255)
                                                                                   tGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                >  3:372602/1‑130 (MQ=255)
                                                                                   tGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                <  4:372602/130‑1 (MQ=255)
                                                                                    gCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGc                                                             >  8:128697/1‑148 (MQ=255)
                                                                                     cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                  >  5:371134/1‑110 (MQ=255)
                                                                                     cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCa                                                                                                  <  6:371134/110‑1 (MQ=255)
                                                                                     cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                >  4:305314/1‑128 (MQ=255)
                                                                                     cGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                <  3:305314/128‑1 (MQ=255)
                                                                                              aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCAccc                                                        >  5:22165/1‑143 (MQ=255)
                                                                                              aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCAccc                                                        <  6:22165/143‑1 (MQ=255)
                                                                                              aCGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAg                                                   >  3:77903/1‑148 (MQ=255)
                                                                                                   ttCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                        >  8:359446/1‑106 (MQ=255)
                                                                                                   ttCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTc                                                                                        <  7:359446/106‑1 (MQ=255)
                                                                                                        gTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGc                                          >  8:106940/1‑147 (MQ=255)
                                                                                                         tCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAg                                       >  2:222918/1‑149 (MQ=255)
                                                                                                            cTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                >  8:336496/1‑105 (MQ=255)
                                                                                                            cTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATc                                                                                <  7:336496/105‑1 (MQ=255)
                                                                                                               tgGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGTCTCGCGCTGGTCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCAACCCAACAGCACCTCTGCCAGATGc                                   >  5:113906/1‑147 (MQ=255)
                                                                                                               tgGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGc                                   >  5:358931/1‑147 (MQ=255)
                                                                                                                 gAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTGGCCGAGTTGctgt                                <  2:227283/148‑4 (MQ=255)
                                                                                                                  aaaCTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGac                              <  7:148793/149‑1 (MQ=255)
                                                                                                                   aaCTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGacc                             <  7:69814/149‑1 (MQ=255)
                                                                                                                        ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCAGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccac                           >  1:290767/1‑147 (MQ=255)
                                                                                                                        ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCaa                                                      <  4:188458/119‑1 (MQ=255)
                                                                                                                        ccGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCaa                                                      >  3:188458/1‑119 (MQ=255)
                                                                                                                                gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccacc                          >  1:167839/1‑139 (MQ=255)
                                                                                                                                gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGaccacc                          <  2:167839/139‑1 (MQ=255)
                                                                                                                                gCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCa                <  4:347840/149‑1 (MQ=255)
                                                                                                                                 ccGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAg               <  7:106940/149‑1 (MQ=255)
                                                                                                                                  cGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAAATCAGCTGTACCAGCACCAGGCCACCCAACAGCaa                                                >  5:350710/1‑115 (MQ=255)
                                                                                                                                  cGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAAATCAGCTGTACCAGCACCAGGCCACCCAACAGCaa                                                <  6:350710/115‑1 (MQ=255)
                                                                                                                                     tCTGAGCCAGTACCTGGCACAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCg                                            >  7:358767/1‑116 (MQ=255)
                                                                                                                                     tCTGAGCCAGTACCTGGCACAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCg                                            <  8:358767/116‑1 (MQ=255)
                                                                                                                                        gtgCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGtt                                     <  7:293119/120‑1 (MQ=255)
                                                                                                                                        gtgCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGtt                                     >  8:293119/1‑120 (MQ=255)
                                                                                                                                              ggTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCAcc  <  6:358931/149‑1 (MQ=255)
                                                                                                                                                tACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGc                                                  <  1:379728/99‑1 (MQ=255)
                                                                                                                                                tACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGc                                                  >  2:379728/1‑99 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGTTCCTCGTCACGCACCGGATGCGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGGCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCAGGCTGGCGCTGGGCGCGCCATACACCAGTCGGTAGATCAGCTGTACCAGCACCAGGCCACCCAACAGCAACTCGGCCAGTTGCCAGACCACCAGGAACACCAGCATCCACAGCACC  >  NC_002947/4710142‑4710432

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: