Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,214 C→T 6.8% intergenic (+20/‑210) mtlA → / → mtlD mannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,774,2140CT 6.8% 158.8 / 3.3 64intergenic (+20/‑210)mtlA/mtlDmannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (38/20);  new base T (2/2);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-02

AACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACG  >  NC_000913/3774072‑3774352
                                                                                                                                              |                                                                                                                                          
aacGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCCGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCCCcac                                                                                                                                            >  2:71435/1‑143 (MQ=255)
 acGAGCGCATTCAGGTTGTCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCAGCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAAGCCCACCCTCTCcaca                                                                                                                                           <  2:201626/143‑1 (MQ=255)
     gCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACCCACCACCAGCGTCGATGAAGTGCTTGAACTGCTGGCAGGTCGTCAGTAATCCAATCCCACCCTCTCCACATCga                                                                                                                                       >  2:162938/1‑143 (MQ=255)
     gCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACACTCTCCACGTGga                                                                                                                                       >  2:429047/1‑143 (MQ=255)
         aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATggcgcag                                                                                                                                   >  2:217116/1‑138 (MQ=255)
         aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                   >  1:159977/1‑143 (MQ=255)
         aTTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTCGCAGGCCGTTTGTTTGCCAATCCCACCCTCTCCACATGGAGAAg                                                                                                                                   >  1:194845/1‑143 (MQ=255)
           tCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                 >  2:415751/1‑143 (MQ=255)
           tCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGt                                                                                                                                 >  2:401324/1‑143 (MQ=255)
               gTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTggga                                                                                                                             <  2:85669/143‑2 (MQ=255)
                    cACCAGCCTGAGCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                        <  1:266481/143‑1 (MQ=255)
                    cACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGCGAAGGTGGGCTTAAt                                                                                                                        >  1:368706/1‑143 (MQ=255)
                    cACCAGCCTGACCAATGCACTGGATGATGAGTCCGCCAGCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAAt                                                                                                                        <  2:96378/143‑1 (MQ=255)
                    cACCAGCCTGACCAATGCACGGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCTTCTCCATGTGGAGAGGGTgggattgga                                                                                                                        <  1:259363/143‑7 (MQ=38)
                              aCCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGCGGGGTTAATGGCCGGATgc                                                                                                              >  1:54275/1‑143 (MQ=255)
                                  aTGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGCGCTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCCCATGGAGGTGGTGGGGTGAATTGCCTGATGCGCtg                                                                                                          >  1:314894/1‑142 (MQ=255)
                                     cattagatggtGAGTGCGTGATCGCGCGTCTGGCACACCCGACCAGCGTGGATGAAGTGCTTGAACTGCTGGCAGGTGGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGc                                                                                                       <  2:180712/138‑1 (MQ=255)
                                              atgaGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTACGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGc                                                                                              >  1:30979/1‑143 (MQ=255)
                                                   tcggTCAGCGCGCGTGTGGCACACACCACCAGCGTGGATGAAGTGCTCGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                         <  2:293788/140‑1 (MQ=39)
                                                   tCCGTCATCGAGCGGCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                                                                         <  1:420752/143‑1 (MQ=255)
                                                   tCCGGCATCGGGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGTGTTAATTGCCTGATTCGCTACGCTTATCAGGCCTAcc                                                                                         <  2:360329/143‑2 (MQ=255)
                                                     cGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAgg                                                                                       >  1:331644/1‑143 (MQ=255)
                                                      ggcagggAGCGTCTGGCCCACCGTCCCAGCGGGTCTGAAGGGCGGGAACTGCTGGCAGGTCGTAAGGAATCCAATCCCGCCCTGTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                      <  2:194689/137‑1 (MQ=25)
                                                      gTCATGGCGCGTCGCGCACACAACACTAGCGTGGATGAAGGGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAGTCCCACCCTCTCCGCATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGa                                                                                      <  2:42384/143‑1 (MQ=38)
                                                       tCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAggct                                                                                     >  1:388408/1‑141 (MQ=255)
                                                        cATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATg                                                                                    <  1:113077/143‑1 (MQ=255)
                                                           cGATCGTCTGGCACACACCACCAGCGTGGATGAAGCGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACAGGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAt                                                                                 <  1:130278/143‑1 (MQ=255)
                                                           cGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCAt                                                                                 <  1:287579/143‑1 (MQ=255)
                                                             aGCCTCTGGCACACACCCCCAGCGTGAACGAACGGCTGGAACTGCTGCCCCGTCGCAACTAATCCAATCCCACCCCCTCCACATGCAGAAGTTGGAGATAACTCCCTAATGCGCTACGCTTATCAGGCCTACAGGATGCATca                                                                               <  2:33262/143‑1 (MQ=255)
                                                                 tCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCGGGCGTCCGGGATGCATCACCAt                                                                           >  2:222266/1‑143 (MQ=255)
                                                                     gCACAGACCACCAGCGTGGATGAAGTGCGGGAACTGCTGGCAGGGCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                       <  1:280103/143‑1 (MQ=255)
                                                                     gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCGGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttg                                                                        >  1:110279/1‑142 (MQ=255)
                                                                     gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCGGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttg                                                                        <  2:110279/142‑1 (MQ=255)
                                                                     gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGCAGGTGGGGTTAATTAACTGATGCGATCAGATTTTCAGCCCGACAGGATGCATCACAATttgt                                                                       >  1:318531/1‑143 (MQ=255)
                                                                     gCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATttgt                                                                       >  1:369248/1‑143 (MQ=255)
                                                                       acacacCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCAACCCTGTCCACATGGAGAAGGTGGGGCTAATTGCCTGACCCGCCCGGTCGGTCTGGGCTACAGGATGCATCACAACttgttg                                                                     >  2:205261/1‑143 (MQ=255)
                                                                                          gAAGTGCTGTAACTGCTGACAGGCCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGACGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAg                                                  <  1:268704/143‑1 (MQ=255)
                                                                                          gAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAg                                                  <  2:159977/143‑1 (MQ=255)
                                                                                            aGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGTAGGTGGGGTTAATTACCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGc                                                >  1:128400/1‑143 (MQ=255)
                                                                                              tGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATCTGTTGAGTTTGCACGTTCTTGTAGGCCg                                              >  2:341912/1‑143 (MQ=255)
                                                                                               gCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGACTTTGCACGTTCTTGTAGGCCgg                                             >  1:76183/1‑143 (MQ=255)
                                                                                               gCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCgg                                             >  1:298719/1‑143 (MQ=255)
                                                                                                 tGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAt                                           >  2:356570/1‑143 (MQ=255)
                                                                                                   gAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACGTGGAGAAGGTGGGGTGAATTGCCTGAGGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATGTGCACGTTCTCGTAGGCCggttga                                         >  1:69630/1‑139 (MQ=38)
                                                                                                       tgcttGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGTTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGgcg                                     <  2:193566/143‑1 (MQ=255)
                                                                                                          tgGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGAGAAGGCGCtt                                  >  1:312309/1‑143 (MQ=255)
                                                                                                           ggCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGTTCTGGCTCTTa                                 >  1:212508/1‑143 (MQ=255)
                                                                                                              aGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGTCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgc                              <  1:98609/143‑1 (MQ=255)
                                                                                                                gTCGTAAGTAATCCAATCCCACCCTCTCCCCATGGATATGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTCGGCCGGATAGGGCACTCAcgccg                            >  2:119965/1‑143 (MQ=37)
                                                                                                                 tCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgccgc                           >  1:84850/1‑143 (MQ=255)
                                                                                                                 tCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgccgc                           >  1:319329/1‑143 (MQ=255)
                                                                                                                   gtgagtaaTCCGATCCCACCCTCTGCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTCGTTGAATGTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCAt                         <  2:148276/140‑1 (MQ=37)
                                                                                                                   gtaagtaaTCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTAcgccgcct                         >  2:10660/1‑141 (MQ=255)
                                                                                                                        taaTCCAATCCCACCCTCTCCACATGGAGAAGGT‑‑GGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCcgccgc                  >  2:82239/1‑143 (MQ=255)
                                                                                                                         aatccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATGTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGgcg                   <  1:42447/143‑1 (MQ=255)
                                                                                                                           tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCt                 >  1:283919/1‑143 (MQ=255)
                                                                                                                           tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCt                 >  2:55608/1‑143 (MQ=255)
                                                                                                                           tccaatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATATGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCt                 >  2:53554/1‑143 (MQ=255)
                                                                                                                              aatccCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGATGAAATTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcc              >  1:468303/1‑143 (MQ=255)
                                                                                                                              aatccCACCCTCTCCACATGGAGAAGGTGGGGGTAATTGCCTGATGCGCTACGCTTGTCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAAGCCGGATAAGGCGCCTACGCCACATCCGGCGGTGcc              >  2:124463/1‑143 (MQ=37)
                                                                                                                               cgcccACCCGCTCGCGATGGAGAGTGTGGGGTTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGccc             <  1:250425/141‑1 (MQ=21)
                                                                                                                                 cccACCTTCTCCATGTGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCCTCCGGCGCTGCCCct           >  2:192064/1‑143 (MQ=25)
                                                                                                                                  ccACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCACGCCTACAGGATGCATCACAGTTTGTTGAATTTGCACGCTCTTGTAAGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGcccccc          >  2:208823/1‑141 (MQ=255)
                                                                                                                                   cACCTTCTCCATGTGGAGAGGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTCACGCCGCATCCGGCGCTGCCCctct         >  2:323301/1‑143 (MQ=25)
                                                                                                                                         ttcCAGATGGAGAAGGTGGGGTTAAGTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCAc   <  1:134758/142‑1 (MQ=255)
                                                                                                                                          tcCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACg  >  1:328440/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                          
AACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACG  >  NC_000913/3774072‑3774352

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