Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,754,208 C→T 5.4% E31K (GAA→AAA)  ratB ← UPF0125 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,754,2080CT 5.4% 241.4 / 2.4 74E31K (GAA→AAA) ratBUPF0125 family protein
Reads supporting (aligned to +/- strand):  ref base C (37/33);  new base T (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.86e-01

GAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCATATTAGCCGCCAGC  >  NC_000913/2754065‑2754351
                                                                                                                                               |                                                                                                                                               
gagaGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCAGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTTAACCg                                                                                                                                            >  6:327616/1‑149 (MQ=255)
  gagGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTc                                                                                                                                          <  1:5509/149‑1 (MQ=255)
        gATAAATCTCCACCCGATCGTCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGAGccc                                                                                                                                    >  5:31870/1‑149 (MQ=255)
          tAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCcctc                                                                                                                                  <  3:212549/149‑1 (MQ=255)
            aaTCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCcctcct                                                                                                                                <  7:115907/149‑1 (MQ=255)
            aaTCTCCACCCGATCGCCATCATGCACACTATCGCTGAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCcctcct                                                                                                                                <  1:266761/149‑1 (MQ=255)
              tctcCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTg                                                                                                                               >  7:382272/1‑148 (MQ=255)
                  cACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCa                                                                                                                             <  5:230040/146‑1 (MQ=255)
                           gCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCt                                                                                                                 <  8:87652/149‑1 (MQ=255)
                            ccATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTg                                                                                                                >  4:23234/1‑149 (MQ=255)
                               tcatGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAg                                                                                                             <  6:31870/149‑1 (MQ=255)
                                   gCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCa                                                                                                                             >  6:136030/1‑129 (MQ=255)
                                   gCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCa                                                                                                                             <  5:136030/129‑1 (MQ=255)
                                         tATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTATTGCTtc                                                                                                   >  1:292917/1‑149 (MQ=255)
                                         tATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTtc                                                                                                   >  2:160997/1‑149 (MQ=255)
                                           tCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTtctc                                                                                                 <  8:382272/149‑1 (MQ=255)
                                             gCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTATTGCTTCTCa                                                                                                <  2:292917/148‑1 (MQ=255)
                                              cTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAgg                                                                                              >  2:355385/1‑149 (MQ=255)
                                                   ttttGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCTTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAgcg                                                                                         >  2:66912/1‑149 (MQ=255)
                                                   ttttGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCTTCACTCGCTGCAGGTAATGCTTCTCAGGTAgcg                                                                                         >  2:109296/1‑149 (MQ=255)
                                                         cGGACGGCTGTATATGCCGACTTTTTTTTTAGTTAAATCGATATCGGTACGCTATTCCAGGAAGCCACTGGCGCTTATTGCTTCTTCAACCGTCGCGGCCTCCTGCAGCGTCTCTCGCTGCAGGTAGTGCTTCTCAGGTAGCGCTTaa                                                                                    <  1:358817/148‑1 (MQ=255)
                                                          ggACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGc                                                                                  >  8:62890/1‑149 (MQ=255)
                                                             cGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTAGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCAc                                                                               <  1:355385/149‑1 (MQ=255)
                                                              ggCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTTAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCa                                                                                                <  4:177074/131‑1 (MQ=255)
                                                              ggCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTTAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCa                                                                                                >  3:177074/1‑131 (MQ=255)
                                                                 tGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTc                                                                            >  5:342657/1‑148 (MQ=255)
                                                                   tAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTTAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCaa                                                                          <  5:327616/148‑1 (MQ=255)
                                                                   aaaaTGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAAc                                                                         >  8:150634/2‑149 (MQ=255)
                                                                        gCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCaa                                                                    >  2:354613/1‑149 (MQ=255)
                                                                         ccGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCaa                                                                    <  8:156216/148‑1 (MQ=255)
                                                                          cGACTTTGTTTTTAGTTAAATCTATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAtt                                                                  >  4:311226/1‑149 (MQ=255)
                                                                          cGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAtt                                                                  >  7:309619/1‑149 (MQ=255)
                                                                          cGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAtt                                                                  >  5:313532/1‑149 (MQ=255)
                                                                           gACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAttt                                                                 <  6:313532/149‑1 (MQ=255)
                                                                           gACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAttt                                                                 >  3:157503/1‑149 (MQ=255)
                                                                           gACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAttt                                                                 >  3:153659/1‑149 (MQ=255)
                                                                            aCTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAAttt                                                                 >  8:290603/1‑148 (MQ=255)
                                                                             cTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTa                                                               <  1:354613/149‑1 (MQ=255)
                                                                             cTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTa                                                               >  3:218325/1‑149 (MQ=255)
                                                                             cTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACCGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTa                                                               <  4:198850/149‑1 (MQ=255)
                                                                              tttGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATaa                                                                                    >  5:225817/1‑127 (MQ=255)
                                                                              tttGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATaa                                                                                    <  6:225817/127‑1 (MQ=255)
                                                                              tttGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTAc                                                              >  1:700/1‑149 (MQ=255)
                                                                              tttGTTTTTAGTTAAATCGATATCGGTACGCAATTCAAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCGCTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTAc                                                              >  8:209296/1‑149 (MQ=255)
                                                                                tGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTAcc                                                             <  3:23234/148‑1 (MQ=255)
                                                                                 gTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTg                                                           >  5:21348/1‑149 (MQ=255)
                                                                                  tttttAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATaa                                                                                    >  4:295684/1‑123 (MQ=255)
                                                                                  tttttAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATaa                                                                                    <  3:295684/123‑1 (MQ=255)
                                                                                  tttttAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTg                                                           <  4:153659/148‑1 (MQ=255)
                                                                                       aGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGAGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTg                                                     >  2:230077/1‑149 (MQ=255)
                                                                                           aaaTCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTaaa                                                 >  7:94955/1‑149 (MQ=255)
                                                                                            aaTCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAAc                                                <  7:290603/149‑1 (MQ=255)
                                                                                             aTCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTgg                                                          <  7:99678/138‑1 (MQ=255)
                                                                                             aTCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTgg                                                          >  8:99678/1‑138 (MQ=255)
                                                                                                  tatCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCACCGGAAAGTTTACCTGGCACTGTAAACCTCttt                                          >  6:250831/1‑149 (MQ=255)
                                                                                                    tCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTc                                         <  3:311226/148‑1 (MQ=255)
                                                                                                    tCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTc                                         <  1:4720/148‑1 (MQ=255)
                                                                                                      ggTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCa                                      >  3:218570/1‑149 (MQ=255)
                                                                                                       gTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCAc                                     >  8:338354/1‑149 (MQ=255)
                                                                                                        tACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCttt                                          <  7:138127/143‑1 (MQ=255)
                                                                                                        tACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCttt                                          >  8:138127/1‑143 (MQ=255)
                                                                                                               ttCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTa                             <  8:309619/149‑1 (MQ=255)
                                                                                                                   aGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTaaaa                          <  5:375036/148‑1 (MQ=255)
                                                                                                                   aGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTaa                            >  6:375036/1‑146 (MQ=255)
                                                                                                                     cAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGcc                       >  3:67139/1‑149 (MQ=255)
                                                                                                                      aaGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCt                      <  7:150634/149‑1 (MQ=255)
                                                                                                                      aaGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCt                      >  8:328838/1‑149 (MQ=255)
                                                                                                                       aGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTg                     <  2:700/149‑1 (MQ=255)
                                                                                                                        gCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTg                     >  8:287492/1‑148 (MQ=255)
                                                                                                                        gCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTg                     <  7:287492/148‑1 (MQ=255)
                                                                                                                         ccACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCATGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGa                   <  6:342657/149‑1 (MQ=255)
                                                                                                                               gcgcGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCatat             <  6:21348/149‑1 (MQ=255)
                                                                                                                                 gcgAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCatat             >  8:219065/1‑147 (MQ=255)
                                                                                                                                 gcgAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCATATTa           <  7:219065/149‑1 (MQ=255)
                                                                                                                                          ttcttcAACCGTCGCGCCCTCCTTCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGACCCGTAAAAGCCTGGACCATATTAGCCGCCAGc  >  7:328145/1‑149 (MQ=255)
                                                                                                                                          ttcttcAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCATATTAGCCGCCAGc  >  4:107799/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GAGAGGGCGATAAATCTCCACCCGATCGCCATCATGCACACTATCGCTTAGTTTTGCCGGACGGCTGTAAATGCCGACTTTGTTTTTAGTTAAATCGATATCGGTACGCAATTCCAGCAAGCCACTGGCGCGAATAGCTTCTTCAACCGTCGCGCCCTCCTGCAGCGTCACTCGCTGCAGGTACTGCTTCTCAGGTAGCGCATAAGCCACCTCAACGGCAATTTTACCTGGCACTGTAAACCTCTTTCGCACGAACCGTAAAAGCCTGGACCATATTAGCCGCCAGC  >  NC_000913/2754065‑2754351

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