Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,484,651 C→T 5.8% R531C (CGC→TGC)  hrpA → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,484,6510CT 5.8% 220.9 / 2.8 69R531C (CGC→TGC) hrpAputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCGCCTATAAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGC  >  NC_000913/1484505‑1484797
                                                                                                                                                  |                                                                                                                                                  
gcgcCTATAAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCcgcgc                                                                                                                                                  >  4:17352/1‑149 (MQ=255)
  gcCTATAAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGa                                                                                                                                                <  8:7423/149‑1 (MQ=255)
   ccTATAAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGaa                                                                                                                                               >  2:71919/1‑149 (MQ=255)
     tataAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACg                                                                                                                                             >  2:95136/1‑149 (MQ=255)
     tataAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACg                                                                                                                                             >  6:42336/1‑149 (MQ=255)
         aaCTGAGGCGGCCCGGCCGCCAGCCCCCGCAGTTCCCGGGCGACCCACGGCGGGCGGGTAGGGGGCGGGAAGCGAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGa                                                                                                                                        <  1:110481/149‑1 (MQ=255)
          aCTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGCCCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCg                                                                                                                                         >  8:141698/1‑148 (MQ=255)
            tGACGCCGCTCGGTCGCCAGCTCTCGCAGATGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATg                                                                                                                                      >  4:172202/1‑149 (MQ=255)
               cgccgcTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGAc                                                                                                                                   <  2:37180/149‑1 (MQ=255)
                gccgcTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGAc                                                                                                                                   >  5:68424/1‑148 (MQ=255)
                gccgcTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGAc                                                                                                                                   >  5:356242/1‑148 (MQ=255)
                   gcgcGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAAc                                                                                                                               <  5:287950/146‑1 (MQ=255)
                   gcTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACaaa                                                                                                                                >  7:336207/1‑148 (MQ=255)
                        gTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTTCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAAcagcag                                                                                                                          <  5:121019/149‑1 (MQ=255)
                        gTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAAcagcag                                                                                                                          <  6:293380/149‑1 (MQ=255)
                        gTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAAcagcag                                                                                                                          <  1:71919/149‑1 (MQ=255)
                           gCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGGGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCa                                                                                                                       >  5:257137/1‑149 (MQ=255)
                                   tcGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGAc                                                                                                                 <  7:211898/147‑1 (MQ=255)
                                        gTTGCCTGTCGACACACGTCTGGCGCGTATGGGGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCAGGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCTGACGAAAAAc                                                                                                          <  6:171975/149‑1 (MQ=255)
                                         ttGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACa                                                                                                         >  6:264891/1‑149 (MQ=255)
                                         ttGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACa                                                                                                         <  1:291613/149‑1 (MQ=255)
                                              tGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGTGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACAtcgtc                                                                                                    >  1:45810/1‑149 (MQ=255)
                                                   aCCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGa                                                                                                                                        >  7:99336/1‑108 (MQ=255)
                                                   aCCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGa                                                                                                                                        <  8:99336/108‑1 (MQ=255)
                                                    cccACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTc                                                                                               >  3:103822/1‑148 (MQ=255)
                                                     ccACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCa                                                                                             >  1:90931/1‑149 (MQ=255)
                                                     ccACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCa                                                                                             >  2:187822/1‑149 (MQ=255)
                                                      cACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCAc                                                                                            <  3:48225/149‑1 (MQ=255)
                                                         gTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGAc                                                                                         <  2:31287/149‑1 (MQ=255)
                                                              gcgcGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAaga                                                                                    <  6:257137/149‑1 (MQ=255)
                                                              gcgcGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAag                                                                                     <  2:90931/148‑1 (MQ=255)
                                                                  gTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCt                                                                                <  5:185530/149‑1 (MQ=255)
                                                                   tATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCt                                                                                >  6:254655/1‑148 (MQ=255)
                                                                     tGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTgcc                                                                             >  5:253136/1‑147 (MQ=255)
                                                                     tGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGAc                                                                             >  8:33461/1‑149 (MQ=255)
                                                                         gCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGTGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTtc                                                                         <  2:45810/149‑1 (MQ=255)
                                                                           tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTtctc                                                                       >  6:73851/1‑149 (MQ=255)
                                                                           tGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTtctc                                                                       >  6:187162/1‑149 (MQ=255)
                                                                               aGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGt                                                                   <  4:67431/149‑1 (MQ=255)
                                                                                gcgcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGtt                                                                  >  3:41258/1‑149 (MQ=255)
                                                                                  gcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTtg                                                                >  5:69702/1‑149 (MQ=255)
                                                                                  gcAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTtg                                                                >  6:11970/1‑149 (MQ=255)
                                                                                     aaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAaga                                                                                    >  8:278738/1‑126 (MQ=255)
                                                                                     aaaaCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAaga                                                                                    <  7:278738/126‑1 (MQ=255)
                                                                                        aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAAt                                                           >  6:307079/1‑148 (MQ=255)
                                                                                        aCATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAAt                                                           <  5:307079/148‑1 (MQ=255)
                                                                                           tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCATGATCCGCGCGTACGTCCGATGGACAAACAGCTGGCATCGGACGAAATACTTCGTCGCTTCGACGACAAAGAGTCTGTCTTTCTCGCGTTTGTGAATCtgt                                                       <  6:125212/149‑1 (MQ=255)
                                                                                           tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCTCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCCTCGATTCCACGGCAAAGAGGATGCCTTTCTCTCGTTTGTGAATCtgt                                                       >  5:147222/1‑149 (MQ=255)
                                                                                           tGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCtgt                                                       <  8:336207/149‑1 (MQ=255)
                                                                                             gCTGCGTGCGTGAGGCGATTATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCtgtg                                                      <  5:222636/148‑1 (MQ=255)
                                                                                             gCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTgg                                                     <  6:356242/149‑1 (MQ=255)
                                                                                             gCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTgg                                                     <  6:68424/149‑1 (MQ=255)
                                                                                                      gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGt                                                                                  >  1:151678/1‑111 (MQ=255)
                                                                                                      gtgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGt                                                                                  <  2:151678/111‑1 (MQ=255)
                                                                                                       tgAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTg                                           <  1:187822/149‑1 (MQ=255)
                                                                                                         aGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTgg                                          <  5:254655/148‑1 (MQ=255)
                                                                                                           gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGTGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCtgtg                                                      <  3:287823/134‑1 (MQ=255)
                                                                                                           gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGTGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCtgtg                                                      >  4:287823/1‑134 (MQ=255)
                                                                                                           gCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGa                                       <  3:17352/149‑1 (MQ=255)
                                                                                                              atgatTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGag                                      >  3:260599/1‑147 (MQ=255)
                                                                                                                     tCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCaaaa                              <  1:95136/148‑1 (MQ=255)
                                                                                                                       aCGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgc                           >  6:308099/1‑149 (MQ=255)
                                                                                                                       aCGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGgc                           >  2:173211/1‑149 (MQ=255)
                                                                                                                           ccGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCtt                       >  3:191796/1‑149 (MQ=255)
                                                                                                                               gcTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTttctt                   <  6:69702/149‑1 (MQ=255)
                                                                                                                                     cATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                    <  6:54420/126‑1 (MQ=255)
                                                                                                                                     cATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGagca                                    >  5:54420/1‑126 (MQ=255)
                                                                                                                                       ttCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGcc           >  3:250325/1‑149 (MQ=255)
                                                                                                                                          aGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCtt         <  6:253136/148‑1 (MQ=255)
                                                                                                                                            gATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCg      >  2:18050/1‑149 (MQ=255)
                                                                                                                                                 gcgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGc  <  7:205511/148‑1 (MQ=255)
                                                                                                                                                 gcgcgAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGc  >  8:205511/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGCCTATAAACTGACGCCGCTCGGTCGCCAGCTCTCGCAGTTGCCTGTCGACCCACGTCTGGCGCGTATGGTGCTGGAAGCGCAAAAACATGGCTGCGTGCGTGAGGCGATGATTATCACGTCCGCGCTCTCCATTCAGGATCCGCGCGAACGTCCGATGGACAAACAGCAGGCATCGGACGAAAAACATCGTCGCTTCCACGACAAAGAGTCTGACTTTCTCGCGTTTGTGAATCTGTGGAATTATCTTGGCGAGCAGCAAAAGGCGCTTTCTTCCAACGCCTTCCGTCGC  >  NC_000913/1484505‑1484797

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