Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,578,875 G→A G407E (GGG→GAG)  tktB → transketolase 2 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,8750GA84.5% 713.0 / 83.7 258G407E (GGG→GAG) tktBtransketolase 2 isozyme
Reads supporting (aligned to +/- strand):  ref base G (26/14);  new base A (110/108);  total (136/122)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 262 total aligned reads displayed.

GTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTT  >  NC_000913/2578785‑2578966
                                                                                          |                                                                                           
gTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACgggg                                                                                           >  1:5987572/1‑93 (MQ=255)
gTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACgagg                                                                                           <  2:3282155/93‑1 (MQ=255)
 ttCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGt                                                                                          <  1:1297688/93‑1 (MQ=255)
 ttCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGt                                                                                          <  1:6860363/93‑1 (MQ=255)
       ggATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGa                                                                                    >  2:1917963/1‑93 (MQ=255)
        gATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGaa                                                                                   <  1:7130529/93‑1 (MQ=255)
        gATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTAAGAGGTGCGTGaa                                                                                   >  1:552241/1‑93 (MQ=255)
         aTCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAt                                                                                  >  1:2842106/1‑93 (MQ=255)
          tCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAtt                                                                                 <  1:735391/93‑1 (MQ=255)
           cTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAttt                                                                                <  2:777817/93‑1 (MQ=255)
           cTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAttt                                                                                <  1:4061097/93‑1 (MQ=255)
           cTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAttt                                                                                <  1:2358664/93‑1 (MQ=255)
           cTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAAttt                                                                                <  1:4428144/93‑1 (MQ=255)
            tGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTg                                                                               >  1:6699384/1‑93 (MQ=255)
            tGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTg                                                                               >  2:5795459/1‑93 (MQ=255)
             ggCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTgg                                                                              >  2:4247163/1‑93 (MQ=255)
              gCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGc                                                                             <  1:2169119/93‑1 (MQ=255)
              gCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGc                                                                             <  1:4893364/93‑1 (MQ=255)
               cTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCa                                                                            <  2:3958639/93‑1 (MQ=255)
                tCCCAGCAACCTGACCATCTGGAAAGTTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTAAATTTGGCAt                                                                           <  2:704454/93‑1 (MQ=255)
                tCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCAt                                                                           >  1:3160228/1‑93 (MQ=255)
                tCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCAt                                                                           >  1:1981927/1‑93 (MQ=255)
                tCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCAt                                                                           >  2:4299731/1‑93 (MQ=255)
                  ccAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGa                                                                         <  1:3550142/93‑1 (MQ=255)
                   cAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGAc                                                                        <  2:4058506/93‑1 (MQ=255)
                   cAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGAc                                                                        <  1:1749073/93‑1 (MQ=255)
                       aaCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCt                                                                    >  2:3871772/1‑93 (MQ=255)
                       aaCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCt                                                                    <  1:5051138/93‑1 (MQ=255)
                       aaCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCt                                                                    >  1:5899019/1‑93 (MQ=255)
                        aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTa                                                                   >  2:3447006/1‑93 (MQ=255)
                        aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTa                                                                   >  1:1033418/1‑93 (MQ=255)
                        aCCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTa                                                                   <  1:5409418/93‑1 (MQ=255)
                          cTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATc                                                                 >  2:6471227/1‑93 (MQ=255)
                          cTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATc                                                                 >  2:1004120/1‑93 (MQ=255)
                          cTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATc                                                                 <  1:1090152/93‑1 (MQ=255)
                           tGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCg                                                                >  2:1617543/1‑93 (MQ=255)
                           tGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCg                                                                >  2:5056814/1‑93 (MQ=255)
                              ccATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCa                                                             >  2:4922502/1‑93 (MQ=255)
                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCaa                                                            <  2:530814/93‑1 (MQ=255)
                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCaa                                                            <  1:1205825/93‑1 (MQ=255)
                               cATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCaa                                                            >  2:3575044/1‑93 (MQ=255)
                                aTCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAAc                                                           <  2:284284/93‑1 (MQ=255)
                                  cTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACgg                                                         <  2:2824823/93‑1 (MQ=255)
                                   tGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGc                                                        <  2:2020517/93‑1 (MQ=255)
                                   tGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGc                                                        <  1:6368947/93‑1 (MQ=255)
                                   tGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGc                                                        >  1:6990722/1‑93 (MQ=255)
                                    ggAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCa                                                       >  2:2405649/1‑93 (MQ=255)
                                     gAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCAt                                                      >  1:5108831/1‑93 (MQ=255)
                                      aaaGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATc                                                     <  2:1331199/93‑1 (MQ=255)
                                      aaaGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATc                                                     <  1:7472237/93‑1 (MQ=255)
                                      aaaGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATc                                                     <  2:6559465/93‑1 (MQ=255)
                                      aaaGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATc                                                     <  2:140971/93‑1 (MQ=255)
                                        aGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATcgc                                                   >  2:412493/1‑93 (MQ=255)
                                          gTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCgcgc                                                 >  1:289500/1‑93 (MQ=255)
                                          gTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCgcgc                                                 >  1:5306726/1‑93 (MQ=255)
                                          gTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCgcgc                                                 >  2:5968504/1‑93 (MQ=255)
                                            tCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGcac                                               <  2:3151300/93‑1 (MQ=255)
                                            tCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGcac                                               >  1:3856661/1‑93 (MQ=255)
                                              tGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGcacca                                             >  1:4544633/1‑93 (MQ=255)
                                              tGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGcacca                                             >  2:4205219/1‑93 (MQ=255)
                                               gTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGcaccac                                            >  1:6170631/1‑93 (MQ=255)
                                                 ttCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcgg                                          <  2:4577991/93‑1 (MQ=255)
                                                  tCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcggc                                         <  2:1144075/93‑1 (MQ=255)
                                                  tCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAcggc                                         <  1:5684663/93‑1 (MQ=255)
                                                      tGAAGGAAGATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCt                                     >  2:4639772/1‑93 (MQ=255)
                                                        aaggaagATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCttt                                   >  1:3320274/1‑93 (MQ=255)
                                                        aaggaagATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCttt                                   >  2:2095548/1‑93 (MQ=255)
                                                         aggaagATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGCCCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTtg                                  <  1:6228319/93‑1 (MQ=255)
                                                           gaagATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCAAGGCGGCTTtgtg                                <  2:6762068/93‑1 (MQ=255)
                                                            aagATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGc                               >  2:6900782/1‑93 (MQ=255)
                                                            aagATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGc                               <  2:341094/93‑1 (MQ=255)
                                                             agATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGcc                              >  1:205129/1‑93 (MQ=255)
                                                             agATCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGcc                              >  2:181803/1‑93 (MQ=255)
                                                               aTCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGt                            <  2:587158/93‑1 (MQ=255)
                                                                tCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGta                           >  1:1928696/1‑93 (MQ=255)
                                                                tCCAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGta                           <  1:4370178/93‑1 (MQ=255)
                                                                  cAGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGtata                         <  2:1529644/93‑1 (MQ=255)
                                                                   aGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATAc                        >  1:3069417/1‑93 (MQ=255)
                                                                   aGCGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATAc                        >  2:6123234/1‑93 (MQ=255)
                                                                     cGGGCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACcg                      >  2:5565440/1‑93 (MQ=255)
                                                                      gggCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACcgc                     >  1:5474693/1‑93 (MQ=255)
                                                                      gggCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACcgc                     >  1:4825295/1‑93 (MQ=255)
                                                                        gCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGa                   <  2:7341328/93‑1 (MQ=255)
                                                                        gCAACTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGa                   <  2:7234631/93‑1 (MQ=255)
                                                                          aaCTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACg                 >  1:674140/1‑93 (MQ=255)
                                                                           aCTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGt                >  1:2655162/1‑93 (MQ=255)
                                                                           aCTACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGt                >  1:5814085/1‑93 (MQ=255)
                                                                             tACATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTc              >  1:2220163/1‑93 (MQ=255)
                                                                              aCATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTcc             <  1:1479578/93‑1 (MQ=255)
                                                                              aCATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTcc             >  2:26491/1‑93 (MQ=255)
                                                                               cATTCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCt            <  1:3957387/93‑1 (MQ=255)
                                                                                 ttCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa          >  2:6386461/1‑93 (MQ=255)
                                                                                 ttCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa          >  1:3440699/1‑93 (MQ=255)
                                                                                 ttCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa          >  2:7450765/1‑93 (MQ=255)
                                                                                 ttCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGa          >  2:317575/1‑93 (MQ=255)
                                                                                  tCACTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAt         <  1:3790225/93‑1 (MQ=255)
                                                                                     cTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgtt      >  2:2287227/1‑93 (MQ=255)
                                                                                     cTACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgtt      >  2:4769151/1‑93 (MQ=255)
                                                                                      tACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgttt     >  1:169685/1‑93 (MQ=255)
                                                                                      tACGAGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGAtgttt     >  2:5459169/1‑93 (MQ=255)
                                                                                          |                                                                                           
GTTCGGCGGATCTGGCTCCCAGCAACCTGACCATCTGGAAAGGTTCTGTTTCGCTGAAGGAAGATCCAGCGGGCAACTACATTCACTACGGGGTGCGTGAATTTGGCATGACCGCTATCGCCAACGGCATCGCGCACCACGGCGGCTTTGTGCCGTATACCGCGACGTTCCTGATGTTTGTT  >  NC_000913/2578785‑2578966

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