Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,910,662 A→T 100% D288E (GAT→GAA pstS ← phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,910,6620AT100.0% 376.9 / ‑6.7 116D288E (GAT→GAApstSphosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (63/52);  minor base C (1/0);  total (64/52)

Only 100 of 116 total aligned reads displayed.

GCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910566‑3910758
                                                                                                |                                                                                                
gCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCt                                                                                               >  2:745661/1‑100 (MQ=255)
 cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  2:2992878/1‑100 (MQ=255)
 cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  2:1380371/1‑100 (MQ=255)
 cTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCAAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCtt                                                                                              >  1:1625955/1‑100 (MQ=255)
  tGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCttt                                                                                             <  2:1746485/100‑1 (MQ=255)
   cttccgaTCTAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACGTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtg                                                                                            <  1:998393/94‑1 (MQ=255)
    gcatAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGACTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           >  1:2479642/4‑100 (MQ=255)
    gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           >  2:1044600/1‑100 (MQ=255)
    gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           >  1:2562974/1‑100 (MQ=255)
    gCGTAATCCAGGTCGTTCGCATGTTTAACCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           >  2:613517/1‑100 (MQ=255)
     cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCTAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          <  1:234156/100‑1 (MQ=255)
     cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          <  2:2554514/100‑1 (MQ=255)
      gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                         <  2:1486928/100‑1 (MQ=255)
      gTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTgg                                                                                         <  2:2221966/100‑1 (MQ=255)
         aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  2:1104837/100‑1 (MQ=255)
           ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    >  2:517090/1‑100 (MQ=255)
           ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTACTTTGTGGATCAg                                                                                    >  1:2735818/1‑100 (MQ=255)
            cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTCCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                   >  2:537398/1‑100 (MQ=255)
               gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATg                                                                                <  1:1830431/100‑1 (MQ=255)
               gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTATGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 >  1:2366276/1‑99 (MQ=255)
                tCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGa                                                                               <  2:262827/100‑1 (MQ=255)
                  gTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                              >  1:1377695/1‑99 (MQ=255)
                     cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          <  1:416814/100‑1 (MQ=255)
                     cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTCCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaaag                                                                            >  2:2095337/1‑96 (MQ=255)
                      gCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTgg                                                                         >  1:1593844/1‑100 (MQ=255)
                       ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCGGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        >  1:745488/1‑100 (MQ=255)
                       ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        <  2:431251/100‑1 (MQ=255)
                       ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        >  2:2172927/1‑100 (MQ=255)
                       ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGACCAGAATGAACGTGGt                                                                        >  1:3004335/1‑100 (MQ=255)
                        cTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTa                                                                       <  1:2512914/100‑1 (MQ=255)
                          gtttgtccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGAATTGTTCTGGTTTCTTCTGTTCTTTTTGGATCAGAATGAACGTGGTaga                                                                     >  2:730989/7‑100 (MQ=255)
                          gTTTAGCCCCGGTTTTGTACGCCCAGTCGAATAATTTCAGTACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     <  1:286091/100‑1 (MQ=255)
                          gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACTTGGTaga                                                                     <  1:1584146/100‑1 (MQ=255)
                          gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     <  1:1923643/100‑1 (MQ=255)
                            ttAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                   <  2:1377682/100‑1 (MQ=255)
                            ttAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAgg                                                                   >  2:1690375/1‑100 (MQ=255)
                             tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  1:1465033/100‑1 (MQ=255)
                             tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  1:1082881/100‑1 (MQ=255)
                              aGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTa                                                                 >  1:1091543/1‑100 (MQ=255)
                               gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTaa                                                                >  1:1327681/1‑100 (MQ=255)
                               gCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTa                                                                 <  2:1090116/99‑1 (MQ=255)
                                ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTAATTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                               >  1:318071/1‑100 (MQ=255)
                                 cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  2:1057525/1‑100 (MQ=255)
                                  ccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAg                                                             >  2:858761/1‑100 (MQ=255)
                                   cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                            <  2:1609151/100‑1 (MQ=255)
                                   cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                            <  2:987672/100‑1 (MQ=255)
                                   cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                            >  1:2475366/1‑100 (MQ=255)
                                      ttttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCa                                                         >  2:273920/1‑100 (MQ=255)
                                      ttttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCa                                                         >  2:455655/1‑100 (MQ=255)
                                         tGTACGCCCAGTCGAAGAATTACAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGCTCTTAGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGc                                                      >  2:1367452/1‑100 (MQ=255)
                                          gTACGCCCAGTCGAAGAATTTCAGCTTTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGTATCAGAATGACCGTGGTAGAGGTAATAGGCCATg                                                       <  2:1235464/98‑1 (MQ=255)
                                           tACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCAt                                                    >  1:2833249/1‑100 (MQ=255)
                                            aCGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATc                                                   <  1:2913270/100‑1 (MQ=255)
                                            aCGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATc                                                   <  1:2069897/100‑1 (MQ=255)
                                              gCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCtt                                                 >  2:1827788/1‑100 (MQ=255)
                                               cccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTc                                                >  1:2147083/1‑100 (MQ=255)
                                                ccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGTTAGAGGTAATAGGCCATGCATCTTCg                                               <  2:2966168/100‑1 (MQ=255)
                                                ccAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTc                                                <  1:613821/99‑1 (MQ=255)
                                                 cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGc                                              >  1:341761/1‑100 (MQ=255)
                                                    tCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCtt                                           <  1:2670334/100‑1 (MQ=255)
                                                     cGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                          >  1:200413/1‑100 (MQ=255)
                                                      gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGTATGAACGTGGGAGAGGTAATAGGCCATGTATCTTCGCCTTTc                                         <  1:2048414/100‑1 (MQ=255)
                                                      gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAGTGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  2:475669/1‑100 (MQ=255)
                                                      gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  2:1154003/1‑100 (MQ=255)
                                                        aggaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATATGCCATGCATCTTCGCCTTTCtg                                       <  2:2454869/97‑1 (MQ=255)
                                                        agaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtg                                       <  2:1091222/100‑1 (MQ=255)
                                                          aaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggn                                     <  1:812069/100‑2 (MQ=255)
                                                           aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    >  2:18987/1‑100 (MQ=255)
                                                           aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    >  1:2634331/1‑100 (MQ=255)
                                                           aTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    >  1:1049313/1‑100 (MQ=255)
                                                            tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                   >  2:232221/1‑100 (MQ=255)
                                                            tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttg                                   <  1:1008727/100‑1 (MQ=255)
                                                            tttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggtt                                    <  1:712822/99‑1 (MQ=255)
                                                             ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  1:2291079/100‑1 (MQ=255)
                                                             ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  >  2:2090808/1‑100 (MQ=255)
                                                               cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                <  1:2878080/100‑1 (MQ=255)
                                                                aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCAGGTTGGTCa                               <  2:925570/100‑1 (MQ=255)
                                                                 gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAg                              <  1:415036/100‑1 (MQ=255)
                                                                 gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAg                              <  1:1154014/100‑1 (MQ=255)
                                                                 gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGATAGAGGTAATAGGCACTGCATCTTCACCTTTCTGGTTGGTCAg                              >  1:1982501/1‑100 (MQ=255)
                                                                  cACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGa                             >  1:1955700/1‑100 (MQ=255)
                                                                   aCTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGAt                            <  1:2328563/100‑1 (MQ=255)
                                                                     ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATcc                          <  1:1631887/100‑1 (MQ=255)
                                                                        tgtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGa                       >  1:3009849/1‑100 (MQ=255)
                                                                         gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  2:1665898/1‑100 (MQ=255)
                                                                         gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATACATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      >  1:250616/1‑100 (MQ=255)
                                                                           gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    >  2:2209086/1‑100 (MQ=255)
                                                                            ccTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGa                   >  2:1874618/1‑100 (MQ=255)
                                                                              ttGTTCTGGTTTCTTCTGTTCTTTGTGGTTCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTTGCCTTTCTGGTTGGTCAGATCCTGAGCGAAg                 <  2:2446275/100‑1 (MQ=255)
                                                                              ttGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGaan                 <  2:2214920/100‑2 (MQ=255)
                                                                               tGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAgg                >  2:1670780/1‑100 (MQ=255)
                                                                                gTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGt               >  1:711729/1‑100 (MQ=255)
                                                                                gTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGt               >  2:2237502/1‑100 (MQ=255)
                                                                                  tCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGttt             >  1:1670002/1‑100 (MQ=255)
                                                                                     ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAGAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          >  2:1675568/1‑100 (MQ=255)
                                                                                      gTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTTGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:2804504/99‑1 (MQ=255)
                                                                                      gTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  2:2530085/99‑1 (MQ=255)
                                                                                        ttcttcTGTTCTTTGTTGATTAGAATTTACGTGGTAGAGGTAATATGCTATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcc       <  1:630873/100‑1 (MQ=255)
                                                                                          cttctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTcc       >  2:1755112/1‑98 (MQ=255)
                                                                                             ctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  2:938312/1‑100 (MQ=255)
                                                                                                |                                                                                                
GCTGGCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910566‑3910758

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