Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,148,802 A→C 6.0% L150R (CTT→CGT)  PP_1007 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,148,8020AC 6.0% 240.0 / 5.2 83L150R (CTT→CGT) PP_1007transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base C (2/3);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

TCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCACCT  >  NC_002947/1148655‑1148949
                                                                                                                                                   |                                                                                                                                                   
tCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAgg                                                                                                                                                   >  1:656606/1‑150 (MQ=255)
  ggggACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtg                                                                                                                                                 <  1:912804/150‑1 (MQ=255)
    ggACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAggtggt                                                                                                                                               <  2:1554917/150‑1 (MQ=255)
     gACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                              <  1:1680665/150‑1 (MQ=255)
     gACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGAGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                              <  1:1253433/150‑1 (MQ=255)
      aCCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCTTACGGTGGTCg                                                                                                                                             <  1:926163/150‑1 (MQ=255)
      aCCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCTTACGGTGGTCg                                                                                                                                             <  1:177701/150‑1 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGc                                                                                                                                           >  2:904313/1‑150 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGc                                                                                                                                           >  2:1125959/1‑150 (MQ=255)
        cTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGc                                                                                                                                           >  1:2257971/1‑150 (MQ=255)
          ggAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCTTACGGTGGTCGGcga                                                                                                                                         <  1:425564/150‑3 (MQ=255)
           gAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAg                                                                                                                                                    <  2:68789/138‑1 (MQ=255)
           gAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAg                                                                                                                                                    >  1:68789/1‑138 (MQ=255)
             aGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCa                                                                                                                                      >  2:969897/1‑150 (MQ=255)
              gTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAg                                                                                                                                     <  1:1684227/150‑1 (MQ=255)
                aaGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGAAACCGGTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgc                                                                                                                                   >  1:643796/1‑150 (MQ=255)
                  gCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGgcgc                                                                                                                                 >  1:1952836/1‑150 (MQ=255)
                    cTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCAt                                                                                                                               <  2:2203756/150‑1 (MQ=255)
                    cTCGCCCTGTAGCAGGCGGGCGAATTAGTGTCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCAt                                                                                                                               <  1:951/150‑1 (MQ=255)
                     tCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATg                                                                                                                              <  1:634466/150‑1 (MQ=255)
                          cTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAg                                                                                                                         <  2:967212/150‑1 (MQ=255)
                           tGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAgg                                                                                                                        >  2:1170683/1‑150 (MQ=255)
                            gTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGt                                                                                                                       >  2:1101288/1‑150 (MQ=255)
                                   gcgggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGAc                                                                                                                >  1:1167923/1‑150 (MQ=255)
                                      ggcgACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGt                                                                                                             >  1:2200606/1‑150 (MQ=255)
                                             tGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGAAGTTTTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGc                                                                                                      <  1:744808/150‑1 (MQ=255)
                                             gggcggcCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGAGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCCGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGc                                                                                                      <  2:1484216/149‑1 (MQ=255)
                                                   ccGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGAcc                                                                                                >  1:320451/1‑150 (MQ=255)
                                                     gTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcag                                                                                              >  1:1266565/1‑150 (MQ=255)
                                                     gTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcag                                                                                              >  2:2106809/1‑150 (MQ=255)
                                                      tCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagc                                                                                             <  1:558206/150‑1 (MQ=255)
                                                        cTGGCGTTCGCTGGAGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagcag                                                                                           <  1:1697193/150‑1 (MQ=255)
                                                            cGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTATTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCagcagcagc                                                                                       <  1:2231060/150‑1 (MQ=255)
                                                            cGTTCGCTGGCCAGTGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCTCCAGCACTCAATTACAATCTCCTACAACCGTAAGGTGGTCGGACTGAAGGCGCATGGGCTGAATACCGACGGTGAATAGCCCGACCagcagcagc                                                                                       <  1:1125959/150‑1 (MQ=255)
                                                                   tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCgg                                                                                >  2:407443/1‑150 (MQ=255)
                                                                   tGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCgg                                                                                >  1:872668/1‑150 (MQ=255)
                                                                        aaTGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGcc                                                                           >  2:1825724/1‑150 (MQ=255)
                                                                           gCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                              <  2:805568/80‑1 (MQ=255)
                                                                           gCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTc                                                                                                                                              >  1:805568/1‑80 (MQ=255)
                                                                           gCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAgg                                                                        <  2:643796/150‑1 (MQ=255)
                                                                            ccGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                       <  2:952207/150‑1 (MQ=255)
                                                                            ccGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAggg                                                                       <  1:1101288/150‑1 (MQ=255)
                                                                               aaTTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGtt                                                                    <  2:2257971/150‑1 (MQ=255)
                                                                               aaTTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGtt                                                                    <  2:1000166/150‑1 (MQ=255)
                                                                                 ttGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTc                                                                  >  2:1563693/1‑150 (MQ=255)
                                                                                    gtgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAg                                                               <  1:1170683/150‑1 (MQ=255)
                                                                                     tgtTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGCGGCAGAACTGCAGGCGCATGGGAAGAACAGCGACGGGGAATAGCCCGACCAGCAGCCGCGCTGCGGTGGCCAGGGGTTCCCAGt                                                              >  2:324368/1‑150 (MQ=255)
                                                                                          agcagcACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTg                                                         <  1:1715056/150‑1 (MQ=255)
                                                                                            cagcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTgcg                                                       >  2:286674/1‑150 (MQ=255)
                                                                                              gcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCa                                                     >  1:233224/1‑150 (MQ=255)
                                                                                              gcaCTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCa                                                     <  2:2200606/150‑1 (MQ=255)
                                                                                                 cTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCAACGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAAGAGCCCGACCAGCAGCAGCACTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgc                                                  <  2:1167923/150‑1 (MQ=255)
                                                                                                  ttttGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGgcg                                                 >  1:185133/1‑150 (MQ=255)
                                                                                                       cgcCGTCGTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCcgc                                            <  1:969897/150‑1 (MQ=255)
                                                                                                              ttCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACt                                                          >  2:1761312/1‑129 (MQ=255)
                                                                                                              ttCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACt                                                          <  1:1761312/129‑1 (MQ=255)
                                                                                                              ttCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGGGCCCGCGCCAt                                       >  2:1296825/1‑148 (MQ=255)
                                                                                                              ttCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCAt                                       <  1:1296825/148‑1 (MQ=255)
                                                                                                              ttCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGa                                     >  1:2227238/1‑150 (MQ=255)
                                                                                                                       cAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGTCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtg                            <  2:872668/150‑1 (MQ=255)
                                                                                                                        aGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc                           >  1:527380/1‑150 (MQ=255)
                                                                                                                        aGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc                           >  1:2280739/1‑150 (MQ=255)
                                                                                                                        aGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc                           >  2:207615/1‑150 (MQ=255)
                                                                                                                        aGGCGCTGGCGCTCGCCGACAACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTTCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgc                           >  1:1086826/1‑150 (MQ=255)
                                                                                                                         ggcgctggcgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgct                          >  2:1561478/1‑150 (MQ=255)
                                                                                                                         ggcgctggcgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgct                          >  2:153837/1‑150 (MQ=255)
                                                                                                                         ggcgctggcgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTtgct                          >  1:1897197/1‑150 (MQ=255)
                                                                                                                              tggcgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTgg                     <  2:233224/150‑1 (MQ=255)
                                                                                                                                gcgctCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGCCCTGCGTTGctgctgctg                   >  2:543459/1‑150 (MQ=255)
                                                                                                                                    tCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtg               >  1:1427412/1‑150 (MQ=255)
                                                                                                                                     cGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtgt              <  1:1825724/150‑1 (MQ=255)
                                                                                                                                      gCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtgtg             >  2:1375463/1‑150 (MQ=255)
                                                                                                                                      gCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCAtgtg             >  1:601655/1‑150 (MQ=255)
                                                                                                                                       ccGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGc            >  2:1802423/1‑150 (MQ=255)
                                                                                                                                        cGACCACCTTACGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCg           >  2:2251729/1‑150 (MQ=255)
                                                                                                                                         gACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTCCCAGTGACTGCGCAGGCGCCCGCGCCATGACCt                                  >  2:2053327/1‑126 (MQ=255)
                                                                                                                                         gACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTCCCAGTGACTGCGCAGGCGCCCGCGCCATGACCt                                  <  1:2053327/126‑1 (MQ=255)
                                                                                                                                          accaccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTg         <  1:1755693/150‑1 (MQ=255)
                                                                                                                                            caccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGc       >  2:1025699/1‑150 (MQ=255)
                                                                                                                                            caccGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGc       >  2:1786247/1‑150 (MQ=255)
                                                                                                                                               cGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagcagcag                                                                                        >  1:2369272/1‑66 (MQ=255)
                                                                                                                                               cGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACcagcagcag                                                                                        <  2:2369272/66‑1 (MQ=255)
                                                                                                                                                 tAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCAcag  <  1:207615/150‑3 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGGGGACCTGGAAGTAAGCCTCGCCCTGTAGCAGGCGGGCGACCTGGCGGCCGTCCTGGCGTTCGCTGGCGAATGCCGAATTGGTGTTGAGCAGCACTTTTGCGCCGTCTTCCAGTTGCAGGCGCTGGCGCTCGCCGACCACCGTAAGGTGGTCGGCCTGCAGGCGCATGGGCAGGTTACCGACGGTGAATAGCCCGACCAGCAGCAGCGCTGCGGTGGCCAGGGGTTTCCAGTGACTGCGCAGGCGCCCGCGCCATGACCTGCGTTGCTGCTGGTGCATGTGCGTGGCCACCT  >  NC_002947/1148655‑1148949

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