Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,700 T→A 5.3% intergenic (+37/+23) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,7000TA 5.3% 454.8 / 11.6 131intergenic (+37/+23)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base T (52/72);  new base A (2/5);  total (54/77)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 136 total aligned reads displayed.

AAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCG  >  NC_000913/681553‑681846
                                                                                                                                                   |                                                                                                                                                  
aaaaaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTa                                                                                                                                                   <  3:671865/149‑1 (MQ=255)
    atgatgAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                                                                               <  1:1241526/149‑1 (MQ=255)
    atgatgAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                                                                               <  3:648717/149‑1 (MQ=255)
      gatgaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccgc                                                                                                                                             <  4:871804/149‑1 (MQ=255)
          aGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                                                                         <  1:1031555/149‑1 (MQ=255)
            ggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                                                                       <  1:309226/149‑1 (MQ=255)
              tGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCa                                                                                                                                     <  2:244486/149‑1 (MQ=255)
               gAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAg                                                                                                                                    >  1:124245/1‑149 (MQ=255)
                aaaaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGc                                                                                                                                   <  3:617817/149‑1 (MQ=255)
                 aaaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGcc                                                                                                                                  <  1:50958/149‑1 (MQ=255)
                  aaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                                                                       <  4:1565238/143‑1 (MQ=255)
                  aaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                                                                       >  3:1565238/1‑143 (MQ=255)
                  aaaCCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCg                                                                                                                                 >  1:769505/1‑149 (MQ=255)
                     ccACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                                                                              <  2:715544/149‑1 (MQ=255)
                      cacaGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGtt                                                                                                                             >  2:364857/1‑149 (MQ=255)
                       acaGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                                                                            >  4:563432/1‑149 (MQ=255)
                        caGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                                                                           <  3:278823/149‑1 (MQ=255)
                         agagCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTAt                                                                                                                          <  2:193085/149‑1 (MQ=255)
                           agCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCgtgt                                                                                                                              <  1:1398942/143‑1 (MQ=255)
                           agCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATgg                                                                                                                        <  3:563432/149‑1 (MQ=255)
                           agCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATgg                                                                                                                        <  1:890967/149‑1 (MQ=255)
                              aTTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                                                                     <  2:1248209/149‑1 (MQ=255)
                              aTTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                                                                     <  1:858189/149‑1 (MQ=255)
                               ttACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtc                                                                                                                    <  2:1099690/149‑1 (MQ=255)
                                tACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtct                                                                                                                   >  2:709703/1‑149 (MQ=255)
                                 aCCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtctc                                                                                                                  >  3:1372514/1‑149 (MQ=255)
                                  ccAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                                                                                 >  4:1135720/1‑149 (MQ=255)
                                   cAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCaa                                                                                                                >  2:75945/1‑149 (MQ=255)
                                    aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                               <  3:1144718/149‑1 (MQ=255)
                                    aGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                               <  1:415584/149‑1 (MQ=255)
                                     gCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                              >  3:171236/1‑149 (MQ=255)
                                     gCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                                                                              <  2:682931/149‑1 (MQ=255)
                                      cATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATgg                                                                                                             >  1:234319/1‑149 (MQ=255)
                                      cATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATgg                                                                                                             <  2:887398/149‑1 (MQ=255)
                                           tCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaaa                                                                                                        >  4:125522/1‑149 (MQ=255)
                                           tCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTaaa                                                                                                        >  3:901113/1‑149 (MQ=255)
                                               aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                                    <  2:769505/149‑1 (MQ=255)
                                               aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                                    >  4:1274776/1‑149 (MQ=255)
                                                aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                                   >  1:920970/1‑149 (MQ=255)
                                                aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                                   <  4:560250/149‑1 (MQ=255)
                                                aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                                   <  2:124245/149‑1 (MQ=255)
                                                   aGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgct                                                                                                <  3:155919/149‑1 (MQ=255)
                                                     tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcg                                                                                              >  1:1074015/1‑149 (MQ=255)
                                                        aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGttt                                                                                                                                                    <  1:589346/92‑1 (MQ=255)
                                                        aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGttt                                                                                                                                                    <  1:497968/92‑1 (MQ=255)
                                                        aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGttt                                                                                                                                                    >  2:497968/1‑92 (MQ=255)
                                                                  gTGATAGGTTATGCTCTCGGCAAAATGGCGATGTTGTTTAGCTGACGATGTGTTTTTGGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                                                                                 <  4:378955/149‑1 (MQ=255)
                                                                       aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCAt                                                                            <  1:1329219/149‑1 (MQ=255)
                                                                       aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCAt                                                                            >  1:223895/1‑149 (MQ=255)
                                                                        ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                           <  1:359504/149‑1 (MQ=255)
                                                                         gTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                          >  1:1258383/1‑149 (MQ=255)
                                                                          ttATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGATGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCa                                                                         >  4:505930/1‑149 (MQ=255)
                                                                           tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                                                                        >  2:480280/1‑149 (MQ=255)
                                                                                tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                                    <  3:283547/116‑1 (MQ=255)
                                                                                tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                                    >  4:283547/1‑116 (MQ=255)
                                                                                tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                   >  3:556764/1‑149 (MQ=255)
                                                                                tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                   >  4:1547870/1‑149 (MQ=255)
                                                                                tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCAGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                   >  2:448500/1‑149 (MQ=255)
                                                                                     gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata                                                              >  2:1097539/1‑149 (MQ=255)
                                                                                          aTCGCGTGGTGGTTTAGCTGATGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCTTCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCTCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt                                                         <  4:1478265/149‑1 (MQ=255)
                                                                                          aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt                                                         <  2:150945/149‑1 (MQ=255)
                                                                                          aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt                                                         <  3:1135720/149‑1 (MQ=255)
                                                                                          aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTGCATCCACGCGGGCCATATCGtt                                                         <  1:968576/149‑1 (MQ=255)
                                                                                          aTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACCTCATCCACGCGGGCCATATCGtt                                                         >  1:700334/1‑149 (MQ=255)
                                                                                           tCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt                                                        <  2:1366409/149‑1 (MQ=255)
                                                                                            cgcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg                                                       <  1:861129/149‑1 (MQ=255)
                                                                                             gcgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGa                                                      <  4:699564/149‑1 (MQ=255)
                                                                                                 tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                                                                     >  2:789626/1‑82 (MQ=255)
                                                                                                 tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAt                                                                                                                     <  1:789626/82‑1 (MQ=255)
                                                                                                   ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa                                                                >  3:1245659/1‑133 (MQ=255)
                                                                                                   ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa                                                                <  4:1245659/133‑1 (MQ=255)
                                                                                                   ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCaa                                                <  4:516744/149‑1 (MQ=255)
                                                                                                      tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                                             >  4:1403735/1‑149 (MQ=255)
                                                                                                      tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                                             <  2:920970/149‑1 (MQ=255)
                                                                                                      tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                                             >  3:886716/1‑149 (MQ=255)
                                                                                                      tttAGCTGACGATGCGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                                             >  4:552142/1‑149 (MQ=255)
                                                                                                       ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgcc                                            <  1:728943/149‑1 (MQ=255)
                                                                                                         aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaat                                                                                                       >  1:1359888/1‑88 (MQ=255)
                                                                                                         aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaat                                                                                                       <  2:1359888/88‑1 (MQ=255)
                                                                                                            tGACGATGTGGTTTTTGCGCCGATGCCGGAGGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCCATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCTACGCGGGCCATATCGTTTGACTGCAACGCCGCATc                                       >  4:248539/1‑149 (MQ=255)
                                                                                                              aCGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCaa                                     >  2:1413561/1‑149 (MQ=255)
                                                                                                               cGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCAGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCaaa                                    <  1:448500/149‑1 (MQ=255)
                                                                                                                   gtgtTTGTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCg                                <  2:591020/149‑1 (MQ=255)
                                                                                                                    tgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCgg                               <  1:201759/149‑1 (MQ=255)
                                                                                                                      ttttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                             >  3:968379/1‑149 (MQ=255)
                                                                                                                      ttttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTg                             >  2:35971/1‑149 (MQ=255)
                                                                                                                         tttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATc                          <  1:946225/149‑1 (MQ=255)
                                                                                                                          ttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGACTAAATAATCGCTTGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCa                         <  4:1331138/149‑1 (MQ=255)
                                                                                                                           tGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCAc                        >  4:894150/1‑149 (MQ=255)
                                                                                                                            gcgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCAcc                       >  3:1442411/1‑149 (MQ=255)
                                                                                                                                cGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGc                   <  3:220558/149‑1 (MQ=255)
                                                                                                                                       ggATGCGGCGTAAACGCTGCATCCGGCATCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtctc            <  1:1555049/149‑1 (MQ=255)
                                                                                                                                         aTGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCtctt          <  3:1412696/149‑1 (MQ=255)
                                                                                                                                          tGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCtcttc         <  3:125522/149‑1 (MQ=255)
                                                                                                                                             aGCGTTTACGCCGCGTCCAGCAGCCGTGGTTATGGGATCTCAATGGCTAAATAATCGCTCTCCCGCCGCCGCACTTCATCCACGCGGGCCATATCGTTTGACTGTAACGCCGCCTCCTAATCGGTGATCACCCGGCCGATCTCTTctcg      >  2:1266252/1‑149 (MQ=255)
                                                                                                                                             aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt                              >  2:1058992/1‑125 (MQ=255)
                                                                                                                                             aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt                              <  1:1058992/125‑1 (MQ=255)
                                                                                                                                                gTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgctc   <  4:1188812/149‑1 (MQ=255)
                                                                                                                                                gTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgctc   <  4:104547/149‑1 (MQ=255)
                                                                                                                                                gTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTctcgctc   <  1:722061/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCG  >  NC_000913/681553‑681846

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