Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,481,515 T→C 10.6% S75G (AGC→GGC)  rhlB ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,481,5150TC10.6% 171.7 / 10.4 66S75G (AGC→GGC) rhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base C (5/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

TAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGT  >  NC_002947/1481367‑1481662
                                                                                                                                                      |                                                                                                                                                   
tAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCt                                                                                                                                                     >  6:231285/1‑149 (MQ=255)
    ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGcc                                                                                                                                                 >  7:314472/1‑149 (MQ=255)
       ggACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCtgt                                                                                                                                              >  4:137396/1‑149 (MQ=255)
         aCCTGCGCCTGGATCGGCGTGCAGTACGGGAATCCCAGGTCGTGTATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCtgtg                                                                                                                                             <  7:346669/148‑1 (MQ=255)
           cTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                                                  <  5:79476/141‑1 (MQ=255)
           cTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                                                  >  6:79476/1‑141 (MQ=255)
            tGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggctt                                                                                                                                         >  5:134986/1‑149 (MQ=255)
                        ggCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGtt                                                                                                                              <  8:314472/148‑1 (MQ=255)
                               aggaCGGTAAGCCCAGGTCGTGGAGGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttg                                                                                                                                        <  8:150994/128‑1 (MQ=255)
                               aGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttt                                                                                                                                        >  7:150994/1‑130 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggctt                                                                                                                                         >  8:287812/1‑128 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTgg                                                                                                                                            <  8:25110/125‑1 (MQ=255)
                                 tACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTgg                                                                                                                                            >  7:25110/1‑125 (MQ=255)
                                    gggAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGttt                                                                                                                 <  4:231932/149‑1 (MQ=255)
                                    gggAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGCGTATTGCCTTACTGCGGCTCGACCCCGAAGTCATCGGGCTTCCACAGGCTCGACTGTGGCTTCGGCTTGCGTTcgcgacgccacct                                                                                                                 >  8:92136/1‑144 (MQ=255)
                                       aaGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttg                                                                                                                                        >  3:339706/1‑123 (MQ=255)
                                       aaGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttg                                                                                                                                        <  4:339706/123‑1 (MQ=255)
                                          cccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCtt                                                                                                           <  8:133217/149‑1 (MQ=255)
                                           ccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACGGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGttt                                                                                                                 >  3:340730/1‑142 (MQ=255)
                                           ccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACGGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGttt                                                                                                                 <  4:340730/142‑1 (MQ=255)
                                           ccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTAGCGGAGAGGTTTGTCCTTg                                                                                                          >  3:72118/1‑149 (MQ=255)
                                                tCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCgg                                                                                                     <  2:337711/149‑1 (MQ=255)
                                                       tGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCg                                                                                              >  8:205891/1‑149 (MQ=255)
                                                        ggCGTGGATCAGCTGGTTGGACTGCTTGATGTCGTGGTAGCGGGTTTTGCCTTCCTGGGGCTGGACCACGAAGTGTTGGGGCTTCCACTGGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATc                                                                                               <  1:103208/148‑1 (MQ=255)
                                                          cctGCATAAGCTCGTAGGACAGCTTGAAGTCGTGGAAGCCGGTTAAGCCTACCTGAGGCTCGAACACTAACTCATCGCGCTTCCACAGGCTGGCCAGTGGCATGGGCTAGCGTTCGCGGCGCGGTTTCTCCTAGGCAGGTTCATCgc                                                                                             >  4:295024/3‑146 (MQ=255)
                                                               aTCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTgg                                                                                      >  7:66730/1‑149 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcgg                                                                              >  7:279351/1‑149 (MQ=255)
                                                                       gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcgg                                                                              <  7:205891/149‑1 (MQ=255)
                                                                           gACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGcggcgg                                                                           >  4:155936/1‑148 (MQ=255)
                                                                           gACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCggcggc                                                                          <  6:310770/149‑1 (MQ=255)
                                                                              aGCTTGAAGTCGTGGAAGTGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAgg                                                                       >  2:137089/1‑149 (MQ=255)
                                                                              aGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAgg                                                                       <  5:231285/149‑1 (MQ=255)
                                                                                cTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGtt                                                                     >  3:210199/1‑149 (MQ=255)
                                                                                   gAAGTCGTGGAAGTGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTg                                                                    <  1:137089/147‑1 (MQ=255)
                                                                                      gTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAg                                                               <  2:164334/149‑1 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCGTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGcc                                                             >  1:204555/1‑149 (MQ=255)
                                                                                         gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCt                                                            >  4:4271/1‑149 (MQ=255)
                                                                                           ggAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTg                                                          <  3:155936/149‑1 (MQ=255)
                                                                                             aaGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggg                                                        >  3:155555/1‑149 (MQ=255)
                                                                                               gCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggggg                                                      >  3:321485/1‑149 (MQ=255)
                                                                                                 gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCg                                                    >  2:270729/1‑149 (MQ=255)
                                                                                                       tGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGtt                                               <  3:295024/148‑1 (MQ=255)
                                                                                                        gCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTg                                              <  2:204555/148‑1 (MQ=255)
                                                                                                           ttCCTGCGGCTCGACCACGAAGTCTGCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCAt                                          >  5:289988/1‑149 (MQ=255)
                                                                                                             ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCt                                                            >  7:318809/1‑129 (MQ=255)
                                                                                                             ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCt                                                            <  8:318809/129‑1 (MQ=255)
                                                                                                              cTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg                                        <  8:279351/148‑1 (MQ=255)
                                                                                                                      tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                                                          >  6:34543/1‑75 (MQ=255)
                                                                                                                         aCCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTAATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggg                                                        <  6:147289/122‑1 (MQ=255)
                                                                                                                         aCCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTAATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTggg                                                        >  5:147289/1‑122 (MQ=255)
                                                                                                                           cACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtt                           >  5:38892/1‑149 (MQ=255)
                                                                                                                            aCGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGttt                          >  3:28459/1‑149 (MQ=255)
                                                                                                                                  tCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcgc                    >  7:277221/1‑149 (MQ=255)
                                                                                                                                    ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcg                     <  7:274788/146‑1 (MQ=255)
                                                                                                                                    ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcg                     >  8:274788/1‑146 (MQ=255)
                                                                                                                                        ggCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccg              <  6:134986/149‑1 (MQ=255)
                                                                                                                                            gccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTgg                                                         >  7:114036/2‑102 (MQ=255)
                                                                                                                                            gccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTgg                                                         <  8:114036/101‑1 (MQ=255)
                                                                                                                                            gccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCt          >  3:78607/2‑149 (MQ=255)
                                                                                                                                            gccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCt          >  2:123450/2‑149 (MQ=255)
                                                                                                                                             ccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCtt         >  5:150249/1‑149 (MQ=255)
                                                                                                                                              cacaGGCTGGCCTGTGGCTTGGGCTTGCGTTAGCGGCGCGGTTTGTCCTTGGCCGGTTGATAGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTc        >  7:90723/1‑149 (MQ=255)
                                                                                                                                                caGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg      <  2:322265/149‑1 (MQ=255)
                                                                                                                                                  aGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCgg     >  1:358254/1‑149 (MQ=255)
                                                                                                                                                   ggCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGACCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGc    >  1:259027/1‑149 (MQ=255)
                                                                                                                                                    gCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCg   <  4:28459/149‑1 (MQ=255)
                                                                                                                                                     cTGGGCGGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCCGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGt  <  8:90723/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TAACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGT  >  NC_002947/1481367‑1481662

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: