Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,055,697 T→C D245G (GAT→GGT)  serA ← D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,055,6970TC99.7% 998.6 / ‑5.1 319D245G (GAT→GGT) serAD‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (146/171);  total (148/171)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

Only 100 of 332 total aligned reads displayed.

TCGCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATATCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGA  >  NC_000913/3055607‑3055785
                                                                                          |                                                                                        
tCGCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATAccc                                                                                        >  2:2542834/1‑93 (MQ=255)
tCGCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATAccc                                                                                        >  2:7505301/1‑93 (MQ=255)
 cGCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACcca                                                                                       <  1:3266431/93‑1 (MQ=255)
  gCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACccac                                                                                      >  2:374792/1‑93 (MQ=255)
     aTTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCCCACAGCGCCGGAATACCcaccac                                                                                   <  2:5405835/93‑1 (MQ=255)
      tgtgTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACcaca                                                                                  <  1:6459373/90‑1 (MQ=255)
         gTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTa                                                                               >  1:4205119/1‑93 (MQ=255)
         gTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTa                                                                               >  1:953936/1‑93 (MQ=255)
           cGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCGGCGCATCACTCCGCGCCGGAATACCCCCCACAGttca                                                                             <  1:7042075/93‑4 (MQ=255)
            gCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcg                                                                            >  1:696724/1‑93 (MQ=255)
             ccGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgc                                                                           >  2:7829212/1‑93 (MQ=255)
               ggTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgc                                                                         <  1:6634428/93‑1 (MQ=255)
               ggTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgc                                                                         >  2:6141823/1‑93 (MQ=255)
                gTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgcg                                                                        >  2:5119005/1‑93 (MQ=255)
                gTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgcg                                                                        >  1:5185082/1‑93 (MQ=255)
                gTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgcg                                                                        >  2:3311993/1‑93 (MQ=255)
                gTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACcgcgcg                                                                        >  1:7884419/1‑93 (MQ=255)
                 ttCCGTCGGGAATACGTCGATTGCCGCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGaa                                                                      >  2:1134328/1‑93 (MQ=255)
                  tCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGaa                                                                      >  1:7399050/1‑93 (MQ=255)
                  tCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGaa                                                                      <  2:3267591/93‑1 (MQ=255)
                  tCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGaa                                                                      <  2:4720318/93‑1 (MQ=255)
                   ccGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAg                                                                     >  1:1661562/1‑93 (MQ=255)
                      tCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCAt                                                                  <  1:264076/93‑1 (MQ=255)
                      tCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCAt                                                                  >  1:1327462/1‑93 (MQ=255)
                       cGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCAtt                                                                 <  2:2179484/93‑1 (MQ=255)
                        gggAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTa                                                                >  2:5574026/1‑93 (MQ=255)
                         ggAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTaa                                                               >  2:1903975/1‑93 (MQ=255)
                         ggAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTaa                                                               >  1:226005/1‑93 (MQ=255)
                         gagaTACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTaa                                                               <  1:2365467/90‑1 (MQ=255)
                          gAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAAt                                                              <  1:2123893/93‑1 (MQ=255)
                            atgcGTCGATTGCCGCCGCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATca                                                            <  1:7314940/90‑1 (MQ=255)
                            aTACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATca                                                            <  2:246956/93‑1 (MQ=255)
                            aTACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATca                                                            <  2:3445488/93‑1 (MQ=255)
                             tACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATcag                                                           >  1:3730002/1‑93 (MQ=255)
                             tACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATcag                                                           <  1:7253874/93‑1 (MQ=255)
                             tACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATcag                                                           >  1:1558185/1‑93 (MQ=255)
                               cGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCCCCACAGTACCGCGCGAAGCATTAATcagca                                                         <  2:659155/93‑1 (MQ=255)
                               cGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATcagca                                                         <  2:1711309/93‑1 (MQ=255)
                               cGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATcagca                                                         <  1:5494907/93‑1 (MQ=255)
                                   gATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGa                                                     >  2:855499/1‑93 (MQ=255)
                                   gATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGa                                                     <  1:1925583/93‑1 (MQ=255)
                                     ttGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGc                                                   <  2:1946711/93‑1 (MQ=255)
                                     ttGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGc                                                   >  1:832423/1‑93 (MQ=255)
                                        ccgccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCgg                                                <  2:3361917/93‑1 (MQ=255)
                                        ccgccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCgg                                                <  1:1133537/93‑1 (MQ=255)
                                        ccgccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCgg                                                <  1:1249279/93‑1 (MQ=255)
                                        ccgccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCgg                                                <  2:2013345/93‑1 (MQ=255)
                                          gccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGATCCGGGc                                              <  1:1824267/93‑1 (MQ=255)
                                          gccCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGc                                              <  1:1423683/93‑1 (MQ=255)
                                          ccccccGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGc                                              <  1:7433699/92‑1 (MQ=255)
                                            ccccGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCtt                                            <  2:1312973/93‑1 (MQ=255)
                                             cccGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCtgc                                           >  1:4951042/1‑91 (MQ=255)
                                             cccGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTc                                           >  2:574030/1‑93 (MQ=255)
                                              ccGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCa                                          <  2:5277757/93‑1 (MQ=255)
                                               cGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATGCCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCAt                                         <  2:922438/93‑1 (MQ=255)
                                               cGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCAt                                         >  1:4120494/1‑93 (MQ=255)
                                                 ccAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTa                                       <  1:4270713/93‑1 (MQ=255)
                                                 ccAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTa                                       <  1:4362382/93‑1 (MQ=255)
                                                  cAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAg                                      <  2:4435179/93‑1 (MQ=255)
                                                   aGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGt                                     <  2:1415019/93‑1 (MQ=255)
                                                   aGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGt                                     >  1:2891245/1‑93 (MQ=255)
                                                   aGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGt                                     <  1:5229158/93‑1 (MQ=255)
                                                   aGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGt                                     <  2:2877795/93‑1 (MQ=255)
                                                     aTGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGa                                   >  2:3646751/1‑93 (MQ=255)
                                                     aTGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGa                                   <  1:5456161/93‑1 (MQ=255)
                                                      tGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGaa                                  >  2:3110193/1‑93 (MQ=255)
                                                       gTTTGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGaaa                                 <  2:7484414/93‑1 (MQ=255)
                                                        tttGCTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAt                                >  1:5101890/1‑93 (MQ=255)
                                                            gtggCCAGCGCATAACACAGCGCAGGAATACCCACCGCAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttct                            <  2:4212781/90‑1 (MQ=255)
                                                            cTCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttct                            <  2:5350021/93‑1 (MQ=255)
                                                             tCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttctt                           <  2:58053/93‑1 (MQ=255)
                                                             tCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttctt                           <  1:5810373/93‑1 (MQ=255)
                                                             tCGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttctt                           >  1:6247012/1‑93 (MQ=255)
                                                              cgacAGCGCATCACACAGCCCCGGAATCCCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGTCTTCATTAGTGAAAtttcttt                          <  1:3075791/90‑1 (MQ=255)
                                                              cGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttcttt                          <  2:2122169/93‑1 (MQ=255)
                                                              cGCCAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAAtttcttt                          <  1:1030583/93‑1 (MQ=255)
                                                                ccAGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCAAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTcg                        <  2:775815/93‑1 (MQ=255)
                                                                  aGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTcgcg                      >  1:4072411/1‑93 (MQ=255)
                                                                  aGCGCATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTcgcg                      >  2:3535045/1‑93 (MQ=255)
                                                                     gcATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGccc                   <  1:7495548/93‑1 (MQ=255)
                                                                     gcATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGccc                   >  1:3591581/1‑93 (MQ=255)
                                                                      cATCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCca                  <  2:4931922/93‑1 (MQ=255)
                                                                        tCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCcatc                <  1:3040839/93‑1 (MQ=255)
                                                                        tCACACAGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCcatc                >  1:2836960/1‑93 (MQ=255)
                                                                         cacacaGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCcatca               <  1:4604065/93‑1 (MQ=255)
                                                                           cacaGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCata             >  2:6974114/1‑93 (MQ=255)
                                                                           cacaGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCata             <  1:7199810/93‑1 (MQ=255)
                                                                           cacaGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCata             <  2:1059142/93‑1 (MQ=255)
                                                                             caGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATAtt           <  1:4816098/93‑1 (MQ=255)
                                                                             caGCGCCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATAtt           <  1:2961715/93‑1 (MQ=255)
                                                                              atctccGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATAttt          <  2:5446593/89‑1 (MQ=255)
                                                                                cgcCGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATAttttt        <  2:2367454/93‑1 (MQ=255)
                                                                                  ccGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTtgg      >  2:7090449/1‑93 (MQ=255)
                                                                                  ccGGAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTtgg      <  1:6538297/93‑1 (MQ=255)
                                                                                   cGGAATATCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTtggt     >  2:668667/1‑93 (MQ=255)
                                                                                     gAATACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTtggtgg   >  2:7059088/1‑93 (MQ=255)
                                                                                      aaTACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGa  >  2:4137039/1‑93 (MQ=255)
                                                                                      aaTACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGa  >  2:3959265/1‑93 (MQ=255)
                                                                                      aaTACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGa  <  2:6468803/93‑1 (MQ=255)
                                                                                      aaTACCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGa  <  2:2777421/93‑1 (MQ=255)
                                                                                          |                                                                                        
TCGCTATTGGTCGCCGGTTCCGTCGGGAATACGTCGATTGCCGCCCCCGCCAGATGTTTGCTCGCCAGCGCATCACACAGCGCCGGAATATCCACCACAGTACCGCGCGAAGCATTAATCAGCAGCGAGCCGGGCTTCATTAGTGAAATTTCTTTCGCGCCCATCATATTTTTGGTGGA  >  NC_000913/3055607‑3055785

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