Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,980,112 C→A 5.1% S120S (TCC→TCA) ‡ uspC → universal stress protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,980,1120CA 5.1% 261.8 / 2.5 78S120S (TCC→TCA) ‡uspCuniversal stress protein C
Reads supporting (aligned to +/- strand):  ref base C (29/45);  new base A (2/2);  total (31/47)
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.34e-01

AAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTA  >  NC_000913/1979963‑1980261
                                                                                                                                                      |                                                                                                                                                     
aaGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATcc                                                                                                                                                       >  2:1187113/1‑150 (MQ=255)
 aGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCt                                                                                                                                                      >  2:759336/1‑150 (MQ=255)
         cAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCg                                                                                                                                              <  2:1512525/150‑1 (MQ=255)
            gATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGaaa                                                                                                                                           <  2:1335161/150‑1 (MQ=255)
                 ggggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGt                                                                                                                                      <  2:1140845/150‑1 (MQ=255)
                  gggTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGtt                                                                                                                                     <  1:1230744/150‑1 (MQ=255)
                      aTCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTg                                                                                                                                 <  1:1167051/150‑1 (MQ=255)
                        cccGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGcc                                                                                                                               <  2:649632/150‑1 (MQ=255)
                          cGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTc                                                                                                                             >  1:1203681/1‑150 (MQ=255)
                           gTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCa                                                                                                                            <  2:1170411/150‑1 (MQ=255)
                             tGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAg                                                                                                                          <  1:1482055/150‑1 (MQ=255)
                               aCAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTg                                                                                                                        <  2:323225/150‑1 (MQ=255)
                                cAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                                       <  2:811334/150‑1 (MQ=255)
                                cAAGACGTTTATTGCCTACGGCGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGa                                                                                                                       <  1:1378634/150‑1 (MQ=255)
                                 aaGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAg                                                                                                                      <  1:1497764/150‑1 (MQ=255)
                                  aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                     >  1:920776/1‑150 (MQ=255)
                                  aGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAgg                                                                                                                     <  1:1177933/150‑1 (MQ=255)
                                        ttATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtg                                                                                                               <  1:405901/150‑1 (MQ=255)
                                           ttGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGc                                                                                                            <  2:355045/150‑1 (MQ=255)
                                             gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                          >  1:775612/1‑150 (MQ=255)
                                             gCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTg                                                                                                          >  2:1100862/1‑150 (MQ=255)
                                               atACTGTGAGTTAAGCTAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTTATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCTAAAATAGTTATTGCCTCAAGTGAGGTCGATGTGCTGtt                                                                                                        <  1:540262/149‑1 (MQ=255)
                                                tACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAag                                                                                                                                         <  2:1466178/116‑1 (MQ=255)
                                                tACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAag                                                                                                                                         >  1:1466178/1‑116 (MQ=255)
                                                 aCGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAg                                                                                                      >  2:624593/1‑150 (MQ=255)
                                                   ggTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGtt                                                                                                    <  2:342813/150‑1 (MQ=255)
                                                      gAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCa                                                                                                 <  2:351753/150‑1 (MQ=255)
                                                      gAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCAGCGAAAAGAGTTATTGCCTCAGGTGAGGTCGATGTGCTGTTAGTTCCa                                                                                                 <  2:1400464/150‑1 (MQ=255)
                                                       aGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCAc                                                                                                <  1:812965/150‑1 (MQ=255)
                                                       aGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAAGCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCAc                                                                                                <  2:949002/150‑1 (MQ=255)
                                                           aaGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCTCAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTAc                                                                                            <  1:422553/150‑1 (MQ=255)
                                                           aaGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTAc                                                                                            <  1:665846/150‑1 (MQ=255)
                                                                     atTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTaa                                                                                  >  2:1044092/1‑150 (MQ=255)
                                                                      ttCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAc                                                                                 <  2:1535259/150‑1 (MQ=255)
                                                                       tCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAAcc                                                                                <  2:1026806/150‑1 (MQ=255)
                                                                           gAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTc                                                                            <  1:759336/150‑1 (MQ=255)
                                                                              gTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTa                                                                         <  1:1065501/150‑1 (MQ=255)
                                                                                aTGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACg                                                                       >  1:1246626/1‑150 (MQ=255)
                                                                                aTGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACg                                                                       <  1:1193933/150‑1 (MQ=255)
                                                                                  gTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCa                                                                     <  2:775612/150‑1 (MQ=255)
                                                                                   tCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCaa                                                                    <  1:1111577/150‑1 (MQ=255)
                                                                                   tCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCaa                                                                    >  1:1305459/1‑150 (MQ=255)
                                                                                        aGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtg                                                                                                               >  2:533096/1‑102 (MQ=255)
                                                                                        aGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGAtg                                                                                                               <  1:533096/102‑1 (MQ=255)
                                                                                            tcatTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaa                                                           <  1:353583/150‑1 (MQ=255)
                                                                                             catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaaa                                                          >  2:1104611/1‑150 (MQ=255)
                                                                                             catTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGaaa                                                          >  1:762863/1‑150 (MQ=255)
                                                                                              atTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAg                                                         <  1:890232/150‑1 (MQ=255)
                                                                                              atTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAg                                                         >  2:1441279/1‑150 (MQ=255)
                                                                                                ttCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                                       <  2:1059666/150‑1 (MQ=255)
                                                                                                 tCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTa                                                      <  1:122186/150‑1 (MQ=255)
                                                                                                  cGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAg                                                     >  1:1106850/1‑150 (MQ=255)
                                                                                                   gATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGc                                                    <  2:1203681/150‑1 (MQ=255)
                                                                                                    aTTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCa                                                   >  2:1329037/1‑150 (MQ=255)
                                                                                                          tGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTa                                             >  2:793140/1‑150 (MQ=255)
                                                                                                            aTTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATc                                           >  2:1523358/1‑150 (MQ=255)
                                                                                                             tttGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATcg                                          <  2:617346/150‑1 (MQ=255)
                                                                                                              ttGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATcgc                                         >  1:1271361/1‑150 (MQ=255)
                                                                                                                 cGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgct                                      >  2:464177/1‑150 (MQ=255)
                                                                                                                  ggTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTACCCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctg                                     >  2:96705/1‑150 (MQ=255)
                                                                                                                     aatcacaatcacAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGctgctg                                  <  1:1370207/150‑1 (MQ=255)
                                                                                                                        cacaatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                               >  1:340446/1‑150 (MQ=255)
                                                                                                                        cacaatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCt                               >  2:789885/1‑150 (MQ=255)
                                                                                                                          caatcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCttt                             >  2:839777/1‑150 (MQ=255)
                                                                                                                            atcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATTTTCTTTTAGTTCCAGTTAAGGTAGATTAACCTCTCCTAATCAATCTTTTTAAATTTAGGAAGTTTATCGAGATGATGGCTTTcg                           <  2:814606/150‑2 (MQ=255)
                                                                                                                            atcacaGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTcc                           >  1:548479/1‑150 (MQ=255)
                                                                                                                                   actCTTTTCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTc                    >  1:1146419/3‑150 (MQ=255)
                                                                                                                                   actCTTTTCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTc                    <  2:456887/148‑1 (MQ=255)
                                                                                                                                     tCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTcgc                  >  1:32753/1‑150 (MQ=255)
                                                                                                                                      cTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgct                                   >  1:1156543/1‑132 (MQ=255)
                                                                                                                                      cTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCgctgct                                   <  2:1156543/132‑1 (MQ=255)
                                                                                                                                       ttttCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                                       >  1:499089/1‑111 (MQ=38)
                                                                                                                                       ttttCGCCGAGCAGGATGCTCGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGt                                                       <  2:499089/111‑1 (MQ=38)
                                                                                                                                          tCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAAc             <  1:1329037/150‑1 (MQ=255)
                                                                                                                                           cgcgAGCATCCCGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaag            <  2:1097816/150‑3 (MQ=255)
                                                                                                                                              cgAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTAt          >  1:1163835/1‑150 (MQ=255)
                                                                                                                                               gggCATACTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTACACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGTAAGGTAGCAACATTATCGCGCTTCTTGCTTTCCGCGCATCGCGGAACTAta         >  1:1362163/3‑149 (MQ=255)
                                                                                                                                                 gCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGAAGCGCGGaaagccagc       <  1:793140/150‑8 (MQ=255)
                                                                                                                                                   aTCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCt     <  1:1405683/150‑1 (MQ=255)
                                                                                                                                                     gcTGCTCGGCGAAAAGAGTTATTGCCTCTTGTGTGGTGGATGTGCTGTTAGTTGCACTTAGGGGAGATTAACCGCTCCTACGCAATGTTTGGAAAGGTAGCAACTTTATCGCGCTGGTGGGTTTCGGCGCTTCGCGGAACTATCTGCttt   <  2:880922/149‑1 (MQ=255)
                                                                                                                                                     ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt   >  2:1318159/1‑150 (MQ=255)
                                                                                                                                                     ccTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCATACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCttt   >  2:1323406/1‑150 (MQ=255)
                                                                                                                                                      cTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTa  <  2:762863/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGTTAATTCAGGATGCGGGGTATCCCGTTGACAAGACGTTTATTGCCTACGGTGAGTTAAGCGAACATATTCTGGAAGTATGTCACAAGCATCATTTCGATTTGGTGATTTGCGGTAATCACAATCACAGTTTCTTTTCGCGAGCATCCTGCTCGGCGAAAAGAGTTATTGCCTCAAGTGAGGTCGATGTGCTGTTAGTTCCACTTACGGGAGATTAACCGCTCCTACGCAAGCTTTGGAAAGGTAGCAACTTTATCGCGCTGCTGGCTTTCCGCGCTTCGCGGAACTATCTGCTTTA  >  NC_000913/1979963‑1980261

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