Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 616,669 C→T 6.2% T838I (ACC→ATC)  entF → enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913616,6690CT 6.2% 207.0 / 3.0 64T838I (ACC→ATC) entFenterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base T (2/2);  total (33/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 = 9.28e-01

GATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCG  >  NC_000913/616523‑616812
                                                                                                                                                  |                                                                                                                                               
gatgatGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATgtaccgtacc                                                                                                                                                <  5:78383/148‑1 (MQ=255)
   gatgCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGtaccgtaccg                                                                                                                                               >  5:336798/1‑146 (MQ=255)
    atgCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGgag                                                                                                                                            <  8:337645/148‑1 (MQ=255)
    atgCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGgag                                                                                                                                            <  6:336798/148‑1 (MQ=255)
           cGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCCGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTg                                                                                                                                      <  5:49628/147‑1 (MQ=255)
            ggTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGcc                                                                                                                                    >  6:214668/1‑148 (MQ=255)
                   ccgGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCgctggc                                                                                                                             <  5:214668/148‑1 (MQ=255)
                   ccgGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCgctggc                                                                                                                             <  8:196061/148‑1 (MQ=255)
                                 tGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGcgcg                                                                                                               <  7:143258/148‑1 (MQ=255)
                                        tATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCGCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTg                                                                                                            <  6:10185/144‑1 (MQ=255)
                                        tATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTgg                                                                                                           <  5:162161/145‑1 (MQ=255)
                                        tATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTgg                                                                                                           <  8:213245/145‑1 (MQ=255)
                                             cACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCACCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTc                                                                                                   >  5:361956/1‑148 (MQ=255)
                                               cTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCgg                                                                                                 >  7:106092/1‑148 (MQ=255)
                                                 ggCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGgc                                                                                               >  5:36297/1‑148 (MQ=255)
                                                      tCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGgag                                                                                                                                            <  8:372282/98‑1 (MQ=255)
                                                      tCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGgag                                                                                                                                            >  7:372282/1‑98 (MQ=255)
                                                      tCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGt                                                                                          >  1:260503/1‑148 (MQ=255)
                                                        aaCTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtga                                                                                        <  7:251374/148‑1 (MQ=255)
                                                            ggCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTATCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTgatgat                                                                                    >  3:147321/1‑148 (MQ=255)
                                                            ggCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTATCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTgatgat                                                                                    >  3:141887/1‑148 (MQ=255)
                                                            ggCGCAGGGCTATCTCGGACGCCACGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTgatgat                                                                                    >  2:230258/1‑148 (MQ=255)
                                                             gcgcAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtgatga                                                                                     <  6:18970/146‑1 (MQ=255)
                                                                    gCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTaa                                                                             >  6:122362/1‑147 (MQ=255)
                                                                       aTCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAAt                                                                         >  2:44971/1‑148 (MQ=255)
                                                                         ctcGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCTCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTc                                                                       >  1:114696/1‑148 (MQ=255)
                                                                           cGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGt                                                                                                             <  2:301809/108‑1 (MQ=255)
                                                                           cGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGt                                                                                                             >  1:301809/1‑108 (MQ=255)
                                                                           cGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCTGTGATGATCAGCTAAAAATTcgc                                                                     >  6:269409/1‑148 (MQ=255)
                                                                            ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcgcg                                                                    >  4:219192/1‑148 (MQ=255)
                                                                            ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcgc                                                                     <  5:122362/147‑1 (MQ=255)
                                                                            ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcgc                                                                     <  5:25180/147‑1 (MQ=255)
                                                                            ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcg                                                                      <  6:165252/146‑1 (MQ=255)
                                                                            ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcg                                                                      <  3:294788/146‑1 (MQ=255)
                                                                                  cccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGc                                                              >  2:179896/1‑148 (MQ=255)
                                                                                   ccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCTCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCg                                                             <  2:114696/148‑1 (MQ=255)
                                                                                   ccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCg                                                             >  7:183512/1‑148 (MQ=255)
                                                                                    cGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTgg                                                                                                           <  5:297897/101‑1 (MQ=255)
                                                                                    cGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTgg                                                                                                           >  6:297897/1‑101 (MQ=255)
                                                                                    cGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtgatga                                                                                     >  3:342023/1‑123 (MQ=255)
                                                                                    cGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtgatga                                                                                     <  4:342023/123‑1 (MQ=255)
                                                                                          gACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGaa                                                      >  4:53884/1‑148 (MQ=255)
                                                                                           aCCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTATCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGt                                                            <  4:147321/141‑1 (MQ=255)
                                                                                           aCCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTATCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGt                                                            <  4:141887/141‑1 (MQ=255)
                                                                                                       tttATTGCCTATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGgc                                                                                                                  <  8:111231/75‑1 (MQ=255)
                                                                                                       tttATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGgc                                                                                                                  >  7:111231/1‑75 (MQ=255)
                                                                                                       tttATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAatc                                          >  2:330542/1‑147 (MQ=255)
                                                                                                              ccGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATcgcg                                   <  8:183512/147‑1 (MQ=255)
                                                                                                                gATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGa                                >  8:67583/1‑148 (MQ=255)
                                                                                                                 aTCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTAt                                                          >  1:286028/1‑121 (MQ=255)
                                                                                                                 aTCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTAt                                                          <  2:286028/121‑1 (MQ=255)
                                                                                                                 aTCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGAt                               >  3:235662/1‑148 (MQ=255)
                                                                                                                    cTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATa                                                                                                                       >  8:292122/1‑57 (MQ=255)
                                                                                                                    cTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATa                                                                                                                       <  7:292122/57‑1 (MQ=255)
                                                                                                                      tttGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGCTCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAgg                          >  7:316169/1‑148 (MQ=255)
                                                                                                                           cccAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTg                     <  6:205328/148‑1 (MQ=255)
                                                                                                                              aGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCg                  <  1:179896/148‑1 (MQ=255)
                                                                                                                                gTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGa                                                       <  5:246951/109‑1 (MQ=255)
                                                                                                                                gTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGa                                                       >  6:246951/1‑109 (MQ=255)
                                                                                                                                 tGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATcgc                                    <  7:118864/127‑1 (MQ=255)
                                                                                                                                 tGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATcgc                                    >  8:118864/1‑127 (MQ=255)
                                                                                                                                 tGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGAt               >  5:147125/1‑148 (MQ=255)
                                                                                                                                  gAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCg                  >  4:189924/1‑144 (MQ=255)
                                                                                                                                  gAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCg                  <  3:189924/144‑1 (MQ=255)
                                                                                                                                  gAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATg              >  8:185521/1‑148 (MQ=255)
                                                                                                                                    aCGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTc            >  4:216903/1‑148 (MQ=255)
                                                                                                                                           taccgtaccgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACaa      <  3:219192/147‑1 (MQ=255)
                                                                                                                                             ccgtaccgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGc    <  4:50668/147‑1 (MQ=255)
                                                                                                                                              cgtaccgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCg  <  1:44971/148‑1 (MQ=255)
                                                                                                                                                 accgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGaa         >  3:163831/1‑138 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCG  >  NC_000913/616523‑616812

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