Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 62,377 A→G 5.4% Y296Y (TAT→TAC 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_00091362,3770AG 5.4% 231.7 / 2.9 74Y296Y (TAT→TACrapARNA polymerase remodeling/recycling factor ATPase; RNA polymerase‑associated, ATP‑dependent RNA translocase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base G (2/2);  total (39/35)
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.46e-01

CTCTTCAACGAACTGCGCAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGCTGTTCGGT  >  NC_000913/62230‑62523
                                                                                                                                                   |                                                                                                                                                  
ctctTCAACGAACTGCGCAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGa                                                                                                                                                    <  6:431082/148‑1 (MQ=255)
       aCGAACTGCGCAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCAc                                                                                                                                             <  7:444701/148‑1 (MQ=255)
        cGAACTGCGCAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACg                                                                                                                                            >  7:112355/1‑148 (MQ=255)
              gcgcAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTgg                                                                                                                                      >  8:119235/1‑148 (MQ=255)
              gcgcAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTgg                                                                                                                                      >  7:212138/1‑148 (MQ=255)
               cgcAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGTGAGCGGTCAGCAGCAGAACGCACGGCACGAGCTATGGCAGTTGTTCTATTGAGTGATATTGACGGCTTGgc                                                                                                                                     >  2:284700/1‑148 (MQ=255)
                gcAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCCGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGgcg                                                                                                                                    >  4:97612/1‑148 (MQ=255)
                   aaaTCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCTGGGTCGCGGGCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCAt                                                                                                                                 <  3:97612/148‑1 (MQ=255)
                     aTCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCt                                                                                                                               <  1:244112/148‑1 (MQ=255)
                     aTCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCt                                                                                                                               <  1:417358/148‑1 (MQ=255)
                        gTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTc                                                                                                                             <  1:179437/147‑1 (MQ=255)
                        gTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTc                                                                                                                             <  1:249291/147‑1 (MQ=255)
                            aaaCGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGGTATTCACGGCTTGGCGCATCTTCGCTc                                                                                                                         >  7:22643/1‑147 (MQ=255)
                            aaaCGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGGTATTCACGGCTTGGCGCATCTTCGCTc                                                                                                                         >  7:16021/1‑147 (MQ=255)
                             aaCGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGa                                                                                                                                                    <  5:407316/119‑1 (MQ=255)
                             aaCGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGa                                                                                                                                                    >  6:407316/1‑119 (MQ=255)
                             aaCGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCACATCTTCGCTCca                                                                                                                       >  7:50911/1‑148 (MQ=255)
                                  ttggggTCCAGCAGACGCAGTCGGGCGAAGTGGCTTTCCATCCGCAGCTGTTGCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTGAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCa                                                                                                                  <  5:226995/145‑1 (MQ=255)
                                      ggTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTg                                                                                                              <  5:212954/148‑1 (MQ=255)
                                       gTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGCCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTg                                                                                                              >  3:6988/1‑147 (MQ=255)
                                              agacgcagacgGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCa                                                                                                                                  >  7:352066/1‑120 (MQ=255)
                                              agacgcagacgGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCa                                                                                                                                  <  8:352066/120‑1 (MQ=255)
                                              agacgcagacgGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCCGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCtt                                                                                                      <  8:13207/148‑1 (MQ=255)
                                                     gacgGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTAAGCAGCAGAAAGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGATCCACACCAGGTCATGCGCTTCATCGac                                                                                               >  5:419994/1‑148 (MQ=255)
                                                         ggCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATATTCGCTCCACACCAGGTGATGCGCTTCATCGaccac                                                                                            >  1:166026/1‑147 (MQ=255)
                                                           cGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACcag                                                                                         <  8:162079/148‑1 (MQ=255)
                                                            gAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACcagc                                                                                        <  8:297471/148‑1 (MQ=255)
                                                                   cTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTccc                                                                                 <  8:112355/148‑1 (MQ=255)
                                                                    tttCCATCCCCAGCTGTTCCGGGGTCGCTGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATgcgc                                                                                                        <  7:111764/124‑1 (MQ=255)
                                                                      tCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCAt                                                                               <  7:119235/147‑1 (MQ=255)
                                                                         aTCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGACACGTGCTCTGCCAGTTGTTCAATGGCCTGAGATTCCCGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCgg                                                                           >  3:103479/1‑148 (MQ=255)
                                                                              cAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTca                                                                      >  8:75217/1‑148 (MQ=255)
                                                                                    ttCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAgaga                                                                <  7:75217/148‑1 (MQ=255)
                                                                                     tCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAGTGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGAt                                                               >  1:91327/1‑148 (MQ=255)
                                                                                     tCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTccc                                                                                 <  1:416319/130‑1 (MQ=255)
                                                                                     tCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTccc                                                                                 >  2:416319/1‑130 (MQ=255)
                                                                                     tCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGAGATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGAt                                                               >  5:171143/1‑148 (MQ=255)
                                                                                          ggTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTc                                                                                                     <  7:233507/105‑1 (MQ=255)
                                                                                          ggTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTc                                                                                                     >  8:233507/1‑105 (MQ=255)
                                                                                           gTCGCGGTCAGCAGCAGAACGCGCGGCACGTGCTCTGCAAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCATGTCCCATTCGGCTTCACAGAGATGTTCCa                                                         <  6:419994/148‑1 (MQ=255)
                                                                                             cgcgGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGacca                                                                                             <  2:316228/110‑1 (MQ=255)
                                                                                             cgcgGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGacca                                                                                             >  1:316228/1‑110 (MQ=255)
                                                                                               cgGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACTGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGgc                                                      <  6:390345/147‑1 (MQ=255)
                                                                                               cgGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGgc                                                      >  5:390345/1‑147 (MQ=255)
                                                                                                 gTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCt                                                   >  3:406114/1‑148 (MQ=255)
                                                                                                   caggagcagAACGCCAGGGACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTGCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGt                                                 <  5:97495/148‑1 (MQ=255)
                                                                                                     gcagcagAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACTGAGATGTTCCAGGCGCTGttt                                               <  4:6988/148‑1 (MQ=255)
                                                                                                     gcagcagAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGttt                                               <  8:212138/148‑1 (MQ=255)
                                                                                                       agcagAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGc                                             >  1:191564/1‑148 (MQ=255)
                                                                                                          agAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGttt                                               >  3:410933/1‑143 (MQ=255)
                                                                                                          agAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGttt                                               <  4:410933/143‑1 (MQ=255)
                                                                                                            aaCGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTcc                                                                                                                        <  2:382365/68‑1 (MQ=255)
                                                                                                            aaCGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTcc                                                                                                                        >  1:382365/1‑68 (MQ=255)
                                                                                                             aCGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTacgac                                       >  6:70675/1‑148 (MQ=255)
                                                                                                                     cACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTacgacg                                      >  3:411010/1‑141 (MQ=255)
                                                                                                                     cACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTacgacg                                      <  4:411010/141‑1 (MQ=255)
                                                                                                                        gTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATc                             >  6:285318/1‑147 (MQ=255)
                                                                                                                        gTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCACATTCGGCTACACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATcc                            >  5:365257/1‑148 (MQ=255)
                                                                                                                         tGCTCTGCCAGTTGTTCAATGGCCTGATATTCACgg                                                                                                                                           >  8:178626/1‑36 (MQ=255)
                                                                                                                         tGCTCTGCCAGTTGTTCAATGGCCTGATATTCACgg                                                                                                                                           <  7:178626/36‑1 (MQ=255)
                                                                                                                         tGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCa                           >  8:167840/1‑148 (MQ=255)
                                                                                                                             ctGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGc                         >  5:33349/1‑146 (MQ=255)
                                                                                                                               gCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCg                        >  3:160968/1‑145 (MQ=255)
                                                                                                                               gCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCg                        <  4:160968/145‑1 (MQ=255)
                                                                                                                               gCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAg                      >  7:334591/1‑147 (MQ=255)
                                                                                                                                ccAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGc                     >  7:2016/1‑147 (MQ=255)
                                                                                                                                 cAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCa                    >  8:419003/1‑147 (MQ=255)
                                                                                                                                     tgttCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATc                <  1:425865/147‑1 (MQ=255)
                                                                                                                                     tgttCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATc                >  2:425865/1‑147 (MQ=255)
                                                                                                                                       ttCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAATTCCAGCGAGAAAATAAc              >  5:167073/1‑147 (MQ=255)
                                                                                                                                        tCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTcc                                                          >  6:209462/1‑102 (MQ=255)
                                                                                                                                        tCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTcc                                                          <  5:209462/102‑1 (MQ=255)
                                                                                                                                          aaTGGCCTGGTATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGc          <  8:16021/148‑1 (MQ=255)
                                                                                                                                          aaTGGCCTGGTATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGc          <  8:22643/148‑1 (MQ=255)
                                                                                                                                              gCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGCTGt       <  4:406114/147‑1 (MQ=255)
                                                                                                                                               ccTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGCTGTTc     <  2:91327/148‑1 (MQ=255)
                                                                                                                                                  gATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGCTGTTCGgt  >  4:408635/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTCTTCAACGAACTGCGCAAAATCGTGGAAACGGTTCGGGTCCAGCAGACGCAGACGGGCGAAGTGGCTTTCCATCCCCAGCTGTTCCGGGGTCGCGGTCAGCAGCAGAACGCCCGGCACGTGCTCTGCCAGTTGTTCAATGGCCTGATATTCACGGCTTGGCGCATCTTCGCTCCACACCAGGTGATGCGCTTCATCGACCACCAGCAGGTCCCATTCGGCTTCACAGAGATGTTCCAGGCGCTGTTTGCTACGACGGGCAAAATCCAGCGAGCAAATCACCAGCTGTTCGGT  >  NC_000913/62230‑62523

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