Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,215 A→G 8.2% intergenic (+21/‑209) 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,2150AG 8.2% 154.6 / 3.8 64intergenic (+21/‑209)mtlA/mtlDmannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (36/22);  new base G (4/2);  total (40/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

ACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACG  >  NC_000913/3774073‑3774352
                                                                                                                                              |                                                                                                                                         
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)
                                                                                                                                              |                                                                                                                                         
ACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCTCCTCACG  >  NC_000913/3774073‑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: ‑