Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 434,687 A→G E109E (GAA→GAG nusB → transcription termination

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913434,6870AG99.6% 910.6 / ‑5.7 264E109E (GAA→GAGnusBtranscription termination
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (135/127);  total (136/128)
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.96e-01

Only 100 of 271 total aligned reads displayed.

GGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAAAA  >  NC_000913/434597‑434776
                                                                                          |                                                                                         
gtcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCt                                                                                         <  1:6237764/90‑1 (MQ=255)
ggcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCt                                                                                         <  2:4049835/90‑1 (MQ=255)
 gcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTg                                                                                        <  2:4612998/91‑1 (MQ=255)
 gcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTg                                                                                        <  2:1454673/91‑1 (MQ=255)
  caGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTgg                                                                                       <  1:2865154/92‑1 (MQ=255)
    gAAAAAGCAGTACTGCGCATTGCTCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCg                                                                                     >  2:725595/1‑93 (MQ=255)
    gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCg                                                                                     >  2:3990708/1‑93 (MQ=255)
    gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCg                                                                                     >  1:4381152/1‑93 (MQ=255)
     aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGa                                                                                    >  1:1057952/1‑93 (MQ=255)
      aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGaa                                                                                   >  2:5895587/1‑93 (MQ=255)
        aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAAt                                                                                 >  2:4689504/1‑93 (MQ=255)
        aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAAt                                                                                 <  1:481475/93‑1 (MQ=255)
        aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAAt                                                                                 >  2:4197111/1‑93 (MQ=255)
         aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATc                                                                                >  2:2731328/1‑93 (MQ=255)
         aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATc                                                                                <  2:5650005/93‑1 (MQ=255)
          gCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCg                                                                               >  1:2104913/1‑93 (MQ=255)
           cAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGt                                                                              <  2:3707734/93‑1 (MQ=255)
            aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGtt                                                                             <  1:6254229/93‑1 (MQ=255)
              tACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCg                                                                           >  2:547721/1‑93 (MQ=255)
              tACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCg                                                                           <  1:5953065/93‑1 (MQ=255)
               aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCgg                                                                          <  1:1362011/93‑1 (MQ=255)
                cTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgc                                                                         <  1:5151023/93‑1 (MQ=255)
                 tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgcg                                                                        <  1:4727921/93‑1 (MQ=255)
                  gcgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgcgc                                                                       <  1:2525607/93‑1 (MQ=255)
                   cgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCa                                                                      >  2:4370420/1‑93 (MQ=255)
                   cgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCa                                                                      <  2:5319360/93‑1 (MQ=255)
                   cgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCa                                                                      >  2:2946890/1‑93 (MQ=255)
                    gcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCag                                                                     >  2:2584014/1‑93 (MQ=255)
                    gcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCag                                                                     >  1:4604012/1‑93 (MQ=255)
                     cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCaga                                                                    <  2:1767093/93‑1 (MQ=255)
                      aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaa                                                                   >  1:2242208/1‑93 (MQ=255)
                       ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaag                                                                  <  2:2458713/93‑1 (MQ=255)
                       ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaag                                                                  <  2:5254338/93‑1 (MQ=255)
                       ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaag                                                                  >  2:640646/1‑93 (MQ=255)
                       ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaag                                                                  >  1:2994324/1‑93 (MQ=255)
                          cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACa                                                               <  2:1048374/93‑1 (MQ=255)
                          cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACa                                                               <  1:5767301/93‑1 (MQ=255)
                           gcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAg                                                              <  1:6460679/93‑1 (MQ=255)
                           gcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAg                                                              >  2:3487184/1‑93 (MQ=255)
                           gcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAg                                                              >  2:5944234/1‑93 (MQ=255)
                             tGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGcc                                                            <  1:5676869/93‑1 (MQ=255)
                             tGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGcc                                                            <  2:4756780/93‑1 (MQ=255)
                              gTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAGGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCa                                                           <  2:3550492/93‑1 (MQ=255)
                                aCGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATa                                                         <  1:1666412/93‑1 (MQ=255)
                                aCGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATa                                                         <  1:5185379/93‑1 (MQ=255)
                                 cGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                                        >  2:5171741/1‑93 (MQ=255)
                                  gAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAg                                                       >  2:4497784/1‑93 (MQ=255)
                                   aaCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGt                                                      >  1:5465856/1‑93 (MQ=255)
                                   aaCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGt                                                      >  2:6039068/1‑93 (MQ=255)
                                    aCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGtt                                                     >  1:5167529/1‑93 (MQ=255)
                                      tGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCg                                                   >  1:4980071/1‑93 (MQ=255)
                                           aaaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAAc                                              >  2:1994690/1‑93 (MQ=255)
                                            aaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                             <  1:183859/93‑1 (MQ=255)
                                            aaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                             <  1:838950/93‑1 (MQ=255)
                                            aaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                             >  1:4302622/1‑93 (MQ=255)
                                            aaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                             >  1:1909058/1‑93 (MQ=255)
                                             aCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACgg                                            <  1:2810193/93‑1 (MQ=255)
                                             aCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACgg                                            <  1:2605039/93‑1 (MQ=255)
                                             aCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACgg                                            >  2:2310405/1‑93 (MQ=255)
                                              cGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGc                                           <  1:1529040/93‑1 (MQ=255)
                                                  gCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTAc                                       >  2:4887770/1‑93 (MQ=255)
                                                  gCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTAc                                       >  1:1011055/1‑93 (MQ=255)
                                                  gCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTAc                                       <  1:3962340/93‑1 (MQ=255)
                                                   cGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACt                                      <  2:5691827/93‑1 (MQ=255)
                                                    gATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                                     >  1:2784587/1‑93 (MQ=255)
                                                    gATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                                     <  1:1140051/93‑1 (MQ=255)
                                                     aTGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCg                                    >  1:997665/1‑93 (MQ=255)
                                                      tgtgCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                                   <  2:3050736/93‑1 (MQ=255)
                                                       gtgCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGAt                                  <  1:445764/93‑1 (MQ=255)
                                                           cATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAag                              >  1:3098818/1‑93 (MQ=255)
                                                            aTACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagc                             >  2:5591989/1‑93 (MQ=255)
                                                            aTACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagc                             <  2:5036152/93‑1 (MQ=255)
                                                            aTACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagc                             <  2:3939446/93‑1 (MQ=255)
                                                             tACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                            <  1:1808901/93‑1 (MQ=255)
                                                             tACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                            <  1:5377730/93‑1 (MQ=255)
                                                              aCAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcag                           <  2:6169448/93‑1 (MQ=255)
                                                                 gaGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAc                        <  2:6201494/92‑1 (MQ=255)
                                                                 aaGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAc                        >  2:1962802/1‑93 (MQ=255)
                                                                  aGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAACAGCAcc                       >  1:3792061/1‑93 (MQ=255)
                                                                   gTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACAGCGTACTCGACAAAACAGCACCt                      >  2:54933/1‑93 (MQ=255)
                                                                       ccATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                  >  2:3356570/1‑93 (MQ=255)
                                                                       ccATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                  <  1:2374110/93‑1 (MQ=255)
                                                                        cATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAt                 >  2:328514/1‑93 (MQ=255)
                                                                         aTTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTTGATAAAGCAGCACCTGTGAtt                <  1:1233826/93‑1 (MQ=255)
                                                                         aTTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt                <  1:1586100/93‑1 (MQ=255)
                                                                         aTTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt                <  2:1815051/93‑1 (MQ=255)
                                                                          ttAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc               <  2:1630983/93‑1 (MQ=255)
                                                                          ttAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc               >  1:3201604/1‑93 (MQ=255)
                                                                          ttAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc               >  1:3183857/1‑93 (MQ=255)
                                                                             aCGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc            >  1:413921/1‑93 (MQ=255)
                                                                             aCGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc            >  2:130104/1‑93 (MQ=255)
                                                                              cGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc           <  1:2328463/93‑1 (MQ=255)
                                                                              cGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc           <  2:5814593/93‑1 (MQ=255)
                                                                               gAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCt          <  2:2194110/93‑1 (MQ=255)
                                                                                aaGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTa         >  2:2996156/1‑93 (MQ=255)
                                                                                aaGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTa         >  1:1375244/1‑93 (MQ=255)
                                                                                aaGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTa         <  1:5841315/93‑1 (MQ=255)
                                                                                    gATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaa     >  1:222957/1‑93 (MQ=255)
                                                                                     aTCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaaa    >  1:3517361/1‑93 (MQ=255)
                                                                                      tCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaaaa   >  1:666437/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAAAA  >  NC_000913/434597‑434776

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