Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,077 G→A 14.3% intergenic (‑113/+15) mfd ← / ← ycfT transcription‑repair coupling factor/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,0770GA14.3% 243.0 / 30.6 91intergenic (‑113/+15)mfd/ycfTtranscription‑repair coupling factor/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base A (7/6);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

GGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACC  >  NC_000913/1173936‑1174224
                                                                                                                                             |                                                                                                                                                   
ggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                <  4:947551/147‑1 (MQ=255)
ggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                >  3:947551/1‑147 (MQ=255)
     cGTATAACGATATTGTTCAGGCATAACGGTTCTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGAt                                                                                                                                          >  4:975794/1‑148 (MQ=255)
              atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGa                                                                                                                                           >  1:959662/1‑138 (MQ=255)
              atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGa                                                                                                                                           <  2:959662/138‑1 (MQ=255)
                atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCg                                                                                                                                            >  1:108664/1‑135 (MQ=255)
                atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCg                                                                                                                                            <  2:108664/135‑1 (MQ=255)
                  tGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAg                                                                                                                             <  3:975794/148‑1 (MQ=255)
                            tAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                <  4:573308/119‑1 (MQ=255)
                            tAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                >  3:573308/1‑119 (MQ=255)
                             taCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                   <  4:1000739/114‑1 (MQ=255)
                             aaCGGTTGTCTGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAAGCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                   >  3:1000739/1‑115 (MQ=255)
                                        gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTg                                                                                                       >  4:392759/1‑148 (MQ=255)
                                          ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                   >  2:553936/1‑102 (MQ=255)
                                          ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                   <  1:553936/102‑1 (MQ=255)
                                                 tAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                        <  3:480588/106‑1 (MQ=255)
                                                 tAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                        >  4:480588/1‑106 (MQ=255)
                                                 tAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCg                                                                                                 >  4:923882/1‑145 (MQ=255)
                                                 tAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCg                                                                                                 <  3:923882/145‑1 (MQ=255)
                                                     atatGCCCCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                        >  1:766406/1‑102 (MQ=255)
                                                     atatGCCCCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATgg                                                                                                                                        <  2:766406/102‑1 (MQ=255)
                                                       atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaaca                                                                                        >  4:959635/1‑148 (MQ=255)
                                                          ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                <  2:562198/137‑1 (MQ=255)
                                                          ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                >  1:562198/1‑137 (MQ=255)
                                                             cAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           >  3:544274/1‑139 (MQ=255)
                                                             cAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           <  4:544274/139‑1 (MQ=255)
                                                             cAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGTCGGGGCGa                                                                                                                                           >  4:650264/1‑91 (MQ=255)
                                                             cAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                           <  3:650264/91‑1 (MQ=255)
                                                                    tGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGaa                                                                           <  3:959635/148‑1 (MQ=255)
                                                                       gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                           >  2:58373/1‑81 (MQ=255)
                                                                       gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                           <  1:58373/81‑1 (MQ=255)
                                                                        ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGAc                                                                                                                     <  3:224660/102‑1 (MQ=255)
                                                                        ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGAc                                                                                                                     >  4:224660/1‑102 (MQ=255)
                                                                        ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                               <  4:906290/124‑1 (MQ=255)
                                                                        ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                               >  3:906290/1‑124 (MQ=255)
                                                                        ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCa                                                                       >  2:525466/1‑148 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         >  4:218979/1‑75 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         >  2:721414/1‑75 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         <  3:218979/75‑1 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         <  1:721414/75‑1 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCtt                                                                                                                                     <  4:437572/79‑1 (MQ=255)
                                                                               ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCtt                                                                                                                                     >  3:437572/1‑79 (MQ=255)
                                                                                tctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGa                                                                                                      >  4:396685/1‑109 (MQ=255)
                                                                                tctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGa                                                                                                      <  3:396685/109‑1 (MQ=255)
                                                                                 ctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      >  4:427956/1‑124 (MQ=255)
                                                                                 ctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      <  3:427956/124‑1 (MQ=255)
                                                                                       aGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                        <  3:392759/148‑1 (MQ=255)
                                                                                           aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           >  4:814255/1‑109 (MQ=255)
                                                                                           aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           <  3:814255/109‑1 (MQ=255)
                                                                                           aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                    <  1:525466/148‑1 (MQ=255)
                                                                                                  tcatcaTCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTg                                                                                                                                                 >  2:345900/1‑48 (MQ=255)
                                                                                                  tcatcaTCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTg                                                                                                                                                 <  1:345900/48‑1 (MQ=255)
                                                                                                  tCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      >  2:735495/1‑107 (MQ=255)
                                                                                                  tCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      <  1:735495/107‑1 (MQ=255)
                                                                                                     ttATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                          >  1:96112/1‑148 (MQ=255)
                                                                                                       aTCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCg                                        <  2:96112/148‑1 (MQ=255)
                                                                                                            tAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         <  1:391418/46‑1 (MQ=255)
                                                                                                            tAAACGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                         >  2:391418/1‑46 (MQ=255)
                                                                                                                            gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      >  3:34216/1‑81 (MQ=255)
                                                                                                                            gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                      <  4:34216/81‑1 (MQ=255)
                                                                                                                               aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                >  4:714944/1‑148 (MQ=255)
                                                                                                                                aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           >  4:942/1‑72 (MQ=255)
                                                                                                                                aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                           <  3:942/72‑1 (MQ=255)
                                                                                                                                aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                   <  4:389108/128‑1 (MQ=255)
                                                                                                                                aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                   >  3:389108/1‑128 (MQ=255)
                                                                                                                                     aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                    >  4:151061/1‑106 (MQ=255)
                                                                                                                                     aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                    <  3:151061/106‑1 (MQ=255)
                                                                                                                                     aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAcc                                                  >  2:93285/1‑108 (MQ=255)
                                                                                                                                     aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAcc                                                  <  1:93285/108‑1 (MQ=255)
                                                                                                                                      tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                   >  3:649514/1‑58 (MQ=255)
                                                                                                                                      tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                   <  4:649514/58‑1 (MQ=255)
                                                                                                                                      tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg         >  3:955980/1‑148 (MQ=255)
                                                                                                                                      tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGCGAACAACAGATCGc                                                                               <  2:612273/78‑1 (MQ=255)
                                                                                                                                      tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGc                                                                               >  1:612273/1‑78 (MQ=255)
                                                                                                                                      tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTg                 >  4:1062897/1‑140 (MQ=255)
                                                                                                                                      tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTg                 <  3:1062897/140‑1 (MQ=255)
                                                                                                                                       cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGa                                                                                                      >  3:392678/1‑54 (MQ=255)
                                                                                                                                       cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGa                                                                                                      <  4:392678/54‑1 (MQ=255)
                                                                                                                                       cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGaa                                                              >  2:574865/1‑94 (MQ=255)
                                                                                                                                       cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGaa                                                              <  1:574865/94‑1 (MQ=255)
                                                                                                                                       cGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCt                                                      >  3:284365/1‑102 (MQ=255)
                                                                                                                                       cGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCt                                                      <  4:284365/102‑1 (MQ=255)
                                                                                                                                       cGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa        >  4:813927/1‑148 (MQ=255)
                                                                                                                                         ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGc                                                                        >  4:1067878/1‑82 (MQ=255)
                                                                                                                                         ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGc                                                                        <  3:1067878/82‑1 (MQ=255)
                                                                                                                                          cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGTATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATaa     >  4:906512/1‑148 (MQ=255)
                                                                                                                                          cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCt                                                      >  4:888339/1‑99 (MQ=255)
                                                                                                                                          cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCt                                                      <  3:888339/99‑1 (MQ=255)
                                                                                                                                          cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATaa     <  4:955980/148‑1 (MQ=255)
                                                                                                                                           ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                <  3:51919/120‑1 (MQ=255)
                                                                                                                                           ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                >  4:51919/1‑120 (MQ=255)
                                                                                                                                             ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGTATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAAcc  <  3:906512/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACC  >  NC_000913/1173936‑1174224

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