Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 233,832 A→T 10.8% intergenic (+20/+43) PP_0180 → / ← PP_0181 cytochrome c family protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947233,8320AT10.8% 263.0 / 22.7 94intergenic (+20/+43)PP_0180/PP_0181cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (43/40);  new base T (4/6);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-01

GACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCC  >  NC_002947/233690‑233981
                                                                                                                                              |                                                                                                                                                     
gACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGAc                                                                                                                                                >  1:487295/1‑150 (MQ=255)
   tggctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACcgc                                                                                                                                             <  2:1622699/150‑1 (MQ=255)
     gctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCgcgc                                                                                                                                           >  2:297400/1‑150 (MQ=255)
      ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  1:1688490/150‑1 (MQ=255)
      ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCCGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  2:1200836/150‑1 (MQ=255)
       tggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCaa                                                                                                                                         <  1:1657944/150‑1 (MQ=255)
         gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAgg                                                                                                                                       >  1:1066077/1‑150 (MQ=255)
             ttCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtctc                                                                                                                                   <  2:726057/150‑1 (MQ=255)
                aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGa                                                                                                                                >  1:267284/1‑150 (MQ=255)
                aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGa                                                                                                                                >  2:1639466/1‑150 (MQ=255)
                 tGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGAt                                                                                                                               >  1:1771286/1‑150 (MQ=255)
                      aTGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTgc                                                                                                                          <  1:540949/150‑1 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                     >  1:1498868/1‑134 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGtc                                                                                                                                     <  2:1498868/134‑1 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                     >  1:1683215/1‑150 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                     >  1:115477/1‑150 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                     >  1:193652/1‑150 (MQ=255)
                           tACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAg                                                                                                                     >  1:293827/1‑150 (MQ=255)
                            aCTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGTCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAACGTCTCAGATCTTGCGCAAGc                                                                                                                    >  1:278255/1‑150 (MQ=255)
                             cTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGcc                                                                                                                   >  1:1317266/1‑150 (MQ=255)
                                 cTGTCGGTGAAAGCGGTGGCCTTGCTGGCAATCGTGTTCCCTTACGGCCGCTCGTGGTTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                               <  2:96010/150‑1 (MQ=255)
                                 cTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                               <  2:1317266/150‑1 (MQ=255)
                                  tGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCagg                                                                                                                                        <  1:501415/124‑3 (MQ=255)
                                  tGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCagg                                                                                                                                        >  2:501415/1‑122 (MQ=255)
                                    tCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCag                                                                                                                                         <  2:447733/121‑2 (MQ=255)
                                    tCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCag                                                                                                                                         >  1:447733/1‑120 (MQ=255)
                                     cGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTc                                                                                                           <  1:297400/150‑1 (MQ=255)
                                      ggCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGCGCCTGCAGGGGGCCATACTGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCa                                                                                                          >  1:1887280/1‑150 (MQ=255)
                                         gcAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggc                                                                                                       >  2:928710/1‑150 (MQ=255)
                                          cAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggct                                                                                                      <  2:115477/150‑1 (MQ=255)
                                                     cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGGGCAGTCCCGCAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGgcgc                                                                                           >  1:336645/1‑150 (MQ=255)
                                                     cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          <  1:692267/103‑1 (MQ=255)
                                                     cTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                          >  2:692267/1‑103 (MQ=255)
                                                      ttGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                          >  1:587639/1‑150 (MQ=255)
                                                        gCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                        <  2:1270965/150‑1 (MQ=255)
                                                          tgcccATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACATGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                      <  1:1813515/147‑1 (MQ=255)
                                                           ggCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTgag                                                                                     <  1:246656/150‑1 (MQ=255)
                                                           ggCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTgag                                                                                     >  2:1851377/1‑150 (MQ=255)
                                                            gCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                   <  2:1094309/87‑1 (MQ=255)
                                                            gCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAgg                                                                                                                                                   >  1:1094309/1‑87 (MQ=255)
                                                                 cGTGTTCCTTTACGGCCACTCGTGTCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                               <  1:591623/150‑1 (MQ=255)
                                                                   tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCt                                                                             <  1:360908/150‑1 (MQ=255)
                                                                    gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                               >  1:1554690/1‑115 (MQ=255)
                                                                    gtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGAt                                                                                                               <  2:1554690/115‑1 (MQ=255)
                                                                      tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTa                                                                          <  1:1253940/150‑1 (MQ=255)
                                                                         tttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                          <  1:875947/131‑1 (MQ=255)
                                                                         tttACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                          >  2:875947/1‑131 (MQ=255)
                                                                               gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGATCGCGCAAGGGCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGa                                                                 <  2:648174/150‑1 (MQ=255)
                                                                               gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGa                                                                 >  1:388934/1‑150 (MQ=255)
                                                                                ccGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGac                                                                <  1:928710/150‑1 (MQ=255)
                                                                                    tcttgggtggcCGAAGAGCTCTGTGCTGTGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACtt                                                            <  2:110631/147‑1 (MQ=255)
                                                                                    tCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCATGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACtt                                                            >  1:151890/1‑150 (MQ=255)
                                                                                      gTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTc                                                          <  2:587639/150‑1 (MQ=255)
                                                                                        ggctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCaa                                                        <  2:1410687/150‑1 (MQ=255)
                                                                                         gctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAAc                                                       >  2:557386/1‑150 (MQ=255)
                                                                                           tggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAActg                                                     <  2:1683215/150‑3 (MQ=255)
                                                                                            ggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGc                                                    <  2:193652/150‑1 (MQ=255)
                                                                                            ggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGc                                                    <  2:1771286/150‑1 (MQ=255)
                                                                                               cGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTc                                                 >  1:398995/1‑150 (MQ=255)
                                                                                               cGAAAAGAGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCCTGAGAACGTTCTGTAGCAATTGGACACTTTCAACGGCGTc                                                 <  2:1086207/150‑1 (MQ=255)
                                                                                                 aaaaGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCAt                                               <  2:398995/150‑1 (MQ=255)
                                                                                                  aaaGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATg                                              >  1:1889739/1‑150 (MQ=255)
                                                                                                       gcAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATATTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCa                                         <  1:1851377/150‑1 (MQ=255)
                                                                                                           tgctgcGGCCAACTGACTGGAGCATGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCt                                     <  2:151890/150‑1 (MQ=255)
                                                                                                             ctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTgg                                   >  2:1716640/1‑150 (MQ=255)
                                                                                                              tgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGt                                  >  1:690024/1‑150 (MQ=255)
                                                                                                                 ggCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                               >  1:241376/1‑150 (MQ=255)
                                                                                                                 ggCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                               <  2:1098332/150‑1 (MQ=255)
                                                                                                                    cAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTa                            >  1:1644364/1‑150 (MQ=255)
                                                                                                                       ctgactgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                               <  1:590364/144‑1 (MQ=255)
                                                                                                                       ctgactgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAg                               >  2:590364/1‑144 (MQ=255)
                                                                                                                          actgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTc                                 >  1:620425/1‑139 (MQ=255)
                                                                                                                          actgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTc                                 <  2:620425/139‑1 (MQ=255)
                                                                                                                             ggAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATc                   >  1:912186/1‑150 (MQ=255)
                                                                                                                              gAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                   <  1:509344/101‑1 (MQ=255)
                                                                                                                              gAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTg                                                                   >  2:509344/1‑101 (MQ=255)
                                                                                                                                gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg                >  1:161334/1‑150 (MQ=255)
                                                                                                                                 cttGCGCGGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                               >  1:833150/3‑86 (MQ=38)
                                                                                                                                 cttGCGCGGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                               <  2:833150/84‑1 (MQ=38)
                                                                                                                                 ccTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGc               <  2:388934/150‑1 (MQ=255)
                                                                                                                                   tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCAt                                               >  1:921132/1‑116 (MQ=255)
                                                                                                                                   tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCAt                                               <  2:921132/116‑1 (MQ=255)
                                                                                                                                   tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg             >  1:827691/1‑150 (MQ=255)
                                                                                                                                   tGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTg             >  1:1516203/1‑150 (MQ=255)
                                                                                                                                      gcgGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATc          >  1:1462820/1‑150 (MQ=255)
                                                                                                                                       cgGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCg         <  2:1887280/150‑1 (MQ=255)
                                                                                                                                          tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATg                                                                                 <  1:1201867/75‑1 (MQ=255)
                                                                                                                                          tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATg                                                                                 >  2:1201867/1‑75 (MQ=255)
                                                                                                                                          tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGATATGTAGCAAATATAGACTAGCTACCGCGGGTTGACGTATTCTTGTCAGGTAGCCGCCATCCGGCTGCTCGGTc      >  2:804952/1‑150 (MQ=255)
                                                                                                                                            ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCt                                                                             <  2:318334/77‑1 (MQ=255)
                                                                                                                                            ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCt                                                                             >  1:318334/1‑77 (MQ=255)
                                                                                                                                            ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTgg                                                                  <  2:406042/88‑1 (MQ=255)
                                                                                                                                            ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTgg                                                                  >  1:406042/1‑88 (MQ=255)
                                                                                                                                            ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg    >  2:1409121/1‑150 (MQ=255)
                                                                                                                                             cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc   >  2:1741403/1‑150 (MQ=255)
                                                                                                                                             cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc   <  2:1644364/150‑1 (MQ=255)
                                                                                                                                             cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc   >  2:1578193/1‑150 (MQ=255)
                                                                                                                                              acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                               <  2:612907/73‑1 (MQ=255)
                                                                                                                                              acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                               >  1:612907/1‑73 (MQ=255)
                                                                                                                                              acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGcc  >  1:1394584/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGCC  >  NC_002947/233690‑233981

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