Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,171,680 A→G V309A (GTA→GCA)  coaA ← pantothenate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,6800AG89.6% 742.8 / 64.1 281V309A (GTA→GCA) coaApantothenate kinase
Reads supporting (aligned to +/- strand):  ref base A (10/19);  new base G (125/125);  total (136/145)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 281 total aligned reads displayed.

GGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTT  >  NC_000913/4171590‑4171770
                                                                                          |                                                                                          
   cGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCa                                                                                       <  1:7473048/93‑1 (MQ=255)
     gTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATg                                                                                     >  1:2944162/1‑93 (MQ=255)
          atatCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAg                                                                                <  1:1348942/93‑1 (MQ=255)
           tatCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGc                                                                               <  2:1211381/93‑1 (MQ=255)
            atCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCa                                                                              >  2:6582572/1‑93 (MQ=255)
             tCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCAc                                                                             <  1:7229755/93‑1 (MQ=255)
              cccTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACt                                                                            >  1:3677712/1‑93 (MQ=255)
               ccTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACtt                                                                           >  1:1203533/1‑93 (MQ=255)
               ccTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACtt                                                                           >  2:1335061/1‑93 (MQ=255)
               ccTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACtt                                                                           >  2:5400266/1‑93 (MQ=255)
                cTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACttt                                                                          >  2:833240/1‑93 (MQ=255)
                 tGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACtttt                                                                         <  1:4735518/93‑1 (MQ=255)
                 tGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACtttt                                                                         >  2:7028011/1‑93 (MQ=255)
                  gCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTc                                                                        <  2:7506180/93‑1 (MQ=255)
                  gCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTc                                                                        >  1:3841515/1‑93 (MQ=255)
                    gTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGt                                                                      >  1:1889405/1‑93 (MQ=255)
                    gTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGt                                                                      >  2:6527503/1‑93 (MQ=255)
                     tAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTc                                                                     >  2:6946115/1‑93 (MQ=255)
                     tAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTc                                                                     >  1:1673788/1‑93 (MQ=255)
                      aGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCa                                                                    <  2:6151028/93‑1 (MQ=255)
                      aGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCa                                                                    <  1:407284/93‑1 (MQ=255)
                       gTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAg                                                                   <  1:5038554/93‑1 (MQ=255)
                       gTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAg                                                                   >  1:3165294/1‑93 (MQ=255)
                       gTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAg                                                                   >  2:7244524/1‑93 (MQ=255)
                         gCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGa                                                                 <  2:555605/93‑1 (MQ=255)
                          cAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGAt                                                                <  2:5925314/93‑1 (MQ=255)
                           aGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGAtt                                                               >  1:7372864/1‑93 (MQ=255)
                           aGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGAtt                                                               <  2:3814408/93‑1 (MQ=255)
                            gAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTa                                                              >  1:3806136/1‑93 (MQ=255)
                             aaaaaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTaa                                                             <  2:1420614/93‑1 (MQ=255)
                             aaaaaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTaa                                                             <  2:1507494/93‑1 (MQ=255)
                             aaaaaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTaa                                                             <  1:3818005/93‑1 (MQ=255)
                                aaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACt                                                          <  1:2062371/93‑1 (MQ=255)
                                aaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACt                                                          >  2:5257544/1‑93 (MQ=255)
                                aaTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTGGCACTTTTCGTCAGGATTAAACt                                                          <  1:222338/93‑1 (MQ=255)
                                 aTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTg                                                         <  2:6038963/93‑1 (MQ=255)
                                 aTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTg                                                         <  1:3201697/93‑1 (MQ=255)
                                 aTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTg                                                         <  2:7048534/93‑1 (MQ=255)
                                 aTAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTg                                                         >  1:5450566/1‑93 (MQ=255)
                                  tAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTgg                                                        <  1:6098796/93‑1 (MQ=255)
                                  tAAGAAAGGGGAGTATTAGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTgg                                                        >  1:4926485/1‑93 (MQ=255)
                                   aaGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgc                                                       <  1:692315/93‑1 (MQ=255)
                                    aGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcg                                                      <  2:2727459/93‑1 (MQ=255)
                                    aGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcg                                                      <  2:3733189/93‑1 (MQ=255)
                                    aGAAAGGGGAGAATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGtcg                                                      >  1:2111010/1‑93 (MQ=255)
                                       aaGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcgcgc                                                   >  2:4805483/1‑93 (MQ=255)
                                       aaGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGgcgcgc                                                   <  1:5845964/93‑1 (MQ=255)
                                        aGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCt                                                  >  2:563974/1‑93 (MQ=255)
                                          gggAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCa                                                <  2:5046088/93‑1 (MQ=255)
                                          gggAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCa                                                <  1:5161012/93‑1 (MQ=255)
                                            gAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACg                                              <  1:4267712/93‑1 (MQ=255)
                                             aGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGa                                             <  2:3334691/93‑1 (MQ=255)
                                              gTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAg                                            <  1:4630845/93‑1 (MQ=255)
                                               tATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgt                                           >  1:4855451/1‑93 (MQ=255)
                                                aTTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgta                                          <  2:2062176/93‑1 (MQ=255)
                                                 ttCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtag                                         >  2:5832208/1‑93 (MQ=255)
                                                 ttCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtag                                         <  1:4137400/93‑1 (MQ=255)
                                                  tCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtagg                                        <  2:1807546/93‑1 (MQ=255)
                                                  tCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtagg                                        <  1:5720209/93‑1 (MQ=255)
                                                    gCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggta                                      >  1:7129089/1‑93 (MQ=255)
                                                    gCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggta                                      <  1:4567735/93‑1 (MQ=255)
                                                    gCTCCCCTGCAAATTATTTGCGTAGGCTAACCTCTTCTGCCGCATGATTCGCCCTTTTCGCCAGGATTAAACTGGCGCCCTCACGAgtaggta                                      <  1:1090845/93‑1 (MQ=255)
                                                     cTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtaggtag                                     <  1:817826/93‑1 (MQ=255)
                                                     cTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAgtagggag                                     >  1:6344281/1‑93 (MQ=255)
                                                      tCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                    <  1:6692657/93‑1 (MQ=255)
                                                      tCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                    >  2:233768/1‑93 (MQ=255)
                                                      tCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                    >  1:4458035/1‑93 (MQ=255)
                                                      tCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGa                                    >  1:4978113/1‑93 (MQ=255)
                                                       ccccTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGaa                                   >  2:4497335/1‑93 (MQ=255)
                                                       ccccTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGaa                                   >  1:4806511/1‑93 (MQ=255)
                                                        cccTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                  >  2:2947545/1‑93 (MQ=255)
                                                        cccTGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                  >  1:3813654/1‑93 (MQ=255)
                                                        cccTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAat                                  <  1:5357923/93‑1 (MQ=255)
                                                           tGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtt                               <  2:2648250/93‑1 (MQ=255)
                                                           tGCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTCTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtt                               >  1:3922074/1‑93 (MQ=255)
                                                            gCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAttt                              >  2:6265254/1‑93 (MQ=255)
                                                            gCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAttt                              >  2:364247/1‑93 (MQ=255)
                                                            gCAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAttt                              <  2:432245/93‑1 (MQ=255)
                                                             cAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtttt                             <  1:1165986/93‑1 (MQ=255)
                                                             cAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtttt                             <  2:7324363/93‑1 (MQ=255)
                                                             cAAATTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATAtttt                             <  1:4183605/93‑1 (MQ=255)
                                                               aaTTATTTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGc                           <  1:3458401/93‑1 (MQ=255)
                                                                 ttattTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCtt                         >  1:2380951/1‑93 (MQ=255)
                                                                  tattTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCttt                        >  1:5287853/1‑93 (MQ=255)
                                                                  tattTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCttt                        >  2:4149271/1‑93 (MQ=255)
                                                                   attTGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTa                       >  2:1016957/1‑93 (MQ=255)
                                                                    tttGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTaa                      >  1:1997831/1‑93 (MQ=255)
                                                                      tGCGTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGt                    <  1:6029972/93‑1 (MQ=255)
                                                                        cGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTc                  <  1:5678699/93‑1 (MQ=255)
                                                                         gTAGTCTGACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCa                 <  1:6526014/93‑1 (MQ=255)
                                                                              cTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAg            <  1:2664183/93‑1 (MQ=255)
                                                                                gACCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGtt          <  1:1580528/93‑1 (MQ=255)
                                                                                 aCCTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTg         >  2:3118696/1‑93 (MQ=255)
                                                                                  ccTCTTCTGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGa        <  2:4716724/93‑1 (MQ=255)
                                                                                    tcttctGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtc      >  1:6243399/1‑93 (MQ=255)
                                                                                      ttctGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtctc    >  2:3580074/1‑93 (MQ=255)
                                                                                      ttctGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtctc    >  2:6714912/1‑93 (MQ=255)
                                                                                      ttctGCCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGAtctc    <  2:3239781/93‑1 (MQ=255)
                                                                                        ctACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCtt  <  1:2734628/93‑1 (MQ=255)
                                                                                        ctACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCtt  >  1:1786260/1‑93 (MQ=255)
                                                                                          |                                                                                          
GGGCGGTGAAATATCCCTGCGTAGTGCAGAAAAATAAGAAAGGGGAGTATTCGCTCCCCTGCAAATTATTTGCGTAGTCTGACCTCTTCTACCGCATGATTAGCACTTTTCGTCAGGATTAAACTGGCGCGCTCACGAGTAGGTAGAATATTTTGCTTTAAGTTCAGCCAGTTGATCTCTT  >  NC_000913/4171590‑4171770

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