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,9610AG98.2% 521.7 / ‑1.0 168S324P (TCC→CCC) speAbiosynthetic arginine decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (101/63);  total (102/66)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 175 total aligned reads displayed.

GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050
                                                                                          |                                                                                         
gCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCt                                                                                         <  2:2998272/93‑1 (MQ=255)
 cGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTg                                                                                        >  2:4056303/1‑93 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGc                                                                                       <  1:1901856/93‑1 (MQ=255)
  ggCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGc                                                                                       >  2:4137292/1‑93 (MQ=255)
    cAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGa                                                                                     >  1:3890484/1‑93 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  1:1595759/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  1:985086/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    <  1:3685425/93‑1 (MQ=255)
     aGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGaa                                                                                    >  1:46622/1‑93 (MQ=255)
      gACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAAc                                                                                   <  2:3801751/93‑1 (MQ=255)
         cGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAg                                                                                <  1:144559/93‑1 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               <  2:2060355/93‑1 (MQ=255)
          gTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGt                                                                               <  1:2072395/93‑1 (MQ=255)
           ttttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                              >  2:1586730/1‑93 (MQ=255)
           ttttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTa                                                                              <  1:3083616/93‑1 (MQ=255)
            tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTAc                                                                             >  1:769541/1‑93 (MQ=255)
            tttCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTAc                                                                             <  2:1589572/93‑1 (MQ=255)
             ttcttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGTCTGCGAACGAGTAcc                                                                            <  2:2301372/93‑1 (MQ=255)
              tcttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCt                                                                           >  1:4208028/1‑93 (MQ=255)
              tattcACACGAGTCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCt                                                                           <  2:4310629/91‑1 (MQ=255)
               cttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCtt                                                                          <  2:1142066/93‑1 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  1:3997867/1‑93 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         <  2:1849465/93‑1 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  1:3687340/1‑93 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  1:2734742/1‑93 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         <  2:2149813/93‑1 (MQ=255)
                ttcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTc                                                                         >  2:4556250/1‑93 (MQ=255)
                 tcACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCa                                                                        <  2:1137942/93‑1 (MQ=255)
                    cacGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATaa                                                                     <  1:2250957/93‑1 (MQ=255)
                         aTCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCAc                                                                <  2:77311/93‑1 (MQ=255)
                          tCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACg                                                               >  2:1242923/1‑93 (MQ=255)
                           cGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGc                                                              <  1:3373406/93‑1 (MQ=255)
                             ccAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGccc                                                            >  1:4354027/1‑93 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           >  1:3340226/1‑93 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           <  2:2957162/93‑1 (MQ=255)
                              cAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCa                                                           <  1:2327587/93‑1 (MQ=255)
                               caatgCCCAGATAATGTTGTTGGCGTAATCATTTTGGCCGTAGTTTTCCCAACAGAGCGGCTGCGTACGAGAACCTTCTTAATCCACGCCCAg                                                          >  2:3186505/5‑93 (MQ=255)
                                 tCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAc                                                        <  1:2749956/93‑1 (MQ=255)
                                    cccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgc                                                     >  1:3069606/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    >  2:4174810/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    >  2:141124/1‑93 (MQ=255)
                                     ccAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCCGCGAACGAGTACCTTCATAATCCACGCCCAGAccgcc                                                    <  2:2346419/93‑1 (MQ=255)
                                     acAGATAATGTTGTTGGCGTATTCATTTAGGGCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGGACGTTCATAATCCACGCCCAGAccgcc                                                    >  2:44318/2‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  >  2:2501721/1‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  >  2:2473522/1‑93 (MQ=255)
                                       aGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGAccgccgc                                                  >  1:3483371/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCCCCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               <  2:894577/93‑1 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               <  2:1493606/93‑1 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               >  1:1013251/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               >  2:2864727/1‑93 (MQ=255)
                                          tAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGa                                               <  1:2054311/93‑1 (MQ=255)
                                            aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACg                                             >  1:3354757/1‑93 (MQ=255)
                                            aTGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACg                                             >  2:4501328/1‑93 (MQ=255)
                                             tgttgtTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGt                                            <  1:432075/93‑1 (MQ=255)
                                              gttgttGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTc                                           <  1:1846952/93‑1 (MQ=255)
                                                tgttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTcca                                         >  1:1015772/1‑91 (MQ=255)
                                                 gttgGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGaa                                        >  1:3463181/1‑93 (MQ=255)
                                                    ggCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCa                                     <  2:1179021/93‑1 (MQ=255)
                                                      cGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACt                                   >  1:1720285/1‑93 (MQ=255)
                                                       gTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg                                  >  2:1021038/1‑93 (MQ=255)
                                                       gTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTg                                  >  2:1768281/1‑93 (MQ=255)
                                                         aTTCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGaa                                <  2:4074259/93‑1 (MQ=255)
                                                          ttCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAat                               >  1:1162375/1‑93 (MQ=255)
                                                           tCATTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAata                              >  1:2550500/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            >  1:943428/1‑93 (MQ=255)
                                                             aTTGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATAtt                            >  2:4127709/1‑93 (MQ=255)
                                                              ttGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTg                           <  2:2988386/93‑1 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          >  2:514521/1‑93 (MQ=255)
                                                               tGAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGa                          <  1:2869158/93‑1 (MQ=255)
                                                                gAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGAc                         <  2:1823074/93‑1 (MQ=255)
                                                                gAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGAc                         >  1:1487608/1‑93 (MQ=255)
                                                                gAGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGAc                         <  2:2227614/93‑1 (MQ=255)
                                                                 aGGCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACg                        >  2:1134744/1‑93 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGc                       >  2:1318139/1‑93 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGc                       <  1:1759607/93‑1 (MQ=255)
                                                                  ggCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCAATACCGCCGCCGACCGCGAAGCACTGAATATGGACGc                       >  1:999214/1‑93 (MQ=255)
                                                                   gCCGTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGcc                      <  2:3321892/93‑1 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                   >  2:2740598/1‑93 (MQ=255)
                                                                      gTAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAg                   <  2:4023437/93‑1 (MQ=255)
                                                                       tAGTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGc                  <  2:4500812/93‑1 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGCCGTCGAAGCACTGACTATTGACGACCCGCtg                >  1:4135088/1‑92 (MQ=255)
                                                                         gTTCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCtt                <  1:2337054/93‑1 (MQ=255)
                                                                          ttCACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtg               >  1:528437/1‑93 (MQ=255)
                                                                            cACCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTtgtg             >  1:4392286/1‑93 (MQ=255)
                                                                             aCCGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGc            >  2:1676969/1‑93 (MQ=255)
                                                                             aCCGAACAGTCGGGCTGCGAACGAGTACATTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGc            >  1:2748174/1‑93 (MQ=255)
                                                                              ccGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCa           <  2:3648930/93‑1 (MQ=255)
                                                                               cGAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAg          <  2:1943121/93‑1 (MQ=255)
                                                                                gAACAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGt         <  1:3221649/93‑1 (MQ=255)
                                                                                 aaCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCCGtt        >  1:4058799/1‑93 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTc       <  1:2200244/93‑1 (MQ=255)
                                                                                  aCAGTCGGGCTGCGAACGAGTACATTCATAATCCACGCCCAGACAGCCGCCGACGTCGAAGCACTGAATATTGACGGCCAGCTTGTGCAGTTc       >  2:4729685/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTacac    >  2:1072226/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    >  1:45524/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    >  1:1732092/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    >  2:1460719/1‑93 (MQ=255)
                                                                                     gTCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcac    >  1:1217834/1‑93 (MQ=255)
                                                                                      tCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcaca   <  2:3355843/93‑1 (MQ=255)
                                                                                      tCGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCcaca   <  1:2916305/93‑1 (MQ=255)
                                                                                       cGGGCTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAt  >  2:565512/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCGGCAGACCGTTTTCTTCACACGCATCGCCAATCGCCCAGATAATGTTGTTGGCGTATTCATTGAGGCCGTAGTTCACCGAACAGTCGGACTGCGAACGAGTACCTTCATAATCCACGCCCAGACCGCCGCCGACGTCGAAGCACTGAATATTGACGCCCAGCTTGTGCAGTTCCACAT  >  NC_000913/3082871‑3083050

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