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)
         gatACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCAc                                                                                       >  2:1664818/1‑100 (MQ=255)
          atACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATTCTTTCTGCCAGTGCCTGGGCTTTTCGCCAc                                                                                       <  1:1056435/99‑1 (MQ=255)
           tACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCg                                                                                     >  2:1131443/1‑100 (MQ=255)
           tACCGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCg                                                                                     >  2:857321/1‑100 (MQ=255)
             ccGATAATATCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATTCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGtt                                                                                   >  1:560327/1‑100 (MQ=255)
                 taataTCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAg                                                                               <  2:1326505/100‑1 (MQ=255)
                   atatCGAACAGGTTCATGGCGAAGGCGATCAGGGTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGTTTTTCGCCACCGTTGGAGTg                                                                             <  1:191928/100‑1 (MQ=255)
                    tatCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTgg                                                                            <  1:992464/100‑1 (MQ=255)
                     atCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGa                                                                           <  2:1757533/100‑1 (MQ=255)
                     atCGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGa                                                                           >  2:916283/1‑100 (MQ=255)
                      tCGAACAGGTTCATGGCGCAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGATTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGTAGTGGAt                                                                          <  1:888791/100‑1 (MQ=255)
                       cGAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGAtt                                                                         >  2:1454900/1‑100 (MQ=255)
                        gAACAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGAttt                                                                        >  2:267024/1‑100 (MQ=255)
                          aCAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTcc                                                                      <  2:394617/100‑1 (MQ=255)
                          aCAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTcc                                                                      <  1:774109/100‑1 (MQ=255)
                           cAGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCca                                                                     <  1:1183817/100‑1 (MQ=255)
                            aGGTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCcag                                                                    <  1:617898/100‑1 (MQ=255)
                             ggTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCcagc                                                                   >  1:916606/1‑100 (MQ=255)
                              gTTCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCcagca                                                                  <  1:1664807/100‑1 (MQ=255)
                                tCATGGCGAAGGCGATCAGGTTTGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCcagcag                                                                 <  1:567418/100‑1 (MQ=255)
                                tCATGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAgg                                                                <  2:1079298/100‑1 (MQ=255)
                                   tGGCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCgg                                                             <  1:1810435/100‑1 (MQ=255)
                                     gCGAAGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTc                                                           >  1:936838/1‑100 (MQ=255)
                                      cGAAGGCGATCAGGTTGGTATCGTCATCCTGCTATTCAATGCTTTCTGCAAGTGCCTGGGCTTTTCGCCACCGTTGGAGAGGATTTCCAGCAGGCGGTCa                                                          >  1:74142/1‑100 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:138204/100‑1 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:137747/100‑1 (MQ=255)
                                        aaGGCGATCAGGTTGGTATCGTCATCCGGCTCTTCAATCCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACg                                                        <  2:896546/100‑1 (MQ=255)
                                          ggCGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAc                                                      <  1:469483/100‑1 (MQ=255)
                                            cGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAccc                                                    >  2:1740570/1‑100 (MQ=255)
                                            cGATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGAccc                                                    <  1:49476/100‑1 (MQ=255)
                                             gATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCt                                                   >  2:1518363/1‑100 (MQ=255)
                                             gATCAGGTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCt                                                   >  2:305907/1‑100 (MQ=255)
                                                  ggTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              >  1:1528300/1‑100 (MQ=255)
                                                  ggTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              >  1:397034/1‑100 (MQ=255)
                                                   gTTGGTATCGTCATCCTGCTCTTCTATGCTTTCTGCCAGTGCCTGGGTTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTcc                                              <  2:1063485/99‑1 (MQ=255)
                                                   gTTGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCa                                             <  1:585944/100‑1 (MQ=255)
                                                    ttGGTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAg                                            >  2:725119/1‑100 (MQ=255)
                                                      ggTATCGTCATCCTTCTCTTCAATGCTTTCTGGTAGTGCCTGGGATTTTCGCCACCGTTGGAGGGGATTTCCAGCAGGCGGTCACGAGCCTGTTCCAGc                                           >  1:1461504/1‑99 (MQ=255)
                                                       gTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCtg                                         <  1:1410781/100‑1 (MQ=255)
                                                       gTATCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTACAGCtg                                         >  1:20476/1‑100 (MQ=255)
                                                         aTCGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAACAGGCGGTCACGACCCTGTTCCAGCtgtg                                       >  2:1488927/1‑100 (MQ=255)
                                                          tCGTCATCCTGCTCTTCAATGCTTTCTGCAAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGc                                      >  1:1042434/1‑100 (MQ=255)
                                                           cGTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCt                                     >  1:892702/1‑100 (MQ=255)
                                                            gTCATCCTTCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATGTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCtt                                    <  1:493362/100‑1 (MQ=255)
                                                            gTCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCtt                                    >  2:1258957/1‑100 (MQ=255)
                                                             tCATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCttt                                   <  2:1147423/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  1:186345/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  2:1674817/100‑1 (MQ=255)
                                                              cATCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTc                                  <  2:709114/100‑1 (MQ=255)
                                                                tCCTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAg                                >  2:95645/1‑100 (MQ=255)
                                                                 ccTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAgc                               <  1:36259/100‑1 (MQ=255)
                                                                  cTGCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAgcg                              >  1:802400/1‑100 (MQ=255)
                                                                    gCTCTTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCt                            >  2:616713/1‑100 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCAGCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         >  1:974561/1‑100 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         <  2:1446995/100‑1 (MQ=255)
                                                                       ctTCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         >  2:732493/1‑100 (MQ=255)
                                                                        ttCACTGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACTCTGTTCCAGCTGTGCTTTCAGCGCTTCa                         <  1:1711670/99‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       <  1:301279/100‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       <  2:916612/100‑1 (MQ=255)
                                                                         tCAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATg                       >  1:507395/1‑100 (MQ=255)
                                                                          cAATGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGt                      <  2:571056/100‑1 (MQ=255)
                                                                          cAATGCTTTCTACCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGt                      >  1:792785/1‑100 (MQ=255)
                                                                           aaTGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGtt                     >  1:1819521/1‑100 (MQ=255)
                                                                           aaTGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGtt                     >  2:1373027/1‑100 (MQ=255)
                                                                            aTGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTg                    >  2:1240008/1‑100 (MQ=255)
                                                                             tGCTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGc                   <  1:559097/100‑1 (MQ=255)
                                                                               cTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCt                  >  1:55063/1‑99 (MQ=255)
                                                                               cTTTCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTc                 >  1:631940/1‑100 (MQ=255)
                                                                                 ttCTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTcgc               >  2:579933/1‑100 (MQ=255)
                                                                                  tCTGCCAGTGCCTGGGTTTTTCGTCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTcgcg              <  2:309060/100‑1 (MQ=255)
                                                                                   cTGCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCgg             <  2:487453/100‑1 (MQ=255)
                                                                                     gCCAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTAATGTTGCTcgcg              >  2:450498/1‑97 (MQ=255)
                                                                                       cAGTGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGt         >  1:962185/1‑100 (MQ=255)
                                                                                          tGCCTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGtttt      <  1:822848/100‑1 (MQ=255)
                                                                                            ccTGGGCTTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGttttt     >  2:925629/1‑100 (MQ=255)
                                                                                             cTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGa   >  1:678958/1‑100 (MQ=255)
                                                                                             cTGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGa   <  2:943273/100‑1 (MQ=255)
                                                                                              tGGGCTTTTCGCCACCGTTGGAGTGGATTTCCAGCAGGCGGTCACGACCCTGTTCCAGCTGTGCTTTCAGCGCTTCATGTTGCTCGCGGCAGTTTTTGAt  >  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: