Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,794,048 T→C intergenic (+33/‑205) gabT → / → gabP 4‑aminobutyrate aminotransferase GabT/4‑aminobutyrate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,0480TC82.6% 720.5 / 136.5 270intergenic (+33/‑205)gabT/gabP4‑aminobutyrate aminotransferase GabT/4‑aminobutyrate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (18/29);  new base C (121/102);  total (139/131)
Fisher's exact test for biased strand distribution p-value = 5.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAG  >  NC_000913/2793950‑2794145
                                                                                                  |                                                                                                 
gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCt                                                                                                 <  1:4546089/101‑1 (MQ=255)
   gACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGAGCCGCTCCTAGGCCGGAGGCGACGCTGCGCGCCTTGt                                                                                              >  1:4546505/1‑101 (MQ=255)
    aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTc                                                                                             <  1:4546270/101‑1 (MQ=255)
    aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTc                                                                                             >  2:4546506/1‑101 (MQ=255)
     cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTcc                                                                                            >  1:4546510/1‑101 (MQ=255)
      gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCg                                                                                           <  2:4546010/101‑1 (MQ=255)
      gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCg                                                                                           >  2:4546511/1‑101 (MQ=255)
      gCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCg                                                                                           <  2:4546088/101‑1 (MQ=255)
         cAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGcc                                                                                        >  2:4546516/1‑101 (MQ=255)
          aGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCt                                                                                       <  2:4546144/101‑1 (MQ=255)
          aGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCt                                                                                       <  2:4546182/101‑1 (MQ=255)
           gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTa                                                                                      >  1:4546520/1‑101 (MQ=255)
           gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTa                                                                                      <  2:4546187/101‑1 (MQ=255)
           gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTa                                                                                      <  1:4546206/101‑1 (MQ=255)
           gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTa                                                                                      <  2:4546207/101‑1 (MQ=255)
            aTCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTAc                                                                                     <  2:4546050/101‑1 (MQ=255)
            aTCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTAc                                                                                     <  1:4546062/101‑1 (MQ=255)
             tCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                    >  1:4546528/1‑101 (MQ=255)
             tCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACg                                                                                    >  2:4546526/1‑101 (MQ=255)
               cGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACggg                                                                                  <  1:4545888/101‑1 (MQ=255)
                gTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACgggg                                                                                 >  1:4546532/1‑101 (MQ=255)
                 tCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGa                                                                                <  2:4546158/101‑1 (MQ=255)
                  cAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGAt                                                                               <  1:4545758/101‑1 (MQ=255)
                   aGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTc                                                                                                  >  2:4546540/1‑81 (MQ=255)
                   aGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATc                                                                              >  1:4546537/1‑101 (MQ=255)
                     ggTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCa                                                                            <  2:4545979/101‑1 (MQ=255)
                     ggTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCa                                                                            >  1:4546541/1‑101 (MQ=255)
                      gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAg                                                                           <  2:4546101/101‑1 (MQ=255)
                        cTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGt                                                                         >  1:4546545/1‑101 (MQ=255)
                        cTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGt                                                                         >  2:4546546/1‑101 (MQ=255)
                         tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTc                                                                        <  1:4545994/101‑1 (MQ=255)
                           gagaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCgg                                                                      >  2:4546551/1‑101 (MQ=255)
                           gagaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCgg                                                                      >  1:4546552/1‑101 (MQ=255)
                           gagaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCgg                                                                      <  2:4546243/101‑1 (MQ=255)
                            agaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGa                                                                     <  2:4546162/101‑1 (MQ=255)
                             gaTCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTg                                                                                               >  2:4546555/1‑74 (MQ=255)
                               tcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATaa                                                                  >  1:4546558/1‑101 (MQ=255)
                                catcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAg                                                                 >  1:4546561/1‑101 (MQ=255)
                                  tcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGc                                                               <  2:4545927/101‑1 (MQ=255)
                                   caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCg                                                              <  2:4546215/101‑1 (MQ=255)
                                   caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCg                                                              >  2:4546565/1‑101 (MQ=255)
                                    aGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAgg                                                                <  2:4546425/98‑1 (MQ=255)
                                      ccAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGttt                                                           >  2:4546570/1‑101 (MQ=255)
                                       cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTa                                                          <  1:4546034/101‑1 (MQ=255)
                                       cAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTa                                                          >  2:4546573/1‑101 (MQ=255)
                                        aGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAc                                                         <  1:4545967/101‑1 (MQ=255)
                                         gtgtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcg                                                        >  1:4546576/1‑101 (MQ=255)
                                           gtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                      <  1:4546135/101‑1 (MQ=255)
                                           gtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgcc                                                      <  2:4546313/101‑1 (MQ=255)
                                             tttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                    >  1:4546580/1‑101 (MQ=255)
                                               tgatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCAt                                                  >  1:4546583/1‑101 (MQ=255)
                                               tgatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCAt                                                  <  1:4545950/101‑1 (MQ=255)
                                                gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAATGCGTTTACGCCGCATc                                                 >  1:4546586/1‑101 (MQ=255)
                                                 atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                <  1:4546288/101‑1 (MQ=255)
                                                 atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                >  1:4546589/1‑101 (MQ=255)
                                                  tgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCg                                               >  2:4546590/1‑101 (MQ=255)
                                                    aGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGc                                                               >  2:4546593/1‑83 (MQ=255)
                                                     ggCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACa                                            <  2:4546233/101‑1 (MQ=255)
                                                        gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACg                                                                                    >  1:4546597/1‑58 (MQ=255)
                                                         aaGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCt                                        >  2:4546599/1‑101 (MQ=255)
                                                           gCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGt                                      >  2:4546604/1‑101 (MQ=255)
                                                            cAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                     <  1:4546242/101‑1 (MQ=255)
                                                            cAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                     >  1:4546605/1‑101 (MQ=255)
                                                             aGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTAc                                    <  1:4546013/101‑1 (MQ=255)
                                                             aGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTAc                                    >  1:4546607/1‑101 (MQ=255)
                                                             aGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTAc                                    <  1:4546019/101‑1 (MQ=255)
                                                             aGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTAc                                    <  1:4546117/101‑1 (MQ=255)
                                                               tAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGt                                  <  1:4546177/101‑1 (MQ=255)
                                                               tAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGt                                  <  2:4546110/101‑1 (MQ=255)
                                                                aGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                 <  1:4546190/101‑1 (MQ=255)
                                                                aGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                 <  2:4546132/101‑1 (MQ=255)
                                                                aGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTg                                 >  2:4546610/1‑101 (MQ=255)
                                                                  cgccgcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACgggg                                                                                 <  1:4546147/51‑1 (MQ=255)
                                                                  cgccgcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGaa                               <  2:4546118/101‑1 (MQ=255)
                                                                   gccgcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAAc                              >  1:4546614/1‑101 (MQ=255)
                                                                     cgcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAg                            >  1:4546619/1‑101 (MQ=255)
                                                                     cgcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAg                            >  2:4546618/1‑101 (MQ=255)
                                                                      gcTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                           >  2:4546622/1‑101 (MQ=255)
                                                                       cTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                          <  1:4546317/101‑1 (MQ=255)
                                                                       cTCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                          >  1:4546624/1‑101 (MQ=255)
                                                                        tCCTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaa                         >  1:4546628/1‑101 (MQ=255)
                                                                         ccTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                        >  1:4546630/1‑101 (MQ=255)
                                                                          cTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                       <  1:4546370/101‑1 (MQ=255)
                                                                          cTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                       >  2:4546635/1‑101 (MQ=255)
                                                                          cTATGCCGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                       >  1:4546634/1‑101 (MQ=255)
                                                                               ccGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCt                  <  2:4546262/101‑1 (MQ=255)
                                                                                cGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                 <  1:4546156/101‑1 (MQ=255)
                                                                                cGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                 >  2:4546641/1‑101 (MQ=255)
                                                                                cGGAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTa                 >  1:4546640/1‑101 (MQ=255)
                                                                                  gAGGCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATg               <  1:4546200/101‑1 (MQ=255)
                                                                                    ggCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGtt             >  1:4546644/1‑101 (MQ=255)
                                                                                    ggCGACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGtt             <  1:4546231/101‑1 (MQ=255)
                                                                                       gACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGc          >  1:4546649/1‑101 (MQ=255)
                                                                                       gACGCTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGc          >  2:4546650/1‑101 (MQ=255)
                                                                                           cTGCGCGCCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCggg      >  1:4546654/1‑101 (MQ=255)
                                                                                            tGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt     <  2:4546036/101‑1 (MQ=255)
                                                                                             gcgcgcCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCg        <  1:4546077/97‑1 (MQ=255)
                                                                                             gcgcgcCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTa    >  2:4546656/1‑101 (MQ=255)
                                                                                               gcgcCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAg  <  2:4546353/101‑1 (MQ=255)
                                                                                               gcgcCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAg  >  1:4546660/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAG  >  NC_000913/2793950‑2794145

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