Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,481,961 T→A 7.0% Q408L (CAG→CTG)  PP_3090 ← OmpA domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,481,9610TA 7.0% 572.5 / 22.7 182Q408L (CAG→CTG) PP_3090OmpA domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (84/85);  new base A (6/7);  total (90/92)
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.79e-01

Only 100 of 194 total aligned reads displayed.

TGCACAGCACCAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAGGGTTG  >  NC_002947/3481812‑3482110
                                                                                                                                                      |                                                                                                                                                     
tGCACAGCACCAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCggct                                                                                                                                                       <  1:3446519/150‑1 (MQ=255)
   acaGCACCAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCggctggc                                                                                                                                                    >  2:1269607/1‑150 (MQ=255)
        aCCAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCg                                                                                                                                               <  1:3412164/150‑1 (MQ=255)
          cAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGc                                                                                                                                             >  1:2379617/1‑150 (MQ=255)
           aGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCATATGCTGGTCGGGAATGATCTCTTTGCGGAACTGGCTGTTGATGAACAGCGGCGCCGGCAGCCAGCCGGcg                                                                                                                                            <  2:3101518/150‑2 (MQ=255)
           aGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAGTGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCAGCCAGCCGGcg                                                                                                                                            <  1:3658630/150‑2 (MQ=255)
           aGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCAGCCAGCCGGcg                                                                                                                                            <  2:3762016/150‑2 (MQ=255)
            gTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCAGCCAGCCGGcgc                                                                                                                                           <  2:716967/150‑3 (MQ=255)
             tGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGt                                                                                                                                          <  1:3848351/150‑1 (MQ=255)
                cgccAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGcc                                                                                                                                       >  1:482660/1‑150 (MQ=255)
                 gccAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGccc                                                                                                                                      >  1:1582344/1‑150 (MQ=255)
                     ggCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTAc                                                                                                                                  <  2:2489176/150‑1 (MQ=255)
                     ggCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTAc                                                                                                                                  <  1:4598958/150‑1 (MQ=255)
                      gCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTAcc                                                                                                                                 <  2:5277537/150‑1 (MQ=255)
                        cggcggCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAACTGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCtt                                                                                                                               <  1:4575914/150‑1 (MQ=255)
                        cggcggCGCCCAACCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGGTGGTGTACAGCGCCACCATATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGACGGCTGGCTGCCGGCATTGCCCGTACCtt                                                                                                                               <  1:1465459/150‑1 (MQ=255)
                          gcggcgCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCt                                                                                                                             >  1:74064/1‑150 (MQ=255)
                          gcggcgCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGGGCCCGTACCTTCt                                                                                                                             <  2:4042923/150‑1 (MQ=255)
                            ggcgCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGc                                                                                                                           <  1:49012/150‑1 (MQ=255)
                            ggcgCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGc                                                                                                                           <  2:3072595/150‑1 (MQ=255)
                              cgcCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCgt                                                                                                                         <  1:2619661/150‑1 (MQ=255)
                               gcCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGATGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCgtg                                                                                                                        <  2:1232907/150‑1 (MQ=255)
                                cccacccaCCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCgtgt                                                                                                                       >  2:2768241/1‑150 (MQ=255)
                                cccacccaCCCGGCCGTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCgtgt                                                                                                                       <  2:443971/150‑1 (MQ=255)
                                 ccacccacCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCgtgt                                                                                                                       <  2:3550677/150‑1 (MQ=255)
                                  cacccaccCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTc                                                                                                                     >  2:1618507/1‑150 (MQ=255)
                                   acccacccGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCg                                                                                                                    >  2:5024629/1‑150 (MQ=255)
                                    cccacccGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGa                                                                                                                   <  1:877382/150‑1 (MQ=255)
                                     ccacccGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAg                                                                                                                  >  2:4867889/1‑150 (MQ=255)
                                     ccacccGGCCTTGCGCCGGCGTTCACGGTGGTTGCAGGTGTACAGAGCCACCAAATGCTGGACGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAg                                                                                                                  >  2:5348797/1‑150 (MQ=255)
                                       acccGGCCGTGCGCCGGCGTTCACGGGGGTTGCGGGGGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTGTGCGGAACAGGCTGGGGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAgcg                                                                                                                >  2:3615906/1‑150 (MQ=255)
                                             ccTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGc                                                                                                          <  1:10173/150‑1 (MQ=255)
                                             ccTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGc                                                                                                          <  2:1510493/150‑1 (MQ=255)
                                               ttGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGt                                                                                                        >  2:3448388/1‑150 (MQ=255)
                                                 gcgcCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTc                                                                                                      <  1:2352463/150‑1 (MQ=255)
                                                       gCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCcac                                                                                                <  2:3684851/150‑1 (MQ=255)
                                                         gTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACAt                                                                                              >  1:2503828/1‑150 (MQ=255)
                                                             aCGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTg                                                                                          >  2:1708410/1‑150 (MQ=255)
                                                                 tggtTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCGGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCg                                                                                      <  2:963771/150‑1 (MQ=255)
                                                                  ggtTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGt                                                                                     <  1:1708410/150‑1 (MQ=255)
                                                                  ggtTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGt                                                                                     <  2:3921210/150‑1 (MQ=255)
                                                                     tGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGccagcat                                                                                  <  2:3483039/150‑6 (MQ=255)
                                                                       cTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAg                                                                                >  1:4082249/1‑150 (MQ=255)
                                                                        tGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGc                                                                               <  2:543518/150‑1 (MQ=255)
                                                                           tgtACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGcccc                                                                            >  1:4690506/1‑150 (MQ=255)
                                                                             tACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCg                                                                                                     >  2:2007021/1‑123 (MQ=255)
                                                                             tACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCg                                                                                                     <  1:2007021/123‑1 (MQ=255)
                                                                              aCAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACt                                                                         >  1:1839169/1‑150 (MQ=255)
                                                                                aGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGa                                                                       >  2:2902793/1‑150 (MQ=255)
                                                                                      accaAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcg                                                                 >  1:2406663/1‑150 (MQ=255)
                                                                                        caAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgtc                                                               >  2:5389033/1‑150 (MQ=255)
                                                                                        caAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgtc                                                               <  2:476677/150‑1 (MQ=255)
                                                                                        caAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTACCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgtc                                                               >  2:4655099/1‑150 (MQ=255)
                                                                                          aaaTGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgtcg                                                              >  1:3092066/1‑150 (MQ=255)
                                                                                          aaaTGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgtcg                                                              >  2:781401/1‑150 (MQ=255)
                                                                                            aTGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGc                                                            <  1:83617/150‑1 (MQ=255)
                                                                                              gCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCt                                                          <  1:4679444/150‑1 (MQ=255)
                                                                                                     gggAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGcc                                                                                                  >  1:2601194/1‑103 (MQ=255)
                                                                                                     gggAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGcc                                                                                                  <  2:2601194/103‑1 (MQ=255)
                                                                                                        aaTGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGcgccgc                                                >  2:1097251/1‑150 (MQ=255)
                                                                                                          tGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGc                                                                                                                                        >  1:4254690/1‑60 (MQ=255)
                                                                                                          tGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGc                                                                                                                                        <  2:4254690/60‑1 (MQ=255)
                                                                                                           gATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGt                                             <  2:2955030/150‑1 (MQ=255)
                                                                                                           gATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGt                                             >  1:5120220/1‑150 (MQ=255)
                                                                                                            aTCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTg                                            <  1:91933/150‑1 (MQ=255)
                                                                                                             tCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGa                                           <  2:5321089/150‑1 (MQ=255)
                                                                                                              cACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGaa                                          >  2:1154735/1‑150 (MQ=255)
                                                                                                               aCTTTGCGGAACAGGCTGTTGATGGACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGaaa                                         <  1:3542205/150‑1 (MQ=255)
                                                                                                               aCTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGaaa                                         <  1:2179243/150‑1 (MQ=255)
                                                                                                               aCTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGaaa                                         <  1:5287351/150‑1 (MQ=255)
                                                                                                                 tttGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATa                                       <  2:4591784/150‑1 (MQ=255)
                                                                                                                 tttGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATa                                       >  1:3065554/1‑150 (MQ=255)
                                                                                                                 tttGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATa                                       <  1:1069136/150‑1 (MQ=255)
                                                                                                                 ttagCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATa                                       >  2:3201452/4‑150 (MQ=255)
                                                                                                                   tGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGa                                     >  1:5029642/1‑150 (MQ=255)
                                                                                                                     cGGAACAGGCTGTGGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGGACCTTCTGCGTGTTCGAGCGCGAAGCGGTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGACGGCCTGCAGCGCCGCGGTGAAATCGaaa                                   <  1:2717468/150‑1 (MQ=255)
                                                                                                                      ggAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAAc                                  >  1:4441901/1‑150 (MQ=255)
                                                                                                                      ggAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAAc                                  <  2:3331870/150‑1 (MQ=255)
                                                                                                                         aCAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACcgc                               <  1:5388960/150‑1 (MQ=255)
                                                                                                                               tGTTGCTGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCCCATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCgcagca                         >  2:3108540/1‑150 (MQ=255)
                                                                                                                               tGTAGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCgcagca                         >  1:381318/1‑150 (MQ=255)
                                                                                                                                 ttGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCAtt                       <  1:1800653/150‑1 (MQ=255)
                                                                                                                                  tgatgaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCAttt                      <  2:4567042/150‑1 (MQ=255)
                                                                                                                                   gatgaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTc                     >  2:870277/1‑150 (MQ=255)
                                                                                                                                   gatgaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTc                     >  2:1843118/1‑150 (MQ=255)
                                                                                                                                   gatgaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTc                     >  2:580035/1‑150 (MQ=255)
                                                                                                                                     tgaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCgg                   >  1:5294837/1‑150 (MQ=255)
                                                                                                                                      gaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCtcgt                                                                <  2:4052123/104‑1 (MQ=255)
                                                                                                                                      gaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGc                  <  2:497493/150‑1 (MQ=255)
                                                                                                                                      gaACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGc                  >  1:4353690/1‑150 (MQ=255)
                                                                                                                                           gcactgcCGGCAGCCAGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCt             >  2:1350457/6‑150 (MQ=255)
                                                                                                                                           gcactgcCGGCAGCCAGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCt             >  1:1257258/6‑150 (MQ=255)
                                                                                                                                             actgcCGGCAGCCAGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTgg           >  2:1462268/4‑150 (MQ=255)
                                                                                                                                               tgcCGGCAGCCAGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGGGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTa         >  2:2311746/2‑150 (MQ=255)
                                                                                                                                                gcCGGCTGGCTGCCGGCAGTGCCCGGACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTaa        <  2:3577762/150‑1 (MQ=255)
                                                                                                                                                  cGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAgg      >  2:1193580/1‑150 (MQ=255)
                                                                                                                                                  cGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAgg      <  1:1132317/150‑1 (MQ=255)
                                                                                                                                                    gctggctgCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAGGGt    <  2:3150773/150‑1 (MQ=255)
                                                                                                                                                      tggctgCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAGGGTTg  <  1:2240278/150‑1 (MQ=255)
                                                                                                                                                      tggctgCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAGGGTTg  <  2:22304/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGCACAGCACCAGTGCCGCCAGGCCGGCGGCGCCCACCCACCCGGCCTTGCGCCGGCGTTCACGGTGGTTGCTGGTGTACAGCGCCACCAAATGCTGGTCGGGAATGATCACTTTGCGGAACAGGCTGTTGATGAACAGCGGCGCCGGCTGGCTGCCGGCAGTGCCCGTACCTTCTGCGTGTTCGAGCGCGAAGCGTTCGGCCACATGCTGGCCGTGGCAGCCCCACTGAGCCTCGTCGGCCTGCAGCGCCGCGGTGAAATAGAAACCGCGCAGCATTTCGGCATTCTGGTAAGGGTTG  >  NC_002947/3481812‑3482110

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