Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,082,961 A→G S324P (TCC→CCC)  speA ← biosynthetic arginine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,9610AG99.2% 797.8 / ‑3.0 256S324P (TCC→CCC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (128/125);  total (130/126)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 269 total aligned reads displayed.

GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050
                                                                                          |                                                                                         
gCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCt                                                                                         <  2:2463954/93‑1 (MQ=255)
gCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCt                                                                                         <  1:4729798/93‑1 (MQ=255)
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                        <  1:6572195/93‑1 (MQ=255)
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                        <  1:500571/93‑1 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       <  1:4571746/93‑1 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       >  1:2340460/1‑93 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       >  1:3011114/1‑93 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGc                                                                                       >  1:6263931/1‑93 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGc                                                                                       >  1:5576397/1‑93 (MQ=255)
    cAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGa                                                                                     >  1:6384590/1‑93 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  1:3160216/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    >  2:6109649/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   <  2:1683/93‑1 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   >  1:5247920/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   <  2:4586835/93‑1 (MQ=255)
       aCCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACg                                                                                  >  2:1085368/1‑93 (MQ=255)
        ccGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGa                                                                                 >  1:6106090/1‑93 (MQ=255)
         cGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAg                                                                                <  1:2272418/93‑1 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               >  1:5801585/1‑93 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               >  1:242876/1‑93 (MQ=255)
           ttttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                              >  1:1873100/1‑93 (MQ=255)
           ttttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                              >  2:1793871/1‑93 (MQ=255)
              tcttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCt                                                                           <  2:5881380/93‑1 (MQ=255)
               cttcACACGCAGCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCtt                                                                          <  2:3351293/93‑1 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  2:2522290/1‑93 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  1:4368191/1‑93 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  2:3175433/1‑93 (MQ=255)
                 tcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGGACCTTCa                                                                        >  1:2391215/1‑93 (MQ=255)
                  cacaCGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCAt                                                                       >  2:3831764/1‑93 (MQ=255)
                  cacaCGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCAt                                                                       <  2:2720948/93‑1 (MQ=255)
                   acacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATa                                                                      >  1:1380948/1‑93 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     <  1:1538586/93‑1 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     <  2:2418735/93‑1 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     <  2:694091/93‑1 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     <  1:6646138/93‑1 (MQ=255)
                     acGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAAt                                                                    >  1:5639089/1‑93 (MQ=255)
                       gCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATcc                                                                  <  2:6589092/93‑1 (MQ=255)
                       gCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATcc                                                                  >  1:3985583/1‑93 (MQ=255)
                        cATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCa                                                                 <  2:6156670/93‑1 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                               >  1:4436425/1‑93 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                               >  1:6010513/1‑93 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCATATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                               <  1:907421/93‑1 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           >  2:66291/1‑93 (MQ=255)
                              cAATCACCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           >  2:1146069/1‑93 (MQ=255)
                               aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAg                                                          >  2:2555501/1‑93 (MQ=255)
                               aaTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAg                                                          >  2:1228286/1‑93 (MQ=255)
                                aTCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGa                                                         <  1:3872390/93‑1 (MQ=255)
                                    cccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgc                                                     >  2:5757999/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    <  1:4949203/93‑1 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    >  2:831210/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    >  2:4007939/1‑93 (MQ=255)
                                      cAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCAAGCCCAGAccgccg                                                   >  2:5309291/1‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  >  2:1840666/1‑93 (MQ=255)
                                         aTAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgccg                                                >  2:2699593/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               <  1:3743810/93‑1 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               >  2:4814393/1‑93 (MQ=255)
                                           aaTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGAc                                              >  1:4773513/1‑93 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGt                                            <  1:1186506/93‑1 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGt                                            <  1:2731318/93‑1 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCCCCGCCGACGt                                            <  2:1543378/93‑1 (MQ=255)
                                              gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTc                                           <  1:4169583/93‑1 (MQ=255)
                                                tgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCTAGACCGCCGCCGACGTCGa                                         <  1:6472729/93‑1 (MQ=255)
                                                tgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa                                         >  1:1308268/1‑93 (MQ=255)
                                                tgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGa                                         >  2:3424232/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                      >  1:1064738/1‑93 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                      <  1:2840566/93‑1 (MQ=255)
                                                   tgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGc                                      <  1:5356373/93‑1 (MQ=255)
                                                        tATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGa                                 >  1:5763387/1‑93 (MQ=255)
                                                          ttCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAat                               >  2:1022916/1‑93 (MQ=255)
                                                          ttCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAat                               <  1:3260632/93‑1 (MQ=255)
                                                          ttCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAat                               <  2:2138385/93‑1 (MQ=255)
                                                           tCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAata                              <  2:3483688/93‑1 (MQ=255)
                                                           tCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAata                              <  1:2076897/93‑1 (MQ=255)
                                                            cATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAatat                             <  2:5028305/93‑1 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            <  1:5066202/93‑1 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            <  2:3768959/93‑1 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                           >  1:2694876/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                           >  1:732177/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                           >  1:1692730/1‑93 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          <  2:4989526/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          <  1:4891657/93‑1 (MQ=255)
                                                                taGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGAc                         <  1:5913086/92‑1 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGcc                      >  2:745882/1‑93 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGGCGTCCAg                   >  1:5203987/1‑93 (MQ=255)
                                                                       tAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGc                  >  2:5539083/1‑93 (MQ=255)
                                                                        aGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCt                 <  2:4917398/93‑1 (MQ=255)
                                                                         gTTCACCGAACAGTCTGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                <  1:2897454/93‑1 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                <  2:410842/93‑1 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                <  2:2908999/93‑1 (MQ=255)
                                                                           tCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgt              <  2:3045633/93‑1 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAGGCACTGAATATTGACGCCCAGCTtggg             >  2:5935774/1‑91 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgtg             >  2:3714636/1‑93 (MQ=255)
                                                                               cGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAg          >  1:3405547/1‑93 (MQ=255)
                                                                                gAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGt         >  2:4257871/1‑93 (MQ=255)
                                                                                gAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGt         >  1:1926360/1‑93 (MQ=255)
                                                                                gAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGt         >  2:1889296/1‑93 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTc       >  2:1973665/1‑93 (MQ=255)
                                                                                   cAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTcc      <  1:2565493/93‑1 (MQ=255)
                                                                                      tCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcaca   <  2:232778/93‑1 (MQ=255)
                                                                                       cGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCTCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAt  <  2:5594875/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050

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