Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,586,561 T→C N97S (AAC→AGC)  ugpC ← ATP‑binding component of sn‑glycerol 3‑phosphate transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,586,5610TC99.5% 688.5 / ‑5.0 219N97S (AAC→AGC) ugpCATP‑binding component of sn‑glycerol 3‑phosphate transport system
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (117/101);  total (117/102)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 227 total aligned reads displayed.

CGCTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGTTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCGC  >  NC_000913/3586471‑3586650
                                                                                          |                                                                                         
cgcTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTc                                                                                         >  2:2261475/1‑93 (MQ=255)
cgcTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTc                                                                                         <  2:3465295/93‑1 (MQ=255)
   tCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTgg                                                                                      <  2:3859320/93‑1 (MQ=255)
       cAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAaca                                                                                  >  2:1460526/1‑93 (MQ=255)
        aaTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacac                                                                                 >  1:2880009/1‑93 (MQ=255)
        aaTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacac                                                                                 >  2:517637/1‑93 (MQ=255)
        aaTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacac                                                                                 >  1:1795343/1‑93 (MQ=255)
        aaTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacac                                                                                 <  1:2615669/93‑1 (MQ=255)
         aTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacaca                                                                                >  2:2740002/1‑92 (MQ=255)
         aTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACAcc                                                                                >  1:1198076/1‑93 (MQ=255)
          tttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCa                                                                               >  2:1715486/1‑93 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              >  2:3327914/1‑93 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              <  2:2454744/93‑1 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              <  2:2330993/93‑1 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              >  2:1006650/1‑93 (MQ=255)
                gctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCaa                                                                         >  2:3171273/1‑93 (MQ=255)
                 ctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAAt                                                                        >  2:3371661/1‑93 (MQ=255)
                 ctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAAt                                                                        <  2:4653198/93‑1 (MQ=255)
                  ttGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATc                                                                       >  1:715668/1‑93 (MQ=255)
                   tGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATcc                                                                      <  1:1617888/93‑1 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACCCCATCGCAATCCcg                                                                    >  2:220041/1‑93 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcg                                                                    <  1:3453604/93‑1 (MQ=255)
                         tGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGAt                                                                <  1:2297628/93‑1 (MQ=255)
                           ccGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCt                                                              <  1:4187077/93‑1 (MQ=255)
                            cgcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtt                                                             >  2:1785405/1‑93 (MQ=255)
                            cgcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtt                                                             <  2:3945531/93‑1 (MQ=255)
                             gcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCttt                                                            <  2:4206559/93‑1 (MQ=255)
                             gcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCttt                                                            <  2:2964254/93‑1 (MQ=255)
                             gcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCttt                                                            <  1:4224357/93‑1 (MQ=255)
                             gcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCttt                                                            >  2:1086549/1‑93 (MQ=255)
                              cgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtttt                                                           <  1:600822/93‑1 (MQ=255)
                               gAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTg                                                          <  1:2660667/93‑1 (MQ=255)
                               gAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTg                                                          >  1:4653280/1‑93 (MQ=255)
                               gAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTg                                                          <  1:1317592/93‑1 (MQ=255)
                                  ttttCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCt                                                       >  1:242877/1‑93 (MQ=255)
                                     tCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCa                                                    >  1:511495/1‑93 (MQ=255)
                                     tCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGCAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCa                                                    <  1:2510398/93‑1 (MQ=255)
                                      cAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAt                                                   <  2:3086222/93‑1 (MQ=255)
                                       aGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAtt                                                  <  1:2255681/93‑1 (MQ=255)
                                         cccccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTc                                                >  2:828247/1‑93 (MQ=255)
                                          ccccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCa                                               >  2:3534974/1‑93 (MQ=255)
                                          ccccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCa                                               >  1:521411/1‑93 (MQ=255)
                                          ccccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCACGCGATCTTTTGGCTCCATTTCa                                               <  2:4342014/93‑1 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              >  2:1609048/1‑93 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              <  1:3625630/93‑1 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              >  2:863575/1‑93 (MQ=255)
                                            ccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGt                                             >  2:631149/1‑93 (MQ=255)
                                             cACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTc                                            >  1:1449729/1‑93 (MQ=255)
                                              aCGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCa                                           <  2:4012804/93‑1 (MQ=255)
                                              aCGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCa                                           >  2:1728027/1‑93 (MQ=255)
                                               cGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAc                                          <  2:2772178/93‑1 (MQ=255)
                                                gCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcg                                         >  2:587061/1‑93 (MQ=255)
                                                 ccATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgc                                        <  1:3889100/93‑1 (MQ=255)
                                                 ccATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgc                                        <  2:858875/93‑1 (MQ=255)
                                                   aTGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACgcgc                                      <  1:210582/93‑1 (MQ=255)
                                                   aTGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACgcgc                                      <  2:673050/93‑1 (MQ=255)
                                                   aTGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACgcgc                                      <  2:2091567/93‑1 (MQ=255)
                                                    tGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCt                                     >  2:1418447/1‑93 (MQ=255)
                                                    tGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCt                                     >  1:3796593/1‑93 (MQ=255)
                                                    tGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCt                                     >  1:4003455/1‑93 (MQ=255)
                                                    tGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCt                                     <  1:1885970/93‑1 (MQ=255)
                                                      ttttCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTgg                                   <  2:4156244/93‑1 (MQ=255)
                                                       tttCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGt                                  >  1:2708860/1‑93 (MQ=255)
                                                         tcttcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGGCg                                >  2:5500/1‑93 (MQ=255)
                                                          cttcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGt                               >  2:204880/1‑93 (MQ=255)
                                                           ttcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGtt                              <  1:371563/93‑1 (MQ=255)
                                                            tcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGTTCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTg                             <  2:1335148/93‑1 (MQ=255)
                                                            tcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTg                             <  1:3093335/93‑1 (MQ=255)
                                                              gACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGAt                           <  1:4318623/93‑1 (MQ=255)
                                                                 acTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATcta                        >  1:1532686/1‑91 (MQ=255)
                                                                 acTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCa                        >  2:3337943/1‑93 (MQ=255)
                                                                 acTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCa                        <  1:28306/93‑1 (MQ=255)
                                                                  cTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAg                       >  1:685201/1‑93 (MQ=255)
                                                                  cTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAg                       >  1:1036215/1‑93 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      <  1:905889/93‑1 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      <  2:3598275/93‑1 (MQ=255)
                                                                    cATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGat                     <  2:4544648/93‑1 (MQ=255)
                                                                     atatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGata                    >  1:3204900/1‑93 (MQ=255)
                                                                     atatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGata                    >  2:4088244/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  2:2372610/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  1:3383126/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  1:192909/1‑93 (MQ=255)
                                                                       atGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATc                  >  1:3563865/1‑93 (MQ=255)
                                                                       atGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATc                  <  1:2956373/93‑1 (MQ=255)
                                                                        tGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCg                 >  1:317738/1‑93 (MQ=255)
                                                                        tGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCg                 >  2:2688606/1‑93 (MQ=255)
                                                                        tGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCg                 >  1:1213811/1‑93 (MQ=255)
                                                                        tGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTCGATCCAGATATCg                 >  2:2104615/1‑93 (MQ=255)
                                                                            gATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCtt             >  2:1672228/1‑93 (MQ=255)
                                                                             aTAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTc            <  1:564754/93‑1 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           <  2:3745053/93‑1 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  2:11208/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           <  2:4217060/93‑1 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           <  2:177579/93‑1 (MQ=255)
                                                                                   cgcgTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCAc      <  2:179445/93‑1 (MQ=255)
                                                                                    gcgTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCAcc     <  1:1141040/93‑1 (MQ=255)
                                                                                     cgTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCAccc    <  2:3239556/93‑1 (MQ=255)
                                                                                      gTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCg   >  1:3017519/1‑93 (MQ=255)
                                                                                      gTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCg   >  1:4486580/1‑93 (MQ=255)
                                                                                      gTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCg   <  2:1776600/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGTTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCGC  >  NC_000913/3586471‑3586650

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: