Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,020,927 G→C 6.7% S13R (AGC→AGG)  dppC ← dipeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,020,9270GC 6.7% 183.0 / 3.9 60S13R (AGC→AGG) dppCdipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base C (2/2);  total (29/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 = 2.84e-01

GTCATGCGGGGCGACCCAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACAGGATGC  >  NC_002947/1020780‑1021070
                                                                                                                                                    |                                                                                                                                                
gTCATGCGGGGCGACCCAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGAAGGCTCTGGTCCACCGGGc                                                                                                                                                 >  8:150124/1‑149 (MQ=255)
        gggCGACCCAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGt                                                                                                                                         <  2:17705/149‑1 (MQ=255)
          gCGACCCAAGGCGCGAACAGCGCGCAGAAC‑GCCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGCCAGCCCGGTgg                                                                                                                                       <  3:185632/149‑2 (MQ=255)
              cccAAGGCGCGAACAGCGCGCAGAACGACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGacac                                                                                                                                    >  1:2886/1‑149 (MQ=255)
               ccAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACg                                                                                                                                   >  1:55656/1‑148 (MQ=255)
                cAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGa                                                                                                                                 <  3:23369/149‑1 (MQ=255)
                cAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGa                                                                                                                                 >  5:308201/1‑149 (MQ=255)
                 aaGGCGCGAACACCCCGCAGAAC‑CCCACCCCGCCCCTGAACGCCCCGCCCACGCCCGCGCCCTTGGTGCGCGCGAAGGCCTGCCACCATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGAt                                                                                                                                <  4:165418/149‑1 (MQ=255)
                     cgcgAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTg                                                                                                                            <  7:225919/149‑1 (MQ=255)
                          aCAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATc                                                                                                                       <  3:251315/149‑1 (MQ=255)
                            aGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCgg                                                                                                                     <  2:55656/149‑1 (MQ=255)
                                cgcAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGt                                                                                                                                         <  3:282107/125‑1 (MQ=255)
                                cgcAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGt                                                                                                                                         >  4:282107/1‑125 (MQ=255)
                                      ac‑acCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCAt                                                                                                           <  1:225189/148‑1 (MQ=255)
                                       c‑accacCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATgg                                                                                                         >  6:211518/1‑149 (MQ=255)
                                           caccaGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGc                                                                                                        >  3:44926/1‑147 (MQ=255)
                                            accaGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAg                                                                                                     >  6:171341/1‑149 (MQ=255)
                                                   acaTGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCt                                                                                                >  5:139552/1‑147 (MQ=255)
                                                       gAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAgcgc                                                                                          >  1:276525/1‑149 (MQ=255)
                                                         aggccaggccCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCt                                                                                         >  3:42563/1‑148 (MQ=255)
                                                           gccaggccCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCtgat                                                                                      >  2:266220/1‑149 (MQ=255)
                                                              aggccCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCggg                                                                                                                    <  4:182606/117‑1 (MQ=255)
                                                              aggccCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCggg                                                                                                                    >  3:182606/1‑117 (MQ=255)
                                                                  ccATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATgg                                                                                                         >  4:82165/1‑123 (MQ=255)
                                                                  ccATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATgg                                                                                                         <  3:82165/123‑1 (MQ=255)
                                                                    aTGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGc                                                                             >  7:216381/1‑149 (MQ=255)
                                                                         cgcgCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGA‑GGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGt                                                                       >  5:145905/1‑149 (MQ=255)
                                                                            gcCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTg                                                                     <  7:150124/149‑1 (MQ=255)
                                                                             cccTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTgg                                                                    <  8:216381/149‑1 (MQ=255)
                                                                              ccTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGc                                                                   <  5:1171/149‑1 (MQ=255)
                                                                                  gttgCGCGCGAAGGCGTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAg                                                                <  6:308201/148‑1 (MQ=255)
                                                                                   ttgCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAgg                                                               <  4:42563/148‑1 (MQ=255)
                                                                                   ttgCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAgg                                                               <  2:2886/148‑1 (MQ=255)
                                                                                   ttgCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGc                                                              >  3:206139/1‑149 (MQ=255)
                                                                                    tgCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GGCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGcc                                                             <  7:206084/149‑1 (MQ=255)
                                                                                     gcgcgcGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGcc                                                             <  4:206139/148‑1 (MQ=255)
                                                                                       gcgcgAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCg                                                            >  6:59083/1‑147 (MQ=255)
                                                                                         gcgAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTaga                                                        <  5:171341/149‑1 (MQ=255)
                                                                                            aaGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGa                                                     <  5:59083/149‑1 (MQ=255)
                                                                                                 ttGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCAACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATTCGTGGGTTGGCCAGGCCGTAGAGGATGTc                                                 >  6:137727/1‑148 (MQ=255)
                                                                                                         aTTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGaa                                        <  5:99852/149‑1 (MQ=255)
                                                                                                               tGTACGGGGATGGGTATAGCTGGCTGTGGTCCACATGGCTGGCCGGTGACACTGATTTTGGAATCGGGCGA‑GTGATGGCAAGGTCCTCAGCGCTGATGTCGGTTGGGTGGGTTGGCATGGGCGTAGAGGATGTCCACGACGAAGTTGa                                   <  4:44926/148‑1 (MQ=255)
                                                                                                                 gaCGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGGCCACGACGAAGTTGAcc                                 <  8:211608/147‑1 (MQ=255)
                                                                                                                   cGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGa                              <  2:276525/149‑1 (MQ=255)
                                                                                                                    ggggAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTAGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGaa                             <  5:211518/149‑1 (MQ=255)
                                                                                                                      ggAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAAt                            >  2:77000/1‑148 (MQ=255)
                                                                                                                       gAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCa                          >  4:132728/1‑149 (MQ=255)
                                                                                                                         ggggTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTTACACGGATTTTGGAATCGTGCTA‑GTCATGGCAAGGGCCTCAGCGCTTATGACGGATGCGTGGGTTGGCCAGGCCATAGAGGATGTCCACGACGAAGTTGACCAGAATCAcc                        <  6:139552/149‑1 (MQ=255)
                                                                                                                           ggTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACcag                      <  1:77000/149‑1 (MQ=255)
                                                                                                                                  cAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCaggcaggc                >  4:199962/1‑148 (MQ=255)
                                                                                                                                   aGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGa              >  5:40341/1‑149 (MQ=255)
                                                                                                                                    ggCTCTGGTCCACCGGGCTGGGCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGAt             >  6:84750/1‑149 (MQ=255)
                                                                                                                                          gtgtcACCGGCCAGCCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACAg       >  3:78297/5‑149 (MQ=255)
                                                                                                                                           tgtcACCGGCCAGCCCGGTGGCACGGATTTTGGAATCGGGCTAGGTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGcc                                                             >  1:214542/4‑95 (MQ=255)
                                                                                                                                           tgtcACCGGCCAGCCCGGTGGCACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGcc                                                             <  2:214542/91‑1 (MQ=255)
                                                                                                                                           gTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACa        <  1:266220/147‑1 (MQ=255)
                                                                                                                                           gTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACAgg      >  5:22859/1‑149 (MQ=255)
                                                                                                                                               aCCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACAGGATGc  >  8:40388/1‑149 (MQ=255)
                                                                                                                                                 cGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCGGTAGAGGATGTCCACGACGAAGTTGAc                                  <  1:150061/115‑1 (MQ=255)
                                                                                                                                                 cGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCAGTAGAGGATGTCCACGACGAAGTTGAc                                  >  2:150061/1‑115 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTCATGCGGGGCGACCCAAGGCGCGAACAGCGCGCAGAAC‑ACCACCAGGCACATGAAGGCCAGGCCCATGACCGCGCCCTTGTTGCGCGCGAAGGCTTGCCAGAATTCTTTGTACGGGGAGGGGTAGAGCAGGCTCTGGTCCACCGGGCTGGCCGGTGACACGGATTTTGGAATCGGGCTA‑GTCATGGCAAGGGCCTCAGCGCTGATGACGGATGCGTGGGTTGGCCAGGCCGTAGAGGATGTCCACGACGAAGTTGACCAGAATCACCAGGCAGGCGATCAACAGGATGC  >  NC_002947/1020780‑1021070

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: ‑