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% 1207.0 / ‑5.4 385N97S (AAC→AGC) ugpCATP‑binding component of sn‑glycerol 3‑phosphate transport system
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (215/168);  total (216/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 390 total aligned reads displayed.

CGCTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGTTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCGC  >  NC_000913/3586471‑3586650
                                                                                          |                                                                                         
  cTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTg                                                                                       <  1:5526240/93‑1 (MQ=255)
  cTCGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTg                                                                                       >  1:3485219/1‑93 (MQ=255)
    cGGCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGa                                                                                     <  1:7017067/93‑1 (MQ=255)
     ggCAATTTGCTGCTTTCCCAGGCCGCGAATATTCAGCCCCCACGCCATGTTTTCTTCGACACGCATATGCGGATAAAGCGCGTAGTTCTGGaa                                                                                    <  1:4181392/93‑1 (MQ=38)
      gCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAac                                                                                   <  1:6402566/93‑1 (MQ=255)
      gCAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAac                                                                                   <  1:303301/93‑1 (MQ=255)
       cAATTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAACGCGCGTAGCTCTGGAaca                                                                                  <  2:3189676/93‑1 (MQ=255)
        aaTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAacac                                                                                 >  2:2656809/1‑93 (MQ=255)
         aTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAGCAcc                                                                                <  2:5128573/93‑1 (MQ=255)
         aTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACAcc                                                                                >  1:1256648/1‑93 (MQ=255)
         aTTTGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACAcc                                                                                >  2:1723145/1‑93 (MQ=255)
          tttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCa                                                                               >  1:6019526/1‑93 (MQ=255)
          tttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCa                                                                               >  2:707016/1‑93 (MQ=255)
          tttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCa                                                                               >  2:4736546/1‑93 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              <  1:3974964/93‑1 (MQ=255)
           ttGCTGCTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCAt                                                                              <  2:4099671/93‑1 (MQ=255)
            tgctgcTTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATc                                                                             <  2:7632192/93‑1 (MQ=255)
                gctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCaa                                                                         >  2:1865334/1‑93 (MQ=255)
                gctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCaa                                                                         <  2:7748794/93‑1 (MQ=255)
                gctTGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCaa                                                                         >  2:5565459/1‑93 (MQ=255)
                  ttGCCCATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATc                                                                       >  1:291253/1‑93 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcg                                                                    <  2:4882492/93‑1 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcg                                                                    <  1:7192990/93‑1 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcg                                                                    >  2:3274434/1‑93 (MQ=255)
                     cccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcg                                                                    <  1:3747069/93‑1 (MQ=255)
                      ccATGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCcgc                                                                   <  2:5970523/93‑1 (MQ=255)
                         tGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGAt                                                                <  2:4570262/93‑1 (MQ=255)
                         tGCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGAGCGTAGCTCTGGAACACCATCGCAATCCCGCGAt                                                                >  2:6732548/1‑93 (MQ=255)
                          gCCGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATc                                                               <  2:5868420/93‑1 (MQ=255)
                           ccGCGAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCt                                                              >  2:7157001/1‑93 (MQ=255)
                            cgcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtt                                                             <  2:4369365/93‑1 (MQ=255)
                            cgcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtt                                                             >  2:3558801/1‑93 (MQ=255)
                             gcgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCttt                                                            >  2:6853471/1‑93 (MQ=255)
                              cgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtttt                                                           <  2:790928/93‑1 (MQ=255)
                              cgAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCtttt                                                           >  1:1562317/1‑93 (MQ=255)
                               gAATTTTCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTg                                                          <  1:3822910/93‑1 (MQ=255)
                                  ttttCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCt                                                       <  2:7450960/93‑1 (MQ=255)
                                    ttCAGCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTcc                                                     <  2:1464160/93‑1 (MQ=255)
                                        gCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAttt                                                 >  1:5048673/1‑93 (MQ=255)
                                        gCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAttt                                                 >  2:3445243/1‑93 (MQ=255)
                                        gCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAttt                                                 >  1:6880705/1‑93 (MQ=255)
                                        gCCCCCACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCAttt                                                 >  2:1026303/1‑93 (MQ=255)
                                         cccccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTc                                                <  2:3896134/93‑1 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              >  1:1594969/1‑93 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              >  2:1233455/1‑93 (MQ=255)
                                           cccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAg                                              <  2:6571393/93‑1 (MQ=255)
                                            ccACGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGt                                             >  2:968551/1‑93 (MQ=255)
                                            ccACGCCATGTTTTCTTCGACACTCATATGCGGATAAACCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGt                                             >  2:2771307/1‑93 (MQ=255)
                                              aCGCCATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCa                                           >  1:6513276/1‑93 (MQ=255)
                                                 ccATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgc                                        >  1:1107530/1‑93 (MQ=255)
                                                 ccATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgc                                        <  2:7594903/93‑1 (MQ=255)
                                                 ccATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgc                                        <  1:2509679/93‑1 (MQ=255)
                                                  cATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgcg                                       >  1:7514222/1‑93 (MQ=255)
                                                  cATGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCAcgcg                                       >  1:5780083/1‑93 (MQ=255)
                                                   aTGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACgcgc                                      >  1:7928475/1‑93 (MQ=255)
                                                   aTGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACgcgc                                      >  2:5374591/1‑93 (MQ=255)
                                                    tGTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCt                                     >  1:6089727/1‑93 (MQ=255)
                                                     gTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTg                                    >  1:1445559/1‑93 (MQ=255)
                                                     gTTTTCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTg                                    >  1:3164490/1‑93 (MQ=255)
                                                       tttCTTCGACACTCATATGCGGCTAACGCGCGTAGTTCTGGAACACCATCGCAATCCCGCGATCTTTTTGCTCCATTTCAGTCACACGCTGGt                                  >  2:4935732/1‑93 (MQ=38)
                                                       tttCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGt                                  >  2:5748816/1‑93 (MQ=255)
                                                       tttCTTCGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGt                                  >  1:2244238/1‑93 (MQ=255)
                                                           ttcGACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGtt                              >  1:45315/1‑93 (MQ=255)
                                                              gACACTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGAt                           >  2:6075775/1‑93 (MQ=255)
                                                                 acTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCa                        >  1:2368459/1‑93 (MQ=255)
                                                                  cTCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAg                       >  1:6694414/1‑93 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      >  2:2512398/1‑93 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      <  1:2699508/93‑1 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      >  1:4008385/1‑93 (MQ=255)
                                                                   tCATATGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGa                      >  2:2438616/1‑93 (MQ=255)
                                                                     atatGGGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGata                    <  1:6239965/93‑1 (MQ=255)
                                                                     atatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGata                    <  1:7818397/93‑1 (MQ=255)
                                                                     atatGCGGATAAAGCGCGCAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACTCGCTGCTCGTTGATCCAGata                    >  2:3842786/1‑93 (MQ=38)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   <  2:7016861/93‑1 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  2:7727155/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  1:521782/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  2:5484894/1‑93 (MQ=255)
                                                                      tatGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGatat                   >  1:1966475/1‑93 (MQ=255)
                                                                       atGCGGATAAAGCGCGTAGCTCTGGAGCACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATc                  >  2:2218749/1‑93 (MQ=255)
                                                                       atGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATc                  <  2:4634428/93‑1 (MQ=255)
                                                                        tGCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCg                 >  1:874617/1‑93 (MQ=255)
                                                                         gCGGATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGc                <  1:5578423/93‑1 (MQ=255)
                                                                           ggATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCt              >  1:2802299/1‑93 (MQ=255)
                                                                           ggATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCt              >  1:407555/1‑93 (MQ=255)
                                                                            gATAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCtt             <  1:5304198/93‑1 (MQ=255)
                                                                             aTAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTc            <  1:7326857/93‑1 (MQ=255)
                                                                             aTAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTc            <  2:1778247/93‑1 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  1:485878/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  2:6920182/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  1:5933892/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  1:2061049/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           >  2:3036635/1‑93 (MQ=255)
                                                                              tAAAGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCt           <  2:2853231/93‑1 (MQ=255)
                                                                               aacgcgcGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGGCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTg          <  2:5251026/90‑1 (MQ=255)
                                                                               aaaGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTg          <  1:5412381/93‑1 (MQ=255)
                                                                                aaGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGt         >  2:6271736/1‑93 (MQ=255)
                                                                                aaGCGCGTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGt         >  2:2326888/1‑93 (MQ=255)
                                                                                    gcgTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCAcc     <  1:5875851/93‑1 (MQ=255)
                                                                                      gTAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCg   >  2:5905081/1‑93 (MQ=255)
                                                                                       tAGCTCTGGAACACCATCGCAATCCCGCGATCTTTTGGCTCCATTTCAGTCACGCGCTGGTCGTTGATCCAGATATCGCCTTCTGTCACCCGc  <  2:2040719/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 < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: