Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 754,103 G→T 15.7% F122L (TTC→TTA)  gltA ← citrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913754,1030GT15.7% 171.2 / 16.9 61F122L (TTC→TTA) gltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base T (4/6);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGCATTTTCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATC  >  NC_000913/753957‑754250
                                                                                                                                                  |                                                                                                                                                   
gTCGGCATTTTCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATgg                                                                                                                                                     >  3:207692/1‑147 (MQ=255)
      aTTTTCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGAGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGAc                                                                                                                                              >  5:169632/1‑148 (MQ=255)
       ttttCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACg                                                                                                                                             >  1:3772/1‑148 (MQ=255)
        tttCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACgg                                                                                                                                            >  6:124981/1‑148 (MQ=255)
             aCAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGCCGGAAAGCATGGAACAGAAGGGTaaa                                                                                                                                       >  6:78411/1‑147 (MQ=255)
                gcagGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTg                                                                                                                                    >  2:127021/1‑148 (MQ=255)
                  agGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACAGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCt                                                                                                                                  >  6:274031/1‑148 (MQ=255)
                     cGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGt                                                                                                                               <  2:2716/148‑1 (MQ=255)
                       gAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTg                                                                                                                              >  5:263034/1‑147 (MQ=255)
                         aCGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGAt                                                                                                                           >  5:242094/1‑148 (MQ=255)
                            cgGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        <  7:209683/148‑1 (MQ=255)
                            cgGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        <  3:205926/148‑1 (MQ=255)
                            cgGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        <  7:78013/148‑1 (MQ=255)
                            cgGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        >  2:240095/1‑148 (MQ=255)
                                             cGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCg                                                                                                       >  7:11752/1‑148 (MQ=255)
                                               aGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCg                                                                                                       >  3:152548/1‑146 (MQ=255)
                                               aGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTa                                                                                                     <  5:242078/148‑1 (MQ=255)
                                                 gATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGt                                                                                                                               <  4:267684/120‑1 (MQ=255)
                                                 gATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGt                                                                                                                               >  3:267684/1‑120 (MQ=255)
                                                 gATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGt                                                                                                                               >  3:271994/1‑120 (MQ=255)
                                                 gATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGt                                                                                                                               <  4:271994/120‑1 (MQ=255)
                                                       tAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGATCAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTaaa                                                                                             >  3:296649/1‑148 (MQ=255)
                                                       tAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTaaa                                                                                             >  2:35260/1‑148 (MQ=255)
                                                        aaCATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTaaa                                                                                             <  8:11752/147‑1 (MQ=255)
                                                         aCATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAAt                                                                                            <  5:124981/147‑1 (MQ=255)
                                                          cATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTc                                                                                          >  3:161905/1‑148 (MQ=255)
                                                             ccAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        <  5:229018/115‑1 (MQ=255)
                                                             ccAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCAt                                                                                                                        >  6:229018/1‑115 (MQ=255)
                                                              cAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCGCACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTc                                                                                       <  2:3122/147‑1 (MQ=255)
                                                                    gTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGt                                                                                <  4:207692/148‑1 (MQ=255)
                                                                     tCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGATCAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGtt                                                                               <  4:296649/148‑1 (MQ=255)
                                                                     tCGGGATAGAAGGCCGCCACGTCGCCGGTATTACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCCTGGAACAGACGGGTAATCTGCTCGTGGTGCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGtt                                                                               <  1:211833/148‑1 (MQ=255)
                                                                       gTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTcc                                                                             <  6:242094/148‑1 (MQ=255)
                                                                            aGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAg                                                                         <  4:152548/147‑1 (MQ=255)
                                                                            aGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGt                                                                        >  1:92236/1‑148 (MQ=255)
                                                                             gAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAGTGGGGGATGGTGTTCCTGAggc                                                                       >  4:231381/1‑146 (MQ=255)
                                                                              aaCGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCg                                                                      <  6:263034/148‑1 (MQ=255)
                                                                                 gccgccAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGGCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGtt                                                                   <  1:240095/148‑1 (MQ=255)
                                                                                     ccAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTc                                                               <  3:180659/148‑1 (MQ=255)
                                                                                         cgcgcCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCa                                                           >  8:250429/1‑148 (MQ=255)
                                                                                             ccGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTcc                                                                             <  1:279431/126‑1 (MQ=255)
                                                                                             ccGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTcc                                                                             >  2:279431/1‑126 (MQ=255)
                                                                                                     aCCACACATGACTGCGATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGt                                                                                                                     >  3:264230/1‑78 (MQ=255)
                                                                                                     aCCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGt                                                                                                                     <  4:264230/78‑1 (MQ=255)
                                                                                                     aCCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATg                                               >  4:106646/1‑148 (MQ=255)
                                                                                                     aCCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAAACTGCTCGTGGATCATGGTATGATGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATg                                               >  4:62817/1‑148 (MQ=255)
                                                                                                     aCCACACATGACTGCCATTGGATGAGAGACGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGACGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACAGTTcc                                                                             >  7:39000/1‑118 (MQ=255)
                                                                                                          acaTGACTGCCATTGGATGCGAGTCGCCACTGAAAGCATTTAACAGACGGGTAATCTGCTCGTGGATCATGGTATCACGGGACACCGTAGTTTTAAATTAGTCATACTGTTCCAGAGTCGGTATATCACCATTCAGCAGGATGTAAc                                           >  4:89406/1‑147 (MQ=255)
                                                                                                                cttcaaTTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGTTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTGTTTCACCATTCAGCAGGATGTAACAAACTTc                                    <  5:78411/143‑1 (MQ=255)
                                                                                                                   ccATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTcc                                   <  3:106646/146‑1 (MQ=255)
                                                                                                                     cattgATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGtt                                                                   >  6:105679/5‑112 (MQ=255)
                                                                                                                     cattgATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGtt                                                                   <  5:105679/108‑1 (MQ=255)
                                                                                                                       tGGATGCGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAg                             >  4:38515/1‑148 (MQ=255)
                                                                                                                             cGAGTCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGa                        <  3:38515/147‑1 (MQ=255)
                                                                                                                             cGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGATGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGa                        <  3:62817/147‑1 (MQ=255)
                                                                                                                               aGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATc                     <  1:127021/148‑1 (MQ=255)
                                                                                                                                 tCGCGACGGAAAGCATGTAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCgg                   <  5:274031/148‑1 (MQ=255)
                                                                                                                                 tCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCgg                   <  3:36631/148‑1 (MQ=255)
                                                                                                                                   gcgACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGGCATACTGTTCCTGAGTCGGttttt                                                                <  4:193178/101‑1 (MQ=255)
                                                                                                                                   gcgACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGACATACTGTTCCTGAGTCGGttttt                                                                >  3:193178/1‑101 (MQ=255)
                                                                                                                                       cGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCTCCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAAGTTCCAGGTAGTTAGAATCGTTCGCCa             <  2:51137/148‑1 (MQ=255)
                                                                                                                                         gAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGc           >  4:297246/1‑148 (MQ=255)
                                                                                                                                                 ggAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcgat   >  4:301445/1‑148 (MQ=255)
                                                                                                                                                  gAACCGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTCAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTTGAATCGGTCGCCAGCTgatcga    >  8:19567/1‑146 (MQ=255)
                                                                                                                                                  gAACAGACTGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTgatcgatc  <  2:3772/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTCGGCATTTTCGACAGCAGGCGGAACGCGGCAATTTCACGGTGACGAGGATTGTTAACATCCAGCGAGTCGTGATAGAACGCCGCCAGCGCGCCGGTAATACCACACATGACTGCCATTGGATGCGAGTCGCGACGGAAAGCATGGAACAGACGGGTAATCTGCTCGTGGATCATGGTATGACGGGTCACCGTAGTTTTAAATTCGTCATACTGTTCCTGAGTCGGTTTTTCACCATTCAGCAGGATGTAACAAACTTCCAGGTAGTTAGAATCGGTCGCCAGCTGATCGATC  >  NC_000913/753957‑754250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑