Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,909,645 T→C D551G (GAT→GGT)  glmS ← L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,909,6450TC84.8% 728.1 / 118.5 289D551G (GAT→GGT) glmSL‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (25/19);  new base C (137/108);  total (162/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 296 total aligned reads displayed.

GCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCA  >  NC_000913/3909553‑3909737
                                                                                            |                                                                                            
gCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCAt                                                                                              >  2:5228090/1‑93 (MQ=255)
  aGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCAccc                                                                                            >  2:5086424/1‑93 (MQ=255)
  aGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATcc                                                                                            >  1:4997533/1‑93 (MQ=255)
    cGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTg                                                                                          <  1:6278099/93‑1 (MQ=255)
     ggAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGa                                                                                         >  2:2246341/1‑93 (MQ=255)
     ggAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGa                                                                                         >  1:1672499/1‑93 (MQ=255)
     ggAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGa                                                                                         <  1:1768499/93‑1 (MQ=255)
     ggAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGa                                                                                         <  2:1299799/93‑1 (MQ=255)
      gAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGAt                                                                                        >  1:4940083/1‑93 (MQ=255)
       aaCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATc                                                                                       <  2:1065458/93‑1 (MQ=255)
       aaCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATc                                                                                       >  2:420949/1‑93 (MQ=255)
       aaCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATc                                                                                       <  1:2380300/93‑1 (MQ=255)
        aCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCg                                                                                      <  1:5075716/93‑1 (MQ=255)
        aCGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCg                                                                                      >  1:2717902/1‑93 (MQ=255)
         cGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCgg                                                                                     <  1:2273780/93‑1 (MQ=255)
            tgtAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGa                                                                                  >  2:4169576/1‑93 (MQ=255)
               agaagaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGa                                                                               >  2:6047433/1‑93 (MQ=255)
                  agaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACAt                                                                            >  1:370684/1‑93 (MQ=255)
                   gaTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGCAGACATa                                                                           >  1:4515015/1‑93 (MQ=255)
                    aTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATAc                                                                          <  2:3229977/93‑1 (MQ=255)
                    aTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATAc                                                                          >  2:1997701/1‑93 (MQ=255)
                    aTCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATAc                                                                          <  2:5947479/93‑1 (MQ=255)
                       ggTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAAc                                                                       >  2:1661831/1‑93 (MQ=255)
                        gTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACt                                                                      >  1:3375350/1‑93 (MQ=255)
                        gTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACt                                                                      >  2:4572634/1‑93 (MQ=255)
                           cAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGAc                                                                   <  2:1171418/93‑1 (MQ=255)
                            aaTCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGAcc                                                                  >  2:3660595/1‑93 (MQ=255)
                            aaTCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGAcc                                                                  <  2:5633850/93‑1 (MQ=255)
                             aTCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCg                                                                 >  1:3567882/1‑93 (MQ=255)
                             aTCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAAATGACCg                                                                 >  2:464639/1‑93 (MQ=255)
                               cACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGcc                                                               <  2:304007/93‑1 (MQ=255)
                                 ccTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAc                                                             >  1:1121838/1‑93 (MQ=255)
                                 ccTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAc                                                             <  1:1630763/93‑1 (MQ=255)
                                  ctctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcg                                                            >  1:2543134/1‑93 (MQ=255)
                                  ctctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcg                                                            >  2:4910979/1‑93 (MQ=255)
                                  ctctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcg                                                            <  1:1380227/93‑1 (MQ=255)
                                   tctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgc                                                           >  1:959388/1‑93 (MQ=255)
                                   tctTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgc                                                           >  1:6190837/1‑93 (MQ=255)
                                     ttCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgc                                                         <  2:607036/93‑1 (MQ=255)
                                     ttCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgc                                                         >  2:1836733/1‑93 (MQ=255)
                                     ttCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgc                                                         >  1:832460/1‑93 (MQ=255)
                                      tCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgcg                                                        >  1:5098790/1‑93 (MQ=255)
                                       ccACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgcgc                                                       <  1:5243003/93‑1 (MQ=255)
                                        cacaTGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgcgcg                                                      >  1:4252552/1‑93 (MQ=255)
                                        cacaTGCGACATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCAcgcgcgcg                                                      >  1:6501167/1‑93 (MQ=255)
                                             gCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAActt                                                 <  1:6437153/93‑1 (MQ=255)
                                               ggCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAActtct                                               >  2:1372254/1‑93 (MQ=255)
                                                gCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAActtctt                                              >  2:4819913/1‑93 (MQ=255)
                                                 cATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACttcttc                                             <  2:4134073/93‑1 (MQ=255)
                                                 cATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACttcttc                                             >  1:2057006/1‑93 (MQ=255)
                                                   tctcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCaa                                           >  2:2424145/1‑93 (MQ=255)
                                                   tctcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCaa                                           >  1:4667824/1‑93 (MQ=255)
                                                    ctcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGTGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAAt                                          <  2:3567006/93‑1 (MQ=255)
                                                     tcGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATg                                         <  2:3753002/93‑1 (MQ=255)
                                                      cGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGt                                        <  1:4078376/93‑1 (MQ=255)
                                                      cGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGt                                        >  1:5383084/1‑93 (MQ=255)
                                                      cGATGATGTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGt                                        >  1:1056744/1‑93 (MQ=255)
                                                          gatgTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGa                                    >  2:4421261/1‑93 (MQ=255)
                                                          gatgTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGa                                    <  1:1203842/93‑1 (MQ=255)
                                                          gatgTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGa                                    <  2:2816373/93‑1 (MQ=255)
                                                           atgTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGAt                                   >  2:2881298/1‑93 (MQ=255)
                                                           atgTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGAt                                   <  1:2846129/93‑1 (MQ=255)
                                                           atgTGCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGAt                                   >  2:5723063/1‑93 (MQ=255)
                                                            tgtgCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCACGCGAACTTCTTCAATGTTGGAtt                                  >  1:3262773/1‑93 (MQ=255)
                                                             gtgCATGTTATAGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGAttt                                 <  2:6254025/93‑1 (MQ=255)
                                                              tgCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTc                                >  1:5849458/1‑93 (MQ=255)
                                                              tgCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTc                                >  1:4342619/1‑93 (MQ=255)
                                                               gCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCa                               <  1:5788701/93‑1 (MQ=255)
                                                               gCATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCa                               >  1:2903859/1‑93 (MQ=255)
                                                                cATGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAg                              >  1:200340/1‑93 (MQ=255)
                                                                  tGTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGtt                            >  2:5574803/1‑93 (MQ=255)
                                                                   gTTATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGttt                           <  1:4463006/93‑1 (MQ=255)
                                                                     tATCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGttttt                         <  2:3545662/93‑1 (MQ=255)
                                                                      aTCGCTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTc                        <  2:6346991/93‑1 (MQ=255)
                                                                        cGCTACTTACAAAGCCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCa                      <  2:3593116/93‑1 (MQ=255)
                                                                        cGCTACTTACAAAACCCGCGCCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCa                      <  1:3045574/93‑1 (MQ=255)
                                                                         gCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAg                     >  2:5427358/1‑93 (MQ=255)
                                                                          cTACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGc                    <  1:5477525/93‑1 (MQ=255)
                                                                           tACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCa                   >  2:821449/1‑93 (MQ=255)
                                                                           tACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCa                   >  2:1439053/1‑93 (MQ=255)
                                                                           tACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCa                   >  1:3666763/1‑93 (MQ=255)
                                                                           tACTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCa                   >  1:5947082/1‑93 (MQ=255)
                                                                            aCTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCaa                  <  2:4342582/93‑1 (MQ=255)
                                                                             cTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAAt                 >  1:2123148/1‑93 (MQ=255)
                                                                             cTTACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAAt                 <  2:653503/93‑1 (MQ=255)
                                                                              ttACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAAtt                >  2:1547139/1‑93 (MQ=255)
                                                                              ttACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAAtt                <  2:390265/93‑1 (MQ=255)
                                                                              ttACAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAAtt                >  2:26477/1‑93 (MQ=255)
                                                                                aCAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCg              >  2:3445085/1‑93 (MQ=255)
                                                                                aCAAAACCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCg              <  2:5855850/93‑1 (MQ=255)
                                                                                  aaaaCCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgtt            <  2:5332159/93‑1 (MQ=255)
                                                                                   aaaCCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttg           >  1:4884405/1‑93 (MQ=255)
                                                                                    accccGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgt          <  2:248961/91‑1 (MQ=255)
                                                                                    aaCCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgt          >  1:6602251/1‑93 (MQ=255)
                                                                                    aaCCCGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCgttgt          <  1:3930930/93‑1 (MQ=255)
                                                                                      cccGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTc        <  1:6034380/93‑1 (MQ=255)
                                                                                      cccGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTc        <  2:1421191/93‑1 (MQ=255)
                                                                                       ccGCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCg       >  1:726297/1‑93 (MQ=255)
                                                                                         gCACCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGt     <  1:5742266/93‑1 (MQ=255)
                                                                                            tCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGca  >  2:6446213/1‑93 (MQ=255)
                                                                                            |                                                                                            
GCAGCGGAACGGTGTAGAAGATCGGTGCAATCACCTCTTCCACATGCGGCATCTCGATGATGTGCATGTTATCGCTACTTACAAAACCCGCATCCTGATCGGCGAAGACATACAACTGACCGCCACGCGCGCGAACTTCTTCAATGTTGGATTTCAGTTTTTCCAGCAATTCGTTGTTCGGTGCA  >  NC_000913/3909553‑3909737

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