Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 121,433 A→G V40A (GTA→GCA)  aroP ← aromatic amino acid transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913121,4330AG86.2% 721.7 / 97.8 298V40A (GTA→GCA) aroParomatic amino acid transport protein
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (152/105);  total (173/125)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-02

Only 100 of 299 total aligned reads displayed.

CACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGC  >  NC_000913/121342‑121525
                                                                                           |                                                                                            
caccaTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTAc                                                                                             >  1:4602097/1‑93 (MQ=255)
  ccaTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCgg                                                                                           >  1:6133658/1‑93 (MQ=255)
  ccaTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCgg                                                                                           >  2:6535728/1‑93 (MQ=255)
   caTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGa                                                                                          >  1:5762041/1‑93 (MQ=255)
   caTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGa                                                                                          <  1:657641/93‑1 (MQ=255)
   caTTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGa                                                                                          >  1:3739229/1‑93 (MQ=255)
     tttCATCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAgg                                                                                        <  1:5588847/93‑1 (MQ=255)
     tttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAgg                                                                                        >  2:4220913/1‑93 (MQ=255)
     tttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAgg                                                                                        >  1:2297012/1‑93 (MQ=255)
     tttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAgg                                                                                        >  2:675991/1‑93 (MQ=255)
      ttCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGgc                                                                                       <  1:1001053/93‑1 (MQ=255)
        cACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGgcgc                                                                                     <  2:5043128/93‑1 (MQ=255)
        cACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGgcgc                                                                                     <  1:4712976/93‑1 (MQ=255)
         aCCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCt                                                                                    >  1:3105522/1‑93 (MQ=255)
         aCCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCt                                                                                    >  1:7034577/1‑93 (MQ=255)
         aCCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCt                                                                                    >  1:4842295/1‑93 (MQ=255)
         aCCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCATGGCCCTGCGGACTGTATTGCGGAGGCGCt                                                                                    >  1:2779544/1‑93 (MQ=255)
          cccAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTa                                                                                   >  2:6290943/1‑93 (MQ=255)
                 gACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGa                                                                            >  1:3895810/1‑93 (MQ=255)
                   cGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAAt                                                                          <  2:2822913/93‑1 (MQ=255)
                    gCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATa                                                                         <  1:909603/93‑1 (MQ=255)
                     cATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCGGGAATaa                                                                        <  2:1924415/93‑1 (MQ=255)
                     cATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATaa                                                                        >  1:1343742/1‑93 (MQ=255)
                      aTGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCTGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAAc                                                                       >  1:3705424/1‑93 (MQ=255)
                      aTGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAAc                                                                       <  1:2087198/93‑1 (MQ=255)
                       tGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAAcc                                                                      >  1:2666319/1‑93 (MQ=255)
                        gATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAAccc                                                                     >  1:1662160/1‑93 (MQ=255)
                        gATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAAccc                                                                     <  2:6083681/93‑1 (MQ=255)
                        gATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAAccc                                                                     >  2:2150430/1‑93 (MQ=255)
                          tCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCgg                                                                   >  2:5690933/1‑93 (MQ=255)
                            aGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTc                                                                 <  2:2387877/93‑1 (MQ=255)
                             gAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTcc                                                                >  2:750144/1‑93 (MQ=255)
                             gAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTcc                                                                >  2:7269364/1‑93 (MQ=255)
                              aaaGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTccc                                                               >  2:4850724/1‑93 (MQ=255)
                              aaaGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTccc                                                               >  1:5084241/1‑93 (MQ=255)
                                aGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTa                                                             <  2:5125177/93‑1 (MQ=255)
                                aGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTa                                                             >  2:3163147/1‑93 (MQ=255)
                                 ggCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTAt                                                            <  1:2119554/93‑1 (MQ=255)
                                 ggCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTAt                                                            >  1:72598/1‑93 (MQ=255)
                                 ggCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTAt                                                            <  2:4581197/93‑1 (MQ=255)
                                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTATCCAGGAATAACCCGGTCCCTATc                                                           >  1:2835025/1‑93 (MQ=255)
                                  gCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATc                                                           <  1:4778217/93‑1 (MQ=255)
                                     aTAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATcgcg                                                        >  1:1777943/1‑93 (MQ=255)
                                      tAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCgcgc                                                       <  2:5476932/93‑1 (MQ=255)
                                       aaaaCCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGcc                                                      <  2:3499612/93‑1 (MQ=255)
                                        aaaCCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCa                                                     >  1:4141938/1‑93 (MQ=255)
                                        aaaCCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCa                                                     >  2:6792591/1‑93 (MQ=255)
                                         aaCCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCAc                                                    >  2:7142483/1‑93 (MQ=255)
                                           ccAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCAccc                                                  >  1:7356269/1‑93 (MQ=255)
                                            cAGCAATGGCGTAACCCAGGATAATCCTTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCa                                                 >  2:1730977/1‑93 (MQ=255)
                                            cAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCa                                                 <  1:1583559/93‑1 (MQ=255)
                                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAg                                                >  1:4162667/1‑93 (MQ=255)
                                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAg                                                <  2:6844494/93‑1 (MQ=255)
                                             aGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAg                                                <  2:2898037/93‑1 (MQ=255)
                                              gCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAgc                                               <  2:4018606/93‑1 (MQ=255)
                                               cAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAgcg                                              <  1:2915242/93‑1 (MQ=255)
                                                aaTGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAgcgc                                             >  1:7134707/1‑93 (MQ=255)
                                                aaTGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAgcgc                                             <  2:5761979/93‑1 (MQ=255)
                                                  tGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGa                                           >  2:5254302/1‑93 (MQ=255)
                                                  tGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGa                                           >  2:4128115/1‑93 (MQ=255)
                                                    gCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATa                                         <  2:6677501/93‑1 (MQ=255)
                                                     cGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATaa                                        <  2:2580563/93‑1 (MQ=255)
                                                     cGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATaa                                        <  1:5260432/93‑1 (MQ=255)
                                                      gTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAg                                       >  1:5906270/1‑93 (MQ=255)
                                                       tAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGc                                      <  2:1612971/93‑1 (MQ=255)
                                                           ccAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGaa                                  >  1:5449385/1‑93 (MQ=255)
                                                            cAGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAat                                 <  1:2520100/93‑1 (MQ=255)
                                                             aGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAata                                >  1:6465556/1‑93 (MQ=255)
                                                             aGGATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAata                                <  2:3782642/93‑1 (MQ=255)
                                                              ggATAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAatat                               >  2:3091718/1‑93 (MQ=255)
                                                                aTAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATgg                             <  2:796751/93‑1 (MQ=255)
                                                                 tAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGc                            <  1:363426/93‑1 (MQ=255)
                                                                 tAATCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGc                            >  2:7044051/1‑93 (MQ=255)
                                                                  aaTCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCg                           <  2:6029815/93‑1 (MQ=255)
                                                                   aTCCCTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCgg                          >  2:2939017/1‑93 (MQ=255)
                                                                     cccTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGtt                        <  2:2679339/93‑1 (MQ=255)
                                                                      ccTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGttt                       <  2:4291602/93‑1 (MQ=255)
                                                                      ccTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGttt                       <  1:1130509/93‑1 (MQ=255)
                                                                       cTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGtttt                      <  1:1723004/93‑1 (MQ=255)
                                                                       cTGGCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGtttt                      >  2:620535/1‑93 (MQ=255)
                                                                          gCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTaa                   >  2:1303298/1‑93 (MQ=255)
                                                                          gCCCTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTaa                   <  2:3682099/93‑1 (MQ=255)
                                                                             cTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGc                >  2:3333995/1‑93 (MQ=255)
                                                                             cTGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGc                <  1:4385959/93‑1 (MQ=255)
                                                                             cTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGc                >  1:663672/1‑93 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTTAgggc               >  1:5313432/1‑91 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGcc               >  2:4466872/1‑93 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGcc               <  1:3246786/93‑1 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGcc               <  1:3443232/93‑1 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGcc               <  2:6585707/93‑1 (MQ=255)
                                                                              tGCGGACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGcc               <  2:4756834/93‑1 (MQ=255)
                                                                               gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCcg              <  2:3131703/93‑1 (MQ=255)
                                                                               gCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCcg              >  1:3418635/1‑93 (MQ=255)
                                                                                 ggACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCcgcg            >  2:1845739/1‑93 (MQ=255)
                                                                                  gACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCgcgc           >  2:3228439/1‑93 (MQ=255)
                                                                                  gACTGTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCgcgc           >  1:1854753/1‑93 (MQ=255)
                                                                                      gTATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTa       >  1:139777/1‑93 (MQ=255)
                                                                                       tATTGCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAg      <  1:1017624/93‑1 (MQ=255)
                                                                                        aTTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAgc     <  1:3763026/93‑1 (MQ=255)
                                                                                           aCGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAgctgc  <  2:6492499/93‑1 (MQ=255)
                                                                                           |                                                                                            
CACCATTTCACCCAGCTGACGCATGATCAGAAAGGCGATAAAACCAGCAATGGCGTAACCCAGGATAATCCCTGGCCCTGCGGACTGTATTACGGAGGCGCTACCCAGGAATAACCCGGTCCCTATCGCGCCACCCAGCGCGATAAGCTGAATATGGCGGTTTTTAAGGCCGCGCTTTAGCTGC  >  NC_000913/121342‑121525

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