Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,774,207 C→T 6.0% intergenic (+13/‑217) 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,2070CT 6.0% 173.1 / 3.1 67intergenic (+13/‑217)mtlA/mtlDmannitol‑specific PTS enzyme IICBA component/mannitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (39/24);  new base T (2/2);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-02

AACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCT  >  NC_000913/3774069‑3774345
                                                                                                                                          |                                                                                                                                          
aacaacGAGCACATTCATGTTATCACCAGCCTGGCCAATGCACTGGGTGATGAGTCCGTCATCGAGCGTCTGGCACACACCAGCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCtctc                                                                                                                                        <  2:117592/143‑1 (MQ=255)
  ccacGACCACATGCAGGTGATCACCAGGCTGACCGTTGCCCTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCACCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAACCCCACCCTCTCca                                                                                                                                      <  2:7606/141‑1 (MQ=255)
  caacGAGAACATTCAGGTTATCACCAGCCTGACCAATGCACTGCACGATGAGTCCGTCATCGAGCGTCTGGCACACAGCACCAGCGTGAATGAAGNGCTGGAACTACTGGCAGGTCGTAAGTAATCCACTCCTACCCTCTCca                                                                                                                                      >  2:18658/1‑143 (MQ=255)
   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)
                                                                                                                                          |                                                                                                                                          
AACAACGAGCACATTCAGGTTATCACCAGCCTGACCAATGCACTGGATGATGAGTCCGTCATCGAGCGTCTGGCACACACCACCAGCGTGGATGAAGTGCTGGAACTGCTGGCAGGTCGTAAGTAATCCAATCCCACCCTCTCCACATGGAGAAGGTGGGGTTAATTGCCTGATGCGCTACGCTTATCAGGCCTACAGGATGCATCACAATTTGTTGAATTTGCACGTTCTTGTAGGCCGGATAAGGCGCTTACGCCGCATCCGGCGCTGCCCCTCT  >  NC_000913/3774069‑3774345

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: