Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,424,935 T→A 6.4% intergenic (+480/‑32) 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,9350TA 6.4% 278.8 / 10.1 94intergenic (+480/‑32)PP_3039/PP_3040hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base A (3/3);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.79e-01

TTCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCG  >  NC_002947/3424786‑3425081
                                                                                                                                                      |                                                                                                                                                  
ttCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCt                                                                                                                                                    <  2:494674/150‑1 (MQ=255)
         cATTCCCCGGACGTTTACCACTTCTTGCCGGTTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCggg                                                                                                                                           <  1:898478/150‑1 (MQ=255)
          aTTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgggg                                                                                                                                          >  1:1261328/1‑150 (MQ=255)
             ccccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttt                                                                                                                                       <  2:1261328/150‑1 (MQ=255)
               ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                     >  1:449059/1‑150 (MQ=255)
               ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                     <  1:1518892/150‑1 (MQ=255)
                     gTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcagtgtt                                                                                                                               >  2:521922/1‑143 (MQ=255)
                       ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                             <  1:589228/150‑1 (MQ=255)
                       ttACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                             <  2:440107/150‑1 (MQ=255)
                        tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                            <  2:449059/150‑1 (MQ=255)
                        tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                            <  2:1360348/150‑1 (MQ=255)
                          ccACTTATTGCCGGGTGGCGTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggggg                                                                                                                          <  1:359414/150‑1 (MQ=255)
                          ccACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggggg                                                                                                                          >  1:834684/1‑150 (MQ=255)
                            aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTGCCTTTGGGGGAtt                                                                                                                       >  2:948489/1‑150 (MQ=255)
                            aCTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGACAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                        >  1:609525/1‑150 (MQ=255)
                             cTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTACAAAGTACCGTGCTGCGAGACTTGCCGCCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAtt                                                                                                                       <  1:628834/150‑1 (MQ=255)
                                aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                    <  2:1884836/150‑1 (MQ=255)
                                aTTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                                    <  1:458744/150‑1 (MQ=255)
                                  tGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATg                                                                                                                  <  2:963792/150‑1 (MQ=255)
                                          ggCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTGCCTTTGGGGGATTCGATGAACAGCGag                                                                                                         <  1:948489/150‑1 (MQ=255)
                                            cATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagc                                                                                                        >  2:1067047/1‑150 (MQ=255)
                                             aTTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                       >  2:905952/1‑150 (MQ=255)
                                              ttAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagcag                                                                                                      <  1:675060/150‑1 (MQ=255)
                                               tAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGaa                                                                                                                >  1:1895493/1‑139 (MQ=255)
                                               tAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGaa                                                                                                                <  2:1895493/139‑1 (MQ=255)
                                                aaaaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAgcagca                                                                                                    >  1:1009877/1‑150 (MQ=255)
                                                  aaCTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGCATCATTTGCCTTTGGGGGATTCGATGAACAGCGagcagagat                                                                                                  >  1:19084/1‑146 (MQ=255)
                                                       gCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgg                                                                                                                                            >  2:1051618/1‑103 (MQ=255)
                                                       gCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgg                                                                                                                                            <  1:1051618/103‑1 (MQ=255)
                                                          tgctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGtt                                                                                          <  1:905952/150‑1 (MQ=255)
                                                            ctgcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTga                                                                                        <  1:343705/150‑1 (MQ=255)
                                                              gcCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                      >  2:878387/1‑150 (MQ=255)
                                                                aaCCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATg                                                                                    >  2:951431/2‑150 (MQ=255)
                                                                     gAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                        <  1:315002/109‑1 (MQ=255)
                                                                     gAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                        >  2:315002/1‑109 (MQ=255)
                                                                     gAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAAt                                                                               <  2:1363749/150‑1 (MQ=255)
                                                                         cGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtgg                                                                           <  1:951431/150‑1 (MQ=255)
                                                                           aCCTGTAAAGACTACCCATCTTCGAAACGGGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAACCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                         <  1:890419/150‑1 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                             <  1:261503/98‑1 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                             >  2:261503/1‑98 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGc                                                                                                            <  1:1484040/115‑1 (MQ=255)
                                                                           aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGc                                                                                                            >  2:1484040/1‑115 (MQ=255)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                                    >  2:1449692/1‑87 (MQ=37)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                                    <  1:1449692/89‑3 (MQ=37)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAtgc                                                                                                                      >  2:1303579/1‑101 (MQ=255)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTc                                                                                                                      <  1:1303579/103‑1 (MQ=255)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCaa                                                                                                   <  2:1503678/122‑1 (MQ=255)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCaa                                                                                                   >  1:1503678/1‑122 (MQ=255)
                                                                             cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttgg                                                                       >  1:1747946/1‑150 (MQ=255)
                                                                              tGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttggt                                                                      <  2:834684/150‑1 (MQ=255)
                                                                              tGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                                        <  2:1581907/148‑1 (MQ=255)
                                                                              tGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                                        >  1:1581907/1‑148 (MQ=255)
                                                                                     aGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                               <  1:545085/150‑1 (MQ=255)
                                                                                          ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                   >  2:837771/1‑93 (MQ=255)
                                                                                          ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                   >  2:837917/1‑93 (MQ=255)
                                                                                          ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                   <  1:837771/93‑1 (MQ=255)
                                                                                          ccATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                                   <  1:837917/93‑1 (MQ=255)
                                                                                           cATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTcc                                                         >  1:467860/1‑150 (MQ=255)
                                                                                            aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                        <  2:1282948/86‑1 (MQ=255)
                                                                                            aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                        >  1:1282948/1‑86 (MQ=255)
                                                                                              cTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTg                                                      <  1:419973/150‑1 (MQ=255)
                                                                                               ttCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTgg                                                     >  2:1662934/1‑150 (MQ=255)
                                                                                               ttCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTgg                                                     >  2:876066/1‑150 (MQ=255)
                                                                                                    gaCGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcggc                                                <  1:876066/150‑1 (MQ=255)
                                                                                                         gCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGgt                                           <  2:232285/150‑1 (MQ=255)
                                                                                                              cgcACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAg                                      <  1:1067047/150‑1 (MQ=255)
                                                                                                                cacaAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTTGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGc                                    <  1:1766061/150‑1 (MQ=255)
                                                                                                                 acaAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCg                                   >  2:1759356/1‑150 (MQ=255)
                                                                                                                 acaAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCg                                   >  2:986577/1‑150 (MQ=255)
                                                                                                                   aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTccc                                                        >  2:1690212/1‑127 (MQ=255)
                                                                                                                   aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCgag                                 <  1:659323/150‑1 (MQ=255)
                                                                                                                   aaaGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTGGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTccc                                                        <  1:1690212/127‑1 (MQ=255)
                                                                                                                        ggCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgtg                            <  2:1747946/150‑1 (MQ=255)
                                                                                                                            cGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCgg                        <  1:1662934/150‑1 (MQ=255)
                                                                                                                                aaaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCACGCGTTAGTTGAGATGCCAATCTGGTTGTTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCg                    >  1:11192/1‑150 (MQ=255)
                                                                                                                                aaaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCg                    >  1:1756632/1‑150 (MQ=255)
                                                                                                                                 aaCACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGa                   <  2:845897/150‑1 (MQ=255)
                                                                                                                                  acacTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGtt                                                                                          >  1:1074001/1‑78 (MQ=255)
                                                                                                                                  acacTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGtt                                                                                          <  2:1074001/78‑1 (MQ=255)
                                                                                                                                      tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg              >  2:933030/1‑150 (MQ=255)
                                                                                                                                      tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg              <  1:1759356/150‑1 (MQ=255)
                                                                                                                                      tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTAAGATGCCAATCTGGTTGGTGGGCTTCGTGTCCCTGGTCGGCGGGGGGTCAGGCGAGATGTGGCGGGCGGGCATgg              >  2:1911451/1‑150 (MQ=255)
                                                                                                                                      caaaaaCCCCGGCGCCAGGGCCGGGGTTTTCGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg              <  2:945117/148‑1 (MQ=255)
                                                                                                                                       aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGc                      >  2:594769/2‑141 (MQ=255)
                                                                                                                                       aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGc                      <  1:594769/140‑1 (MQ=255)
                                                                                                                                       aaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGAGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGc             >  2:21293/2‑150 (MQ=255)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGAt                                                                                     <  1:1775973/77‑1 (MQ=255)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGAt                                                                                     >  2:1775973/1‑77 (MQ=255)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGcc            >  1:1137772/1‑150 (MQ=255)
                                                                                                                                        aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGcc            >  2:30987/1‑150 (MQ=255)
                                                                                                                                         aaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCg           >  1:1265268/1‑150 (MQ=255)
                                                                                                                                           cccccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc         <  1:1911451/149‑1 (MQ=255)
                                                                                                                                                 gCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGc   <  1:986577/150‑1 (MQ=255)
                                                                                                                                                    ccTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCg  >  2:83183/1‑149 (MQ=255)
                                                                                                                                                    ccTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCg  <  1:83183/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TTCAACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCCACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGTTAGCG  >  NC_002947/3424786‑3425081

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: