Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 6.6% intergenic (+131/+14) rraB → / ← yjgM ribonuclease E inhibitor protein B/putative acetyltransferase YjgM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0200AT 6.6% 156.9 / 2.2 59intergenic (+131/+14)rraB/yjgMribonuclease E inhibitor protein B/putative acetyltransferase YjgM
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

GCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAG  >  NC_000913/4478882‑4479162
                                                                                                                                          |                                                                                                                                              
gCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGGATGCTCGCGAGTTGTAcgcgcg                                                                                                                                            >  1:395609/1‑143 (MQ=255)
 ccACTAGTTAAGCGTAGCGCATCAGGCCCAGGTAGTGCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATACACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCgg                                                                                                                                           <  2:413151/143‑1 (MQ=255)
  cgcTAATCCAGCGTAGCGCATCAGGCCCAGGTTGTTCCTGCATCCTGGTGAATTTGTAGCCCGCATAAGGCATTCATGTCACCTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGTGTTCAACTCGCggg                                                                                                                                          <  2:877433/141‑1 (MQ=255)
  cACTTACTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCggg                                                                                                                                          <  2:1072807/143‑1 (MQ=255)
        ttAAGCGTAGCGTATCAGGCCTAGGTAGTGCCTGCAAGCTGGCGAATTTGTAGCCCGGATAAGGCATTCATGCCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                    <  1:223029/143‑1 (MQ=255)
        ttAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAAGAAAACGTACGTTGTTAGTGTTCGACTGCCCGCGGGCTGAACCCGCGGGCGtttt                                                                                                                                    >  1:192735/1‑143 (MQ=255)
        ttAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                    >  1:138471/1‑143 (MQ=255)
         tAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGCTAAGGCATCCCTGTCACATCCGGCATAAACAACACGTACGTTGTTAATGTTCGAATACCCGCGTGCTGAACAGGCGGCGAttttt                                                                                                                                   >  1:863534/1‑143 (MQ=255)
              gTAGCGCATCAGGCCCAGGTAGTTCCTGCAAGCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTaga                                                                                                                              <  2:866975/143‑1 (MQ=255)
                 gcgcATGACGCCCAGTGAGTTCGTGCCACCTGGGGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                                                           <  2:819919/143‑1 (MQ=255)
                 gcgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGACAAGGCAGTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGtt                                                                                                                           <  2:979017/143‑1 (MQ=255)
                    cATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                                                                                                                        >  1:720169/1‑143 (MQ=255)
                      tCAGGCCTAGGTCGTTCCTGCAACCTGGTGAATTTGGAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCgcgc                                                                                                                      <  2:313659/143‑1 (MQ=255)
                       cAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGGTGAACACGCGCGCAGTTTTTTAGAGGTCGCGCa                                                                                                                     >  2:6524/1‑143 (MQ=255)
                       cAGGCCCAGGTAGTTCCTGCAACCAGGTGGATTTTTAGCCCGGATAAGGGATTCAAACAACAGCCGGCATCATCACATCTCGCGTTGTTGATGTTCGAATGCCCGCGAGTTTAACACGCGGGCATTCTTGTAGAGCTCGCGCa                                                                                                                     >  2:1103744/1‑143 (MQ=255)
                         ggCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGGATTCATGTCACATCCGGCATAAACAAATCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAAc                                                                                                                   <  2:566915/143‑1 (MQ=255)
                           cccAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                                                                 >  1:760726/1‑143 (MQ=255)
                              aGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCAGATCCAGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACACCCg                                                                                                              >  2:152274/1‑143 (MQ=255)
                              aGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                                                                              >  1:894725/1‑143 (MQ=255)
                              aGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTCTAGAGTTCGCACAACATCCg                                                                                                              >  2:637639/1‑143 (MQ=255)
                                gTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                                                            <  1:219971/143‑1 (MQ=255)
                                gTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTGGGTAATGGTCGAATGCCCGCGAGTTGAGCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                                                            >  1:778035/1‑143 (MQ=255)
                                gTAGTTCCTGCAACCGGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCa                                                                                                            <  1:338920/143‑1 (MQ=255)
                                 tAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCAGcgc                                                                                                           <  2:196447/143‑3 (MQ=255)
                                    ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                                                                        >  1:538675/1‑143 (MQ=255)
                                        tacatcatGGTGAAGTTGTAGCCCGGATACGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTGAATTTCGCTCAACATCCGCACTTCACaa                                                                                                    <  2:403318/136‑1 (MQ=255)
                                        tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                                                                    >  1:974382/1‑143 (MQ=255)
                                            cccTGGTGGATTTGGAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGa                                                                                                <  2:448827/142‑1 (MQ=255)
                                                       ttGTAGCGTGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                     <  2:102839/143‑1 (MQ=255)
                                                        cgtaccgcGGATAGGCCAGTCATGTCCCATCCTGCATAAAAAATCCGTACGTTGTTAATGTTCCAATGCCCGCGAGTTTAACACGCGGGCATTGTTTTAGAGTCCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                                                                    <  2:107562/136‑1 (MQ=255)
                                                          tAGCCCGGATGAGGCAGCCGAGTCAGATCCGGCATACGCAAAACGTACGTTGTTAATTGTCGAATGCGCGCGAGTTGACCACGCGGGCATTTTTTTAGATTTTGCGCAACTTCTGAACTTCACTATCGTCATGGCCCGTGCAg                                                                                  <  2:880336/143‑1 (MQ=255)
                                                                ggATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGCTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCCTTTTTTTAGAGTTCTCGCAACATCAGCCCTGCGCAATCGACATAGCCCGTGCAGACaacg                                                                            >  2:182426/1‑141 (MQ=255)
                                                                ggAGAACGCATTCATGTCACATCCGGCATAACCAAAAGGAACGTACCAAAAATAGGTATGCCCACGAGTTCAGCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                                            <  1:474394/143‑1 (MQ=255)
                                                                 taTAAGGCAGTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgcg                                                                           <  1:166562/142‑1 (MQ=255)
                                                                       gCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                                     >  1:868068/1‑143 (MQ=255)
                                                                         aTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCg                                                                   >  2:497614/1‑143 (MQ=255)
                                                                          ttCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCAGCGTTTTACGACCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                  >  2:796262/1‑143 (MQ=255)
                                                                               gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGc                                                             <  1:840919/143‑1 (MQ=255)
                                                                                     gccGGCAGAAAGAAAACGGACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                       <  2:325180/142‑1 (MQ=255)
                                                                                     gccGGAAGAAACAAAACGTACGTGGTTAAGGTGCGAAGGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                       <  1:122873/142‑1 (MQ=255)
                                                                                      ccGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                      <  2:894725/143‑1 (MQ=255)
                                                                                        ggCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGccc                                                    >  1:433498/1‑143 (MQ=255)
                                                                                              aaCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAACCCAAATGc                                              >  2:31209/1‑143 (MQ=255)
                                                                                               aCAAAACGTACGTTGTTAATGTTCGAATGCCCGAGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                                             >  2:1107867/1‑143 (MQ=255)
                                                                                                cAAAACGTACGTTGTTAGTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                                            <  1:605302/143‑1 (MQ=255)
                                                                                                   aaCGTACGTTGTTAATGTTCGAATGCCAGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATa                                         <  2:240990/143‑1 (MQ=255)
                                                                                                     cctacgtTGTTATTGTTCGATTGTCCGCGAGTTGATCTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATaaa                                       <  1:199693/141‑1 (MQ=255)
                                                                                                           tcgttAATGTTCGCATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCaa                                 <  1:229729/141‑1 (MQ=255)
                                                                                                                      cgcaTGCCCGCGAGTTGAACCCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTACCttt                      <  2:460784/140‑1 (MQ=255)
                                                                                                                            cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaaaaaa                <  1:497614/143‑1 (MQ=255)
                                                                                                                             ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg               >  1:310394/1‑143 (MQ=255)
                                                                                                                             ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg               >  1:966761/1‑143 (MQ=255)
                                                                                                                             ccGCGAGTTTAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAACGCGCCATGGCCTCCTTTCAAAAAg               >  2:1055579/1‑143 (MQ=255)
                                                                                                                             ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCCGCCAAGCCCCTAGTCGATAAGCTCAAAGCCCAAATGCTCGTAAAGCGCAATGGCTTCCTCTaacggag               >  2:839395/1‑138 (MQ=255)
                                                                                                                               gcgAGTTGAACGCCCGGGCACTGTTGTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGGGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATAGCTTCCTTTAAAAAAGCg             <  2:192735/143‑1 (MQ=255)
                                                                                                                                cgAGTTGAACACGCGTGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCgg            <  2:760726/143‑1 (MQ=255)
                                                                                                                                cgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCGCAATGGAGATGGCCAGTGCGGCCAAGCGCATAGTCGGTAGGCTCAAAGCCCAAATGCTCATAAAGCGGAGTTGCCTCCGTTAAAAGAGCgg            >  2:835182/1‑143 (MQ=255)
                                                                                                                                    ttGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGGCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGt        >  1:1051588/1‑143 (MQ=255)
                                                                                                                                     tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTGACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTAAAAGCCCCAATGCTCATCAAGCGCAATGGCTTCCTTTAGAAAAGCgtgtgtt       >  2:325229/1‑137 (MQ=255)
                                                                                                                                      gAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCTACATGGCCCGTGCAGCCAAGCGCATATTCGATATTCTCAACGCGCAAATGCTCCTCAAGCGCAATTGCTTCCTTTAAAGCAGCGGTCGttt      >  2:724436/1‑143 (MQ=255)
                                                                                                                                        acacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGAATTCATTTAAAAAAGCGGACGattaa    >  1:1109122/1‑141 (MQ=255)
                                                                                                                                          acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACCTGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAACGGCTTCCTTTAAAAAAGCGGTCGATTcccg  >  2:919605/1‑141 (MQ=255)
                                                                                                                                          acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTGGTTTCCAg  >  1:377221/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
GCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAG  >  NC_000913/4478882‑4479162

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