Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,172,320 T→C S96G (AGT→GGT)  coaA ← pantothenate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,172,3200TC89.6% 723.5 / 63.4 252S96G (AGT→GGT) coaApantothenate kinase
Reads supporting (aligned to +/- strand):  ref base T (11/15);  new base C (129/96);  total (141/111)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 261 total aligned reads displayed.

TGTAGTGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCTTATA  >  NC_000913/4172228‑4172411
                                                                                            |                                                                                           
tGTAGTGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACt                                                                                             <  1:988845/93‑1 (MQ=255)
   aGTGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGcc                                                                                          >  2:6001392/1‑93 (MQ=255)
   aGTGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGcc                                                                                          <  1:548820/93‑1 (MQ=255)
     tGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAg                                                                                        <  1:1044936/93‑1 (MQ=255)
      gATCAGTTCAAAACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGc                                                                                       >  1:1747831/1‑93 (MQ=255)
       aTCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCa                                                                                      >  2:3729019/1‑93 (MQ=255)
        tCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCaa                                                                                     >  1:4185197/1‑93 (MQ=255)
        tCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCaa                                                                                     <  2:2886576/93‑1 (MQ=255)
        tCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCaa                                                                                     <  2:318302/93‑1 (MQ=255)
           gTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATAc                                                                                  >  1:2121809/1‑93 (MQ=255)
            ttCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACt                                                                                 >  2:2110845/1‑93 (MQ=255)
             tCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTg                                                                                <  2:6277614/93‑1 (MQ=255)
               aaCACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGAt                                                                              >  2:3297380/1‑93 (MQ=255)
                acacGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATa                                                                             >  1:4503844/1‑93 (MQ=255)
                acacGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATa                                                                             >  2:4818964/1‑93 (MQ=255)
                 cacGGCGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATaa                                                                            <  1:2933693/93‑1 (MQ=255)
                 cacGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATaa                                                                            <  1:3722035/93‑1 (MQ=255)
                  acgacgATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAAt                                                                           <  2:3778970/93‑1 (MQ=255)
                      cgaTGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTaa                                                                       <  2:1263011/93‑1 (MQ=255)
                       gaTGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAg                                                                      >  2:5455542/1‑93 (MQ=255)
                          gTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGaa                                                                   >  1:5024866/1‑93 (MQ=255)
                          gTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGaa                                                                   >  1:4874260/1‑93 (MQ=255)
                           ttCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAAt                                                                  <  1:2319906/93‑1 (MQ=255)
                            tCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATg                                                                 <  2:2748760/93‑1 (MQ=255)
                             ccGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGc                                                                >  2:5912519/1‑93 (MQ=255)
                              cGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCg                                                               <  1:4588415/93‑1 (MQ=255)
                               ggCCAACGGCTTAATAGCGCCTGCAATGCACGGGCTGTTGTCCTTTCCCCCACCGCGACGCCGCCAGCAATACTGATAATGTAAGGAATGCGt                                                              <  2:4608267/93‑1 (MQ=255)
                                gCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGtt                                                             >  1:4056512/1‑93 (MQ=255)
                                gCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGtt                                                             <  2:491649/93‑1 (MQ=255)
                                gCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGtt                                                             >  1:4812404/1‑93 (MQ=255)
                                 ccAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTg                                                            >  1:5905661/1‑93 (MQ=255)
                                 ccAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTg                                                            >  1:1232318/1‑93 (MQ=255)
                                   aaCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGcc                                                          >  2:2190791/1‑93 (MQ=255)
                                   aaCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGcc                                                          >  1:2849546/1‑93 (MQ=255)
                                   aaCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGcc                                                          >  2:5377154/1‑93 (MQ=255)
                                    aCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGccc                                                         >  1:3284337/1‑93 (MQ=255)
                                    aCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGccc                                                         <  1:1683300/93‑1 (MQ=255)
                                    aCGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGccc                                                         <  2:894215/93‑1 (MQ=255)
                                     cGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCg                                                        >  1:5619790/1‑93 (MQ=255)
                                      ggCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACGGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGt                                                       >  2:671745/1‑93 (MQ=255)
                                      ggCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGt                                                       >  2:4127009/1‑93 (MQ=255)
                                         ttAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTgg                                                    >  2:2038371/1‑93 (MQ=255)
                                         ttAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTgg                                                    >  1:1397001/1‑93 (MQ=255)
                                            aTAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTAc                                                 >  1:3356836/1‑93 (MQ=255)
                                             tAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTAcc                                                <  1:428065/93‑1 (MQ=255)
                                              aGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCa                                               <  2:4653101/93‑1 (MQ=255)
                                                 gcCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGa                                            <  2:1219887/93‑1 (MQ=255)
                                                 gcCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGa                                            <  1:3934097/93‑1 (MQ=255)
                                                    tGCAATACACGGGCGGTTGTACTTTCCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAAc                                         <  1:1582547/93‑1 (MQ=255)
                                                     gCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACt                                        >  1:1629419/1‑93 (MQ=255)
                                                     gCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACt                                        >  1:3023935/1‑93 (MQ=255)
                                                      cAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTg                                       >  2:3593876/1‑93 (MQ=255)
                                                       aaTACACGGGCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGt                                      >  1:2177500/1‑93 (MQ=255)
                                                       aaTACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAATTGt                                      >  1:6064697/1‑93 (MQ=255)
                                                       aaTACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGt                                      >  1:6528604/1‑93 (MQ=255)
                                                        aTACACGGGCGTTTGTGCTTTTCCCCACCGCGCCACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGtt                                     <  1:2606473/93‑1 (MQ=255)
                                                         tACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTc                                    >  2:674787/1‑93 (MQ=255)
                                                         tACACGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTc                                    >  2:1587805/1‑93 (MQ=255)
                                                          acacGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTcc                                   >  2:5614482/1‑93 (MQ=255)
                                                          acacGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTcc                                   >  2:2629125/1‑93 (MQ=255)
                                                           cacGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCa                                  >  2:2658939/1‑93 (MQ=255)
                                                            acGGGCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAg                                 >  2:1168240/1‑93 (MQ=255)
                                                             cGGGCGGTTGTACTTTTCCCAACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGa                                >  1:3776712/1‑93 (MQ=255)
                                                               ggCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAAc                              >  1:3685968/1‑93 (MQ=255)
                                                               ggCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAAc                              <  1:6109125/93‑1 (MQ=255)
                                                               ggCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAAc                              >  1:6261421/1‑93 (MQ=255)
                                                               ggCGGTTGTACTTTTCCCCACCACGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAAc                              <  2:5682202/93‑1 (MQ=255)
                                                                gCGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACt                             <  2:1617182/93‑1 (MQ=255)
                                                                 cGGTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTg                            >  2:5787039/1‑93 (MQ=255)
                                                                  ggTTGTACTTTTCCCCCCCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGc                           <  2:4419440/93‑1 (MQ=255)
                                                                   gTTGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGcc                          >  2:338242/1‑93 (MQ=255)
                                                                    ttGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCt                         <  2:1389153/93‑1 (MQ=255)
                                                                    ttGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCt                         <  1:5749722/93‑1 (MQ=255)
                                                                    ttGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCt                         >  2:6518795/1‑93 (MQ=255)
                                                                     tGTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTg                        >  1:5111401/1‑93 (MQ=255)
                                                                      gTACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGa                       >  2:6035464/1‑93 (MQ=255)
                                                                       tACTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGAc                      <  2:5139134/93‑1 (MQ=255)
                                                                         cTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACgg                    <  1:6346168/93‑1 (MQ=255)
                                                                         cTTTTCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACgg                    >  2:3037122/1‑93 (MQ=255)
                                                                          ttttCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgc                   <  1:4432189/93‑1 (MQ=255)
                                                                          ttttCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgc                   <  2:3443104/93‑1 (MQ=255)
                                                                          ttttCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgc                   >  2:2906050/1‑93 (MQ=255)
                                                                           tttCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgcg                  >  2:1130253/1‑93 (MQ=255)
                                                                            ttCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgcgc                 >  1:3440720/1‑93 (MQ=255)
                                                                            ttCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgcgc                 <  1:774598/93‑1 (MQ=255)
                                                                            ttCCCCACCGCGACACCGCCAGCAACACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGgccc                 <  2:6423055/93‑3 (MQ=255)
                                                                             tCCCCACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCa                <  1:581237/93‑1 (MQ=255)
                                                                               cccACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGa              <  1:2657969/93‑1 (MQ=255)
                                                                               cccACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGa              <  2:1874759/93‑1 (MQ=255)
                                                                               cccACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGa              >  1:2420033/1‑93 (MQ=255)
                                                                                ccACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGAt             >  2:4559270/1‑93 (MQ=255)
                                                                                 cACCGCGACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGAtt            <  1:5327405/93‑1 (MQ=255)
                                                                                     gcgACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGACACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAg        >  2:176879/1‑93 (MQ=255)
                                                                                     gcgACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAg        >  1:2971887/1‑93 (MQ=255)
                                                                                     gcgACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAg        >  2:4211481/1‑93 (MQ=255)
                                                                                      cgACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGc       >  1:11596/1‑93 (MQ=255)
                                                                                       gACACCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCt      <  1:1947340/93‑1 (MQ=255)
                                                                                        acacCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCtt     <  1:904733/93‑1 (MQ=255)
                                                                                         cacCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCTta    >  1:1003597/1‑93 (MQ=255)
                                                                                         cacCGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCTta    <  2:521673/93‑1 (MQ=255)
                                                                                            |                                                                                           
TGTAGTGATCAGTTCAACACGACGATGTTCCGGCCAACGGCTTAATAGCGCCTGCAATACACGGGCGGTTGTACTTTTCCCCACCGCGACACTGCCAGCAATACTGATAATGTAAGGAATGCGTTGCCCGTTGGTACCAAGAAACTGTTCCAGAACTGCCTGACGGCGCAGATTCGAGCTTATA  >  NC_000913/4172228‑4172411

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: