Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,533,129 T→C L994P (CTG→CCG)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,533,1290TC98.4% 218.1 / ‑4.4 62L994P (CTG→CCG) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (27/33);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

GTGGGCCGGCAACGGCCTCAAACCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCTGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAACGC  >  NC_002947/3532996‑3533261
                                                                                                                                     |                                                                                                                                    
gtggGCCGGCAACGGCCTCAAACCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAgccgc                                                                                                                                    >  2:1067488/1‑136 (MQ=255)
      cGGCAACGGCCTCAAACCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAgg                                                                                                                              >  1:151587/1‑136 (MQ=255)
          aaCGGCCTCAATCCGGCTGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCTTTGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCTGCGACGGGCGCTGGTGTGTTCGCAGCAGCTCAAGGTCAg                                                                                                                          <  1:1368270/136‑1 (MQ=255)
                    aaCCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTg                                                                                                                <  1:1021425/136‑1 (MQ=255)
                    aaCCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTg                                                                                                                <  1:10529/136‑1 (MQ=255)
                    aaCCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTg                                                                                                                <  1:586600/136‑1 (MQ=255)
                      cggCCAGCGGTGTGACCTTCGTTTCTGTGGAAGACGAACATGGCTTGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGa                                                                                                              <  2:1019161/134‑1 (MQ=255)
                      ccGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGa                                                                                                              <  2:1036157/136‑1 (MQ=255)
                      ccGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGa                                                                                                              <  2:885075/136‑1 (MQ=255)
                            gCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAggt                                                                                                        >  2:1267254/1‑135 (MQ=255)
                            gCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGga                                                                                                        <  1:1192301/136‑1 (MQ=255)
                              ggTGTGACCTTCGTTACTCTGGAAGACGAACCTGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCCGCGACGGGCGCTGGTGGGTGCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGgaga                                                                                                      <  2:1154215/136‑1 (MQ=255)
                                            tACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCAc                                                                                        >  1:1053827/1‑136 (MQ=255)
                                                   gAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTg                                                                                 <  2:790072/136‑1 (MQ=255)
                                                    aaGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgc                                                                                >  1:998226/1‑136 (MQ=255)
                                                     aGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               <  2:272037/136‑1 (MQ=255)
                                                     aGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               <  1:347833/136‑1 (MQ=255)
                                                       aCGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGcgcg                                                                             <  1:354046/136‑1 (MQ=255)
                                                         gAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGTACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCTCGTGATTGCgcggc                                                                           <  1:42755/136‑1 (MQ=255)
                                                         gAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAgccgc                                                                                                                                    <  2:218900/79‑1 (MQ=255)
                                                         gAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAgccgc                                                                                                                                    >  1:218900/1‑79 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               >  2:629558/1‑130 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               >  2:1145474/1‑130 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               >  2:635785/1‑130 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               <  1:1145474/130‑1 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               <  1:635785/130‑1 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTgcg                                                                               <  1:629558/130‑1 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGcggcgg                                                                         >  2:886188/1‑136 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGcggcgg                                                                         <  1:425201/136‑1 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGcggcgg                                                                         >  1:1268883/1‑136 (MQ=255)
                                                           aCATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGcggcgg                                                                         <  2:88575/136‑1 (MQ=255)
                                                             atggcatggTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGtt                                                                       <  2:933542/136‑1 (MQ=255)
                                                                  atggTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTggagg                                                                  >  2:91339/1‑136 (MQ=255)
                                                                        aaCGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGa                                                            >  1:1040431/1‑136 (MQ=255)
                                                                         aCGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAg                                                           >  2:1144072/1‑136 (MQ=255)
                                                                          cGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGt                                                          <  1:1019275/136‑1 (MQ=255)
                                                                                  tgtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  >  1:955782/1‑136 (MQ=255)
                                                                                  tgtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  <  1:257653/136‑1 (MQ=255)
                                                                                  tgtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  >  1:230575/1‑136 (MQ=255)
                                                                                  tgtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGCGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  >  1:967448/1‑136 (MQ=255)
                                                                                  ggtgGCGTGCGCTGGCCGATCGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  <  2:74229/135‑1 (MQ=255)
                                                                                  ggtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  <  1:408098/135‑1 (MQ=255)
                                                                                  ggtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  <  1:349120/135‑1 (MQ=255)
                                                                                  ggtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCgctgct                                                  <  2:1196991/135‑1 (MQ=255)
                                                                                   gtgGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGctgctg                                                 >  2:1345533/1‑136 (MQ=255)
                                                                                     ggCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCa                                               >  2:1206305/1‑136 (MQ=255)
                                                                                       cGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCTGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAg                                             <  2:21296/136‑1 (MQ=255)
                                                                                       cGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAg                                             >  2:1340033/1‑136 (MQ=255)
                                                                                       cGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAg                                             >  2:964612/1‑136 (MQ=255)
                                                                                         tGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAggg                                           >  1:1321996/1‑136 (MQ=255)
                                                                                         tGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAggg                                           >  1:1157343/1‑136 (MQ=255)
                                                                                                 ccGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGAtgt                                   <  1:1010115/136‑1 (MQ=255)
                                                                                                   gAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTgc                                 >  2:1149180/1‑136 (MQ=255)
                                                                                                        gCAGCGACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGc                            >  1:1336344/1‑136 (MQ=255)
                                                                                                             gACGGGCGCTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGa                       >  1:225594/1‑136 (MQ=255)
                                                                                                                    gcTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCAc                >  2:1287178/1‑136 (MQ=255)
                                                                                                                    gcTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCAc                <  1:837181/136‑1 (MQ=255)
                                                                                                                    gcTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCAc                <  1:673890/136‑1 (MQ=255)
                                                                                                                    gcTGGTGGGTTCGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCAc                <  1:655921/136‑1 (MQ=255)
                                                                                                                              gcGCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCa      <  1:173893/135‑1 (MQ=255)
                                                                                                                                gCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAAc    <  1:908708/136‑1 (MQ=255)
                                                                                                                                gCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAAc    >  1:635844/1‑136 (MQ=255)
                                                                                                                                gCAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAAc    <  1:1346517/136‑1 (MQ=255)
                                                                                                                                 cAGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAACg   <  2:1240434/136‑1 (MQ=255)
                                                                                                                                  aGCCGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGAGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAACGc  >  1:1107113/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                    
GTGGGCCGGCAACGGCCTCAAACCGCCAGCGGTGTGACCTTCGTTACTCTGGAAGACGAACATGGCATGGTCAACGTGGTGGTGTGGCGTGCGCTGGCCGAACGGCAGCGACGGGCGCTGGTGGGTTCGCAGCTGCTCAAGGTCAGTGGGCGGTTGGAGCAGGAGAATGGGGTGCGGCACCTGATTGCGCGGCGGTTGGAGGATGTGAGTCCGCTGCTGCAAGGGCTGGATGTGCGCAGCCGGGATTTTCACTGAAACCCCAACGC  >  NC_002947/3532996‑3533261

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