Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,342,306 T→C L241P (CTC→CCC)  yhbG → putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,342,3060TC88.2% 787.2 / 82.9 279L241P (CTC→CCC) yhbGputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base T (14/19);  new base C (120/126);  total (134/145)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

Only 100 of 292 total aligned reads displayed.

GTCAGGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCCAACAACT  >  NC_000913/3342216‑3342395
                                                                                          |                                                                                         
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         <  2:3381089/93‑1 (MQ=255)
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         >  1:5445743/1‑93 (MQ=255)
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         >  2:5962007/1‑93 (MQ=255)
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         <  2:5344677/93‑1 (MQ=255)
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         <  1:2823454/93‑1 (MQ=255)
gtcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGAccct                                                                                         <  2:2273974/93‑1 (MQ=255)
 tcagGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTg                                                                                        >  1:5815877/1‑93 (MQ=255)
 gcagGCGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTg                                                                                        <  1:2092908/92‑1 (MQ=255)
   agGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGAt                                                                                      >  2:3787032/1‑93 (MQ=255)
   agGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGAt                                                                                      <  2:2863494/93‑1 (MQ=255)
      ggCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAgg                                                                                   <  2:435931/93‑1 (MQ=255)
      ggCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAgg                                                                                   <  2:4194388/93‑1 (MQ=255)
       gCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAggg                                                                                  >  2:1605709/1‑93 (MQ=255)
        cATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGt                                                                                 <  1:4806235/93‑1 (MQ=255)
         aTTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTa                                                                                >  1:4900872/1‑93 (MQ=255)
         aTTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTa                                                                                >  2:900161/1‑93 (MQ=255)
           ttGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGa                                                                              <  1:4750256/93‑1 (MQ=255)
           ttGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGa                                                                              <  1:3358039/93‑1 (MQ=255)
           ttGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGa                                                                              <  1:4102510/93‑1 (MQ=255)
            tGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGaa                                                                             >  1:2522945/1‑93 (MQ=255)
             gATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACTCTGATAGGGTAGAAg                                                                            <  2:2577269/93‑1 (MQ=255)
              aTCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGt                                                                           >  1:945076/1‑93 (MQ=255)
                cGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGttt                                                                         >  1:1620389/1‑93 (MQ=255)
                cGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGttt                                                                         >  1:2774752/1‑93 (MQ=255)
                   ccACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGCATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCg                                                                      <  1:5423299/93‑1 (MQ=255)
                    cACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGa                                                                     <  2:2009006/93‑1 (MQ=255)
                     aCGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGAc                                                                    >  2:3588039/1‑93 (MQ=255)
                     aCGGCACGCCTACAGAAATCCTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGAc                                                                    >  2:6643928/1‑93 (MQ=255)
                      cGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACg                                                                   <  2:2641993/93‑1 (MQ=255)
                       ggCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGt                                                                  <  1:3281618/93‑1 (MQ=255)
                       ggCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGt                                                                  <  2:3282494/93‑1 (MQ=255)
                        gCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGtt                                                                 >  2:4638849/1‑93 (MQ=255)
                        gCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGtt                                                                 <  1:1788262/93‑1 (MQ=255)
                          aCGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGtttt                                                               <  2:3482041/93‑1 (MQ=255)
                          aCGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGtttt                                                               <  2:5116835/93‑1 (MQ=255)
                          aCGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGtttt                                                               >  2:1450447/1‑93 (MQ=255)
                           cGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTa                                                              >  1:4488762/1‑93 (MQ=255)
                           cGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTa                                                              <  2:4498255/93‑1 (MQ=255)
                            gCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAg                                                             <  2:3901478/93‑1 (MQ=255)
                            gCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAg                                                             <  2:4449754/93‑1 (MQ=255)
                            gCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAg                                                             >  1:1999430/1‑93 (MQ=255)
                             ccTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGt                                                                  <  1:6556139/87‑1 (MQ=255)
                              cTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCa                                                           <  2:507464/93‑1 (MQ=255)
                              cTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCa                                                           >  1:6146481/1‑93 (MQ=255)
                              cTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCa                                                           >  1:1075822/1‑93 (MQ=255)
                               tACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAg                                                          >  2:5827327/1‑93 (MQ=255)
                                aCAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAgg                                                         >  1:3835521/1‑93 (MQ=255)
                                  aGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGgag                                                       >  2:2947046/1‑93 (MQ=255)
                                   gAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGgaga                                                      <  2:2835842/93‑1 (MQ=255)
                                    aaaTCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGagag                                                     >  1:1604066/1‑93 (MQ=255)
                                     aaTCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGt                                                    <  1:1435648/93‑1 (MQ=255)
                                     aaTCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGt                                                    <  2:3861625/93‑1 (MQ=255)
                                      aTCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTa                                                   >  1:6250410/1‑93 (MQ=255)
                                      aTCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTa                                                   >  2:1290806/1‑93 (MQ=255)
                                       tCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTAc                                                  <  1:3601012/93‑1 (MQ=255)
                                       tCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTAc                                                  >  1:2913716/1‑93 (MQ=255)
                                       tCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTAc                                                  <  2:6454901/93‑1 (MQ=255)
                                         ttACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGa                                                <  1:2358964/93‑1 (MQ=255)
                                         ttACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGa                                                >  2:1149398/1‑93 (MQ=255)
                                         ttACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGa                                                <  1:4592724/93‑1 (MQ=255)
                                         ttACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGa                                                >  1:5588380/1‑93 (MQ=255)
                                          tACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGAt                                               >  2:5440009/1‑93 (MQ=255)
                                           aCAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGAtt                                              <  2:1538597/93‑1 (MQ=255)
                                            cAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTc                                             <  1:2654158/93‑1 (MQ=255)
                                             aaGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCt                                            >  1:5198818/1‑93 (MQ=255)
                                               gACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGa                                          <  2:4395613/93‑1 (MQ=255)
                                                aCGTACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGaa                                         >  2:6117428/1‑93 (MQ=255)
                                                aCGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGATAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGaa                                         <  1:1563199/93‑1 (MQ=255)
                                                 cGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAAc                                        >  2:3940184/1‑93 (MQ=255)
                                                      acGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGaa                                   <  2:239543/93‑1 (MQ=255)
                                                       cGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAg                                  <  2:219673/93‑1 (MQ=255)
                                                        gTTAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGc                                 <  2:1802949/93‑1 (MQ=255)
                                                        gTTAAGCGTGTATACCGTGGGGAAGGCTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGc                                 <  2:6197775/93‑1 (MQ=255)
                                                          tAAGCGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCaa                               <  1:5262335/93‑1 (MQ=255)
                                                              cGTGTATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGtt                           <  1:5539390/93‑1 (MQ=255)
                                                               gtgtATGCCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGttt                          <  1:1121554/93‑1 (MQ=255)
                                                               gtgtATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGttt                          >  1:5390658/1‑93 (MQ=255)
                                                                 gtATACCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGc                        <  2:4113590/93‑1 (MQ=255)
                                                                  tataCCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCa                       <  2:4879549/93‑1 (MQ=255)
                                                                   ataCCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCaa                      >  2:661215/1‑93 (MQ=255)
                                                                    taCCTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAAc                     <  1:3624547/93‑1 (MQ=255)
                                                                      ccTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTc                   >  1:2843910/1‑93 (MQ=255)
                                                                       cTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCa                  >  1:1148967/1‑93 (MQ=255)
                                                                       cTTGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCa                  <  1:3469922/93‑1 (MQ=255)
                                                                        ttGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAg                 >  2:6539544/1‑93 (MQ=255)
                                                                        ttGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAg                 >  2:2760479/1‑93 (MQ=255)
                                                                         tGGGGAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAgg                >  1:769804/1‑93 (MQ=255)
                                                                          ggggAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGc               <  1:4353022/93‑1 (MQ=255)
                                                                           gggAAGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAAGTTTGCAACTCAGGCt              >  1:5737899/1‑93 (MQ=255)
                                                                             gAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTa            <  1:4276586/93‑1 (MQ=255)
                                                                             gAAGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTa            >  2:5603871/1‑93 (MQ=255)
                                                                               aGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTaac          >  2:714094/1‑91 (MQ=255)
                                                                               aGACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGc          >  2:3620453/1‑93 (MQ=255)
                                                                                gACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGcc         <  1:3908877/93‑1 (MQ=255)
                                                                                gACTTCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGca         >  2:5188433/1‑92 (MQ=255)
                                                                                    tCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCcaaca     >  2:6378134/1‑93 (MQ=255)
                                                                                    tCAGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCcaaca     <  2:1892662/93‑1 (MQ=255)
                                                                                      aGACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCCaacaac   <  1:4091958/93‑1 (MQ=255)
                                                                                       gACCCTGCTAGGGAAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCCAACAACt  >  1:5977667/1‑93 (MQ=255)
                                                                                       gACCCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCCAACAACt  <  2:777237/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCAGGGGCATTTGATCGCCCACGGCACGCCTACAGAAATCTTACAAGACGAACACGTTAAGCGTGTATACCTTGGGGAAGACTTCAGACTCTGATAGGGTAGAAGTTTGCGACGTTTTAGCAGGAGAGTACGATTCTGAACATGAAGCAAGGTTTGCAACTCAGGCTTAGCCAACAACT  >  NC_000913/3342216‑3342395

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