Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,424,895 C→G 5.7% intergenic (+483/‑29) PP_3039 → / → PP_3040 hypothetical protein/holin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,424,8950CG 5.7% 221.5 / 3.6 70intergenic (+483/‑29)PP_3039/PP_3040hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base G (2/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

AGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCG  >  NC_002947/3424749‑3425029
                                                                                                                                                  |                                                                                                                                      
aGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGgcgc                                                                                                                                      >  6:23849/1‑149 (MQ=255)
   cATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgg                                                                                                                                   <  4:102381/149‑1 (MQ=255)
    aTTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCggg                                                                                                                                  >  7:221599/1‑149 (MQ=255)
     ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgggg                                                                                                                                 >  3:165097/1‑149 (MQ=255)
          cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                            <  6:134428/149‑1 (MQ=255)
          cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                            <  8:221599/149‑1 (MQ=255)
            gaggTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGAGCTGTAATAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTTAAAACCCCGGCCCTTGCGCCGGGGTTTTTGc                                                                                                                          <  8:138855/146‑1 (MQ=255)
             aCGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                           <  5:220315/147‑3 (MQ=255)
                tttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCttt                                                                                                                      <  3:336524/149‑1 (MQ=255)
                   aCCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                   >  4:166371/1‑149 (MQ=255)
                   aCCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                    <  2:156548/148‑1 (MQ=255)
                    ccACTTATTGCCGTGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                   <  3:223685/148‑1 (MQ=255)
                    ccACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCATCGAGACGTGCGTGCGCACAAAGCGTCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                   <  2:85545/148‑1 (MQ=255)
                     cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggggg                                                                                                                 >  5:286312/1‑149 (MQ=255)
                     cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCGGGCGCCGGGGGTTTTGGCTTTggggg                                                                                                                 >  8:294894/1‑149 (MQ=255)
                      aCTTATTGCCGGGTGGCATTAAATCGCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCGATCTTCGAGACGGGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGa                                                                                                                <  5:53381/149‑1 (MQ=255)
                       cTTATTGCCGGGGGGCATTAACACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAACTACCCATCTTCGCGACGGGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAACCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                               <  7:87457/149‑1 (MQ=255)
                              ccGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATg                                                                                                         >  6:168975/1‑148 (MQ=255)
                              ccGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATg                                                                                                         <  3:166371/148‑1 (MQ=255)
                              ccGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGa                                                                                                        <  4:165097/149‑1 (MQ=255)
                                   tGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCATCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGc                                                                                                   <  7:340657/149‑1 (MQ=255)
                                     gCATTAAAACTCGCTTGCTGCCACCGTAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGa                                                                                                 <  7:231938/149‑1 (MQ=255)
                                       aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGag                                                                                                >  5:317349/1‑148 (MQ=255)
                                       aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGag                                                                                                >  2:1678/1‑148 (MQ=255)
                                        ttAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttc                                                                                                                            >  6:241552/1‑118 (MQ=39)
                                        ttAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttc                                                                                                                            <  5:241552/119‑2 (MQ=39)
                                         tAAAACTCGCTTGCTGCCACCGGAAGCGCCCTGTAAAAAGTCCCCACCTTCGAGACGGGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagcag                                                                                             <  6:360838/149‑1 (MQ=255)
                                           aaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAgcagca                                                                                           >  5:310509/1‑149 (MQ=255)
                                                 gCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGTGCGCAAAAAACGGCCCGCCCAACACTTAAAAACCCCGCCCGGTCGCCGGGTTTTTTTCCGTTGGGGGGTTCGGTGAACACCGCGAAGCAAgaggt                                                                                     >  6:228748/1‑145 (MQ=255)
                                                 gCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCCCAAACCGGCCCGCCAAACACAGAAAACCCCGGCCCAGGCGCCGGGGGTTTTTGCAGTGGGGGATTAGATGAACAGCGCGCAGCAACCgg                                                                                      >  6:233946/1‑147 (MQ=255)
                                                  ctcgctgcCACCGGAATCGACCTGAAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTa                                                                                    <  6:83401/146‑1 (MQ=255)
                                                     gctgcCACCGGAATCGACCAGTAAAAAGTACCCATCTTCGAGACGTGAGCGCGCACAAAGCGCCCCCCCAAACACTGAAAACCCACGGCCTGGCTCCCGGGGTTTTGCACTTGCGGGATTCGAAGAACATCGGGCAGCAACCGGTAGtt                                                                                 >  6:320052/1‑149 (MQ=255)
                                                         ccACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                             <  7:294894/149‑1 (MQ=255)
                                                           aCCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATg                                                                           >  3:75279/1‑149 (MQ=255)
                                                                aaTCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAAt                                                                      >  2:308142/1‑149 (MQ=255)
                                                                      ccTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                <  5:168975/149‑1 (MQ=255)
                                                                       cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTTCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                               >  7:209303/1‑149 (MQ=255)
                                                                       cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCCGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                               >  7:220345/1‑149 (MQ=255)
                                                                        tGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                <  8:143288/147‑1 (MQ=255)
                                                                        tGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                >  7:143288/1‑147 (MQ=255)
                                                                            aaaaGTACGCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCGAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGTATTCGATGAACAGCGAGCAGCAAGCGTTAGGTGAGATGCCAATCTGGTTGGTGa                                                           <  5:320052/148‑1 (MQ=255)
                                                                                gTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                      <  4:75279/149‑1 (MQ=255)
                                                                                gTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                      <  4:102023/149‑1 (MQ=255)
                                                                                gTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                      <  6:286312/149‑1 (MQ=255)
                                                                                gTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                      <  8:209303/149‑1 (MQ=255)
                                                                                   cccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCt                                                    <  8:220345/148‑1 (MQ=255)
                                                                                    ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGt                                                  <  6:310509/149‑1 (MQ=255)
                                                                                    ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGt                                                  >  5:217598/1‑149 (MQ=255)
                                                                                      aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTcc                                                >  1:310543/1‑149 (MQ=255)
                                                                                       gcTTCGAGCCGGGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAACCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATCAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTccc                                               <  2:132460/148‑1 (MQ=255)
                                                                                            gagaCGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGc                                                                          <  3:263216/117‑1 (MQ=255)
                                                                                             agaCGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcg                                         >  2:243947/1‑149 (MQ=255)
                                                                                                   gCGGGCGAACAAAGCGGCCCGCCAAACACTTAAAACCCCGGCCCTGGCGCCGTGGTTTTTGCCTTTGGGGGTTTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCgggg                                   <  8:343398/149‑1 (MQ=255)
                                                                                                       gcgcACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgtgt                                >  2:133862/1‑148 (MQ=255)
                                                                                                             aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATc                                                                     <  8:63812/105‑1 (MQ=255)
                                                                                                             aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATc                                                                     >  7:63812/1‑105 (MQ=255)
                                                                                                             aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCga                         <  5:23849/149‑1 (MQ=255)
                                                                                                                 cGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCCTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCt                                                    >  1:214271/1‑118 (MQ=255)
                                                                                                                  ggCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgt                    <  1:1678/149‑1 (MQ=255)
                                                                                                                     ccGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGc                 >  4:127214/1‑149 (MQ=255)
                                                                                                                         cAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGc                                                                          <  6:213332/88‑1 (MQ=255)
                                                                                                                         cAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGc                                                                          >  5:213332/1‑88 (MQ=255)
                                                                                                                          aaaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGcc            <  6:217598/149‑1 (MQ=255)
                                                                                                                              acTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACa        <  1:133862/149‑1 (MQ=255)
                                                                                                                                tGAAAACCCCGGCCCTGGCGCCGTGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGAc         <  7:189848/146‑1 (MQ=255)
                                                                                                                                tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGAc         >  8:189848/1‑146 (MQ=255)
                                                                                                                                tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg      >  3:301643/1‑149 (MQ=255)
                                                                                                                                tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg      >  2:217182/1‑149 (MQ=255)
                                                                                                                                  aaaaCCTCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc    >  7:33783/1‑149 (MQ=255)
                                                                                                                                  aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc    >  5:36405/1‑149 (MQ=255)
                                                                                                                                    aaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCg  <  6:317349/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
AGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCG  >  NC_002947/3424749‑3425029

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