Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,993,134 G→C 5.6% P525A (CCC→GCC)  uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,993,1340GC 5.6% 246.7 / 5.4 72P525A (CCC→GCC) uvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base C (2/2);  total (34/38)
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.23e-01

CAACAACATTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAGTTCGGCGAAGACATTTTT  >  NC_000913/1992988‑1993281
                                                                                                                                                  |                                                                                                                                                   
caacaaCATTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCggg                                                                                                                                                     <  2:17163/147‑1 (MQ=255)
      caTTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCt                                                                                                                                              >  3:35048/1‑148 (MQ=255)
      caTTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCt                                                                                                                                              >  6:233109/1‑148 (MQ=255)
      caTTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCt                                                                                                                                              >  8:21681/1‑148 (MQ=255)
        tttGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTgg                                                                                                                                            <  5:187952/148‑1 (MQ=255)
                    ttGGCCCGACGCCTTCAATGGTATCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaa                                                                                                                                 <  1:5104/147‑1 (MQ=255)
                    ctGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTaa                                                                                                                                 <  6:105066/146‑1 (MQ=255)
                           gACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCa                                                                                                                         <  7:21681/148‑1 (MQ=255)
                           gACGCCTTCAATGGTTTCCAGGGAACTGGTATATTTGACCTTCGCCCGTTTTTTACGGTGCCAGCCAATCGCGTGATCATGTGATTCATCGCGAATATTCTCGATAACATGCAGCGCGGGTGAATCTCGCGGCAAACTAAATCCCTCa                                                                                                                         >  6:354798/1‑148 (MQ=255)
                                cTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTATTACGTTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGAAAACATGCAGCGCGGGTTAATCTGGCGGCAAACTAAATCCCTCACCTTc                                                                                                                    <  3:40833/148‑1 (MQ=255)
                                   cAATTGTTTCCAGGGAACTGGTATATTTGACCTTCGCCCGTTTTTTACTTTGCACGCCAATCGAGTGATCATGTGATTAATCGCGAATATGCTGTATAAAATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCg                                                                                                                  >  1:302460/1‑147 (MQ=255)
                                       ggTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTc                                                                                                              >  7:319261/1‑147 (MQ=255)
                                             cttggAACTGGTATTTATGACCTTCGCCGTTTTTTTAGGGTGCCCGCCTATCGGGTGATCCTGTGATTCTTCGCGAATATGCTGGATAACATGCAGCGCGGGTTAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAAc                                                                                                       <  6:277241/145‑1 (MQ=255)
                                              aGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCa                                                                                                                                       <  1:297574/115‑1 (MQ=255)
                                              aGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCa                                                                                                                                       >  2:297574/1‑115 (MQ=255)
                                              aGGGAACTGGTAGTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGaa                                                                                                        >  8:87354/1‑146 (MQ=255)
                                               gggAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAg                                                                                                     >  3:24170/1‑148 (MQ=255)
                                               gggAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGATGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAg                                                                                                     >  2:58838/1‑148 (MQ=255)
                                                ggAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGc                                                                                                    >  2:211059/1‑148 (MQ=255)
                                                  aaCTGGTATTTTTGACCTTCGCCCGTTTTGTACGGTGCCCGCCAATCGCGTGATCATGTGATGCCGCCCGAATATGCTGGATAACATGCCGCGGGGGGGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCGCAAAGAACAGCGt                                                                                                  <  2:248743/148‑1 (MQ=255)
                                                    cTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGttt                                                                                                <  8:377115/148‑1 (MQ=255)
                                                    cTGGTATTTTTGACCTTCGCCCGTGTTTTACGGTGCCCGCCTATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGttt                                                                                                <  2:374546/148‑1 (MQ=255)
                                                      ggTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                                                                                               >  3:23800/1‑147 (MQ=255)
                                                       gTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcc                                                                                              <  1:58838/147‑1 (MQ=255)
                                                       gTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtc                                                                                               >  1:135429/1‑146 (MQ=255)
                                                                   ttCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGAGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCCAAGAACAGCGTTTCCAGTCCCGCCTTAc                                                                                 >  6:201448/1‑148 (MQ=255)
                                                                    tCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc                                                                                 >  7:202794/1‑147 (MQ=255)
                                                                    tCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc                                                                                 >  7:364411/1‑147 (MQ=255)
                                                                    tCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc                                                                                 >  5:175547/1‑147 (MQ=255)
                                                                    tCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTAc                                                                                 >  5:268682/1‑147 (MQ=255)
                                                                     cGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGa                                                                               >  7:161698/1‑148 (MQ=255)
                                                                      gCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGAt                                                                              <  4:23800/148‑1 (MQ=255)
                                                                      gCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGAt                                                                              <  1:211059/148‑1 (MQ=255)
                                                                        ccGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCt                                                                            <  4:35048/148‑1 (MQ=255)
                                                                           ttttttACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCt                                                                         >  6:301892/1‑148 (MQ=255)
                                                                               ttACGGTGCCCGCCAATCGCGTGATCAGGTGATTCATCGCGAATATGCTGGATACCATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCt                                                                      <  2:365734/147‑1 (MQ=255)
                                                                                    gTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGGCCAGCCTTACGATCTGGTCCTTTGGCa                                                                <  1:179579/148‑1 (MQ=255)
                                                                                       cccGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAAGCGGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACg                                                             <  3:71676/148‑1 (MQ=255)
                                                                                          gCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCa                                                          <  3:95542/148‑1 (MQ=255)
                                                                                          gCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCa                                                          <  4:97288/148‑1 (MQ=255)
                                                                                            cAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAg                                                        <  2:287917/148‑1 (MQ=255)
                                                                                             aaTCGCGTGATCATGTGATTAATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTACGGCTCAAAGAACAGCGATTCCAGTACCGCCTTACGCTCTGCTCCTTTGGCAACGCCAAg                                                        >  1:337940/1‑147 (MQ=255)
                                                                                              aTCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAg                                                        >  5:145276/1‑146 (MQ=255)
                                                                                                cgcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTag                                                     >  4:302259/1‑147 (MQ=255)
                                                                                                 gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTagca                                                   <  5:301892/148‑1 (MQ=255)
                                                                                                 gcgTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTagc                                                    <  6:175547/147‑1 (MQ=255)
                                                                                                   gTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTagcagc                                                 >  6:154128/1‑148 (MQ=255)
                                                                                                      aTCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGga                                              >  4:86912/1‑148 (MQ=255)
                                                                                                        cATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGCTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcc                                                                                              <  8:329224/98‑1 (MQ=255)
                                                                                                        cATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGCTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcc                                                                                              <  8:334352/98‑1 (MQ=255)
                                                                                                        cATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGCTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcc                                                                                              >  7:334352/1‑98 (MQ=255)
                                                                                                        cATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGCTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTtcc                                                                                              >  7:329224/1‑98 (MQ=255)
                                                                                                           gtgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAtt                                         >  8:344689/1‑148 (MQ=255)
                                                                                                           gtgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAtt                                         >  8:344158/1‑148 (MQ=255)
                                                                                                            tgATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAttt                                        <  8:161698/148‑1 (MQ=255)
                                                                                                              aTTCATCGCGAATATGCTGGATAACATGCAGCGCGGGGGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAtttt                                       <  3:86912/147‑1 (MQ=255)
                                                                                                                tCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTAt                                    >  5:321500/1‑148 (MQ=255)
                                                                                                                   tCGCTAATATGCTGGATATCATGCAGCGCGTGTTAATCTGGCGTCAAAATAAATCCCTCTCCTTGAGGCTCACAGAACATCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGTCAACGCCAAGTAGCATCTTATGATTTGTATcac                                 <  2:337940/148‑3 (MQ=255)
                                                                                                                   tCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATccc                                 <  8:202794/148‑1 (MQ=255)
                                                                                                                    cgcgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAttt                                        <  7:372637/140‑1 (MQ=255)
                                                                                                                    cgcgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGAttt                                        >  8:372637/1‑140 (MQ=255)
                                                                                                                    cgcgAATATGCTGGATAACATGCAGAGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGgatgat                                          >  6:157246/1‑138 (MQ=255)
                                                                                                                    cgcgAATATGCTGGATAACATGCAGAGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGgatgat                                          <  5:157246/138‑1 (MQ=255)
                                                                                                                      cgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATg                              <  3:351810/148‑1 (MQ=255)
                                                                                                                      cgAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTACAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATg                              <  8:319261/148‑1 (MQ=255)
                                                                                                                             gCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCaa                                                         >  7:176632/1‑114 (MQ=255)
                                                                                                                             gCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCaa                                                         <  8:176632/114‑1 (MQ=255)
                                                                                                                              cTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCATCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATcc                      <  5:233109/148‑1 (MQ=255)
                                                                                                                                 gATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAg                    <  4:24170/147‑1 (MQ=255)
                                                                                                                                   tAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAGTTc                 <  7:27927/148‑1 (MQ=255)
                                                                                                                                        tGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAGTTCGGCGa            <  2:135429/148‑1 (MQ=255)
                                                                                                                                              gcgGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAg                                                        >  6:218707/1‑98 (MQ=255)
                                                                                                                                              gcgGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAg                                                        <  5:218707/98‑1 (MQ=255)
                                                                                                                                                  gTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAGTTCGGCGAAGACAttttt  <  7:87354/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAACAACATTTGCCGACGTTTTGGCCCGACGCCTTCAATGGTTTCCAGGGAACTGGTATTTTTGACCTTCGCCCGTTTTTTACGGTGCCCGCCAATCGCGTGATCATGTGATTCATCGCGAATATGCTGGATAACATGCAGCGCGGGTGAATCTGGCGGCAAACTAAATCCCTCACCTTCCGGCTCAAAGAACAGCGTTTCCAGTCCAGCCTTACGATCTGCTCCTTTGGCAACGCCAAGTAGCAGCGGATGATTTTTATCCCATGAGACATCCAGTTCGGCGAAGACATTTTT  >  NC_000913/1992988‑1993281

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