Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 61,175 (T)5→4 coding (2090/2907 nt) rapA ← RNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091361,1710T.92.3% 297.9 / 22.7 79coding (2094/2907 nt)rapARNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base . (38/34);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGAT  >  NC_000913/61072‑61266
                                                                                                   |                                                                                               
gTCCTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCt                                                                                                 >  2:163013/1‑100 (MQ=255)
  ccTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGTGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCttt                                                                                               <  1:545931/100‑1 (MQ=255)
  ccTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCttt                                                                                               >  1:415940/1‑100 (MQ=255)
   cTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCtttt                                                                                              >  1:1614897/1‑100 (MQ=255)
   cTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCtttt                                                                                              >  1:236305/1‑100 (MQ=255)
     gattgatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCttttcg                                                                                            <  1:1065502/100‑3 (MQ=255)
     gattgatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTTc                                                                                            <  1:329968/100‑1 (MQ=255)
      attgatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCttttcgc                                                                                           >  1:840432/1‑97 (MQ=255)
      attgatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCttttcgc                                                                                           <  2:579174/100‑4 (MQ=255)
         gatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCAc                                                                                       >  2:1664818/1‑100 (MQ=255)
          atACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATTCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCAc                                                                                       <  1:1056435/99‑1 (MQ=255)
           tACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCg                                                                                     >  2:1131443/1‑100 (MQ=255)
           tACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCg                                                                                     >  2:857321/1‑100 (MQ=255)
             ccGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATTCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGtt                                                                                   >  1:560327/1‑100 (MQ=255)
                 taataTCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAg                                                                               <  2:1326505/100‑1 (MQ=255)
                   atatCGAACAGGTTCATGGCGAAGGCGATCAGGGTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGG‑TTTTTCGCCACCGTTGGAGTg                                                                             <  1:191928/100‑1 (MQ=255)
                    tatCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTgg                                                                            <  1:992464/100‑1 (MQ=255)
                     atCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGa                                                                           <  2:1757533/100‑1 (MQ=255)
                     atCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGa                                                                           >  2:916283/1‑100 (MQ=255)
                      tCGAACAGGTTCATGGCGCAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGATTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGTAGTGGAt                                                                          <  1:888791/100‑1 (MQ=255)
                       cGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGAtt                                                                         >  2:1454900/1‑100 (MQ=255)
                        gAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGAttt                                                                        >  2:267024/1‑100 (MQ=255)
                          aCAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTcc                                                                      <  2:394617/100‑1 (MQ=255)
                          aCAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTcc                                                                      <  1:774109/100‑1 (MQ=255)
                           cAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCca                                                                     <  1:1183817/100‑1 (MQ=255)
                            aGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCcag                                                                    <  1:617898/100‑1 (MQ=255)
                             ggTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCcagc                                                                   >  1:916606/1‑100 (MQ=255)
                              gTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCcagca                                                                  <  1:1664807/100‑1 (MQ=255)
                                tCATGGCGAAGGCGATCAGGTTTGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCcagcag                                                                 <  1:567418/100‑1 (MQ=255)
                                tCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAgg                                                                <  2:1079298/100‑1 (MQ=255)
                                   tGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCgg                                                             <  1:1810435/100‑1 (MQ=255)
                                     gCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTc                                                           >  1:936838/1‑100 (MQ=255)
                                      cGAAGGCGATCAGGTTGGTATCGTCATCCTGCTATTCAATGCTTTCTGCAAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGAGGATTTCCAGCAGGCGGTCa                                                          >  1:74142/1‑100 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:138204/100‑1 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:137747/100‑1 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCGGCTCTTCAATCCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:896546/100‑1 (MQ=255)
                                          ggCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAc                                                      <  1:469483/100‑1 (MQ=255)
                                            cGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAccc                                                    >  2:1740570/1‑100 (MQ=255)
                                            cGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAccc                                                    <  1:49476/100‑1 (MQ=255)
                                             gATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCt                                                   >  2:1518363/1‑100 (MQ=255)
                                             gATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCt                                                   >  2:305907/1‑100 (MQ=255)
                                                  ggTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              >  1:1528300/1‑100 (MQ=255)
                                                  ggTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              >  1:397034/1‑100 (MQ=255)
                                                   gTTGGTATCGTCATCCTGCTCTTCTATGCTTTCTGCCAGTGCCTGGG‑TTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              <  2:1063485/99‑1 (MQ=255)
                                                   gTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCa                                             <  1:585944/100‑1 (MQ=255)
                                                    ttGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAg                                            >  2:725119/1‑100 (MQ=255)
                                                      ggTATCGTCATCCTTCTCTTCAATGCTTTCTGGTAGTGCCTGGG‑ATTTTCGCCACCGTTGGAGGGGATTTCCAGCAGGCGGTCACGAGCCTGTTCCAGc                                           >  1:1461504/1‑99 (MQ=255)
                                                       gTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCtg                                         <  1:1410781/100‑1 (MQ=255)
                                                       gTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTACAGCtg                                         >  1:20476/1‑100 (MQ=255)
                                                         aTCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAACAGGCGGTCACGACCCTGTTCCAGCtgtg                                       >  2:1488927/1‑100 (MQ=255)
                                                          tCGTCATCCTGCTCTTCAATGCTTTCTGCAAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGc                                      >  1:1042434/1‑100 (MQ=255)
                                                           cGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCt                                     >  1:892702/1‑100 (MQ=255)
                                                            gTCATCCTTCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATGTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCtt                                    <  1:493362/100‑1 (MQ=255)
                                                            gTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCtt                                    >  2:1258957/1‑100 (MQ=255)
                                                             tCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCttt                                   <  2:1147423/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  1:186345/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  2:1674817/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  2:709114/100‑1 (MQ=255)
                                                                tCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAg                                >  2:95645/1‑100 (MQ=255)
                                                                 ccTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAgc                               <  1:36259/100‑1 (MQ=255)
                                                                  cTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAgcg                              >  1:802400/1‑100 (MQ=255)
                                                                    gCTCTTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCt                            >  2:616713/1‑100 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCAGCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         >  1:974561/1‑100 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         <  2:1446995/100‑1 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         >  2:732493/1‑100 (MQ=255)
                                                                        ttCACTGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACTCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         <  1:1711670/99‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       <  1:301279/100‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       <  2:916612/100‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       >  1:507395/1‑100 (MQ=255)
                                                                          cAATGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGt                      <  2:571056/100‑1 (MQ=255)
                                                                          cAATGCTTTCTACCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGt                      >  1:792785/1‑100 (MQ=255)
                                                                           aaTGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGtt                     >  1:1819521/1‑100 (MQ=255)
                                                                           aaTGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGtt                     >  2:1373027/1‑100 (MQ=255)
                                                                            aTGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTg                    >  2:1240008/1‑100 (MQ=255)
                                                                             tGCTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGc                   <  1:559097/100‑1 (MQ=255)
                                                                               cTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCt                  >  1:55063/1‑99 (MQ=255)
                                                                               cTTTCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTc                 >  1:631940/1‑100 (MQ=255)
                                                                                 ttCTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTcgc               >  2:579933/1‑100 (MQ=255)
                                                                                  tCTGCCAGTGCCTGGG‑TTTTTCGTCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTcgcg              <  2:309060/100‑1 (MQ=255)
                                                                                   cTGCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCgg             <  2:487453/100‑1 (MQ=255)
                                                                                     gCCAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTAATGTTGCTcgcg              >  2:450498/1‑97 (MQ=255)
                                                                                       cAGTGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGt         >  1:962185/1‑100 (MQ=255)
                                                                                          tGCCTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGtttt      <  1:822848/100‑1 (MQ=255)
                                                                                            ccTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGttttt     >  2:925629/1‑100 (MQ=255)
                                                                                             cTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGa   >  1:678958/1‑100 (MQ=255)
                                                                                             cTGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGa   <  2:943273/100‑1 (MQ=255)
                                                                                              tGGGC‑TTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGAt  >  2:377054/1‑100 (MQ=255)
                                                                                                   |                                                                                               
GTCCTGATTGATACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGAT  >  NC_000913/61072‑61266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑