Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,726,458 T→G 10.4% F356L (TTT→TTG PP_3291 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,726,4580TG10.4% 226.0 / 18.3 87F356L (TTT→TTGPP_3291metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base G (6/3);  total (41/46)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

GGACGGTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATG  >  NC_002947/3726310‑3726601
                                                                                                                                                    |                                                                                                                                                
ggACGGTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGGTCTCCGGCTTTc                                                                                                                                                 >  1:390245/1‑150 (MQ=255)
   cGGTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAg                                                                                                                                              <  1:1597854/150‑1 (MQ=255)
    ggTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGCCGCCACTGGTCCCTAACCCGCGCAACAGTAGTCGTCGCTCCGGCTTTCAAGc                                                                                                                                             >  2:1647071/1‑150 (MQ=255)
    ggTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGc                                                                                                                                             <  1:637690/150‑1 (MQ=255)
       cGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCAAGCGGAATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGCCACTGGCCCCAAACCCGGGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgg                                                                                                                                          <  2:1162893/150‑1 (MQ=255)
         aTGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggg                                                                                                                                        <  2:1373298/150‑1 (MQ=255)
          tGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggggg                                                                                                                                       <  2:390245/150‑1 (MQ=255)
          ggAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggggg                                                                                                                                       <  2:775164/149‑1 (MQ=255)
           tggTCCAGGTGCCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACTGGTGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggggg                                                                                                                                      >  1:541933/3‑150 (MQ=255)
           gAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggggg                                                                                                                                      >  1:122251/1‑150 (MQ=255)
             gTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgg                                                                                                                                          >  1:427741/1‑144 (MQ=255)
             gTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgg                                                                                                                                          <  2:427741/144‑1 (MQ=255)
             gTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGAc                                                                                                                                    >  1:763023/1‑150 (MQ=255)
                cAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTGAAAGCCggag                                                                                                                                        <  2:1606121/143‑3 (MQ=255)
                cAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTGAAAGCCggag                                                                                                                                        >  1:1606121/1‑141 (MQ=255)
                    tgccTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACTGGTGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGt                                                                                                                             <  2:541933/148‑1 (MQ=255)
                    cACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGt                                                                                                                             <  2:271196/150‑1 (MQ=255)
                       cTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGcc                                                                                                                          >  2:1249457/1‑150 (MQ=255)
                       cTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGcc                                                                                                                          >  1:1684915/1‑150 (MQ=255)
                           aCCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGata                                                                                                                      <  1:1651874/150‑1 (MQ=255)
                           aCCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGata                                                                                                                      <  1:1886649/150‑1 (MQ=255)
                              aGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATAtcg                                                                                                                   <  2:724414/150‑1 (MQ=255)
                                  gcgcGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTcgc                                                                                                               <  2:1325756/150‑1 (MQ=255)
                                    gcgAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTATCGATATCGTCGCgg                                                                                                             <  1:1499419/150‑1 (MQ=255)
                                         cccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGcc                                                                                                        >  2:578239/1‑150 (MQ=255)
                                          ccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGccc                                                                                                       <  1:1321772/150‑1 (MQ=255)
                                            gCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCcac                                                                                                     >  1:1044055/1‑150 (MQ=255)
                                             cAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCcaca                                                                                                    >  2:1462840/1‑150 (MQ=255)
                                                tGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGGCGCGGGGCACCGATATCGTCGCGGGCGCCCACAgtg                                                                                                 >  1:1435642/1‑150 (MQ=255)
                                                  atcatcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGtgtg                                                                                               >  1:516376/1‑150 (MQ=255)
                                                  atcatcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGtgtg                                                                                               <  2:1111677/150‑1 (MQ=255)
                                                   tcatcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTATGATACGGGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGc                                                                                              <  2:731646/150‑1 (MQ=255)
                                                    ccatgccgctGGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCg                                                                                             >  2:669328/5‑150 (MQ=255)
                                                      tcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTTAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGt                                                                                           <  1:1198861/150‑1 (MQ=255)
                                                       cgccgcCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGtt                                                                                          <  2:1684915/150‑1 (MQ=255)
                                                       cgccgcCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGtt                                                                                          <  1:195781/150‑1 (MQ=255)
                                                                cGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTGCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGa                                                                                <  1:1382047/150‑1 (MQ=255)
                                                                cGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGa                                                                                <  1:1380321/150‑1 (MQ=255)
                                                                      ggCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGAtgt                                                                          <  1:461134/150‑1 (MQ=255)
                                                                         gACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGcc                                                                       >  1:1506206/1‑150 (MQ=255)
                                                                              gcggACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATcc                                                                  >  1:154280/1‑150 (MQ=255)
                                                                                ggACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGc                                                                        >  1:711143/1‑142 (MQ=255)
                                                                                ggACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGc                                                                        <  2:711143/142‑1 (MQ=255)
                                                                                  aCGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGc                                                              <  2:122251/150‑1 (MQ=255)
                                                                                    gggTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGcac                                                            >  1:1191260/1‑150 (MQ=255)
                                                                                        gCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGt                                                        <  2:763023/150‑1 (MQ=255)
                                                                                        gCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGt                                                        <  1:774956/150‑1 (MQ=255)
                                                                                        gCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGt                                                        <  1:345494/150‑1 (MQ=255)
                                                                                        acttacacaccGGAATGCACTGGACACAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGt                                                        >  2:550867/8‑150 (MQ=255)
                                                                                          ttCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGc                                                                                  <  2:1530575/122‑1 (MQ=255)
                                                                                          ttCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGc                                                                                  >  1:1530575/1‑122 (MQ=255)
                                                                                          ttCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGcaca                                                           >  2:742259/1‑145 (MQ=255)
                                                                                          ttCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGcaca                                                           <  1:742259/145‑1 (MQ=255)
                                                                                               caccTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCAAAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGAACAGGTGCATGCg                                                 <  2:1506206/150‑1 (MQ=255)
                                                                                                  cTGAAGGCACTGGCACCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATg                                              <  2:474291/150‑1 (MQ=255)
                                                                                                     aaGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGAAAGAGAATTGCCTATCCGTGCACAGGTGCATGCGATGgag                                           >  2:1096374/1‑150 (MQ=255)
                                                                                                       ggCACTGGCCCTAAACCCGCGCAATAGCATGCTGTTCTCCGGCTTTCAATCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGaa                                         <  2:607364/150‑1 (MQ=255)
                                                                                                       ggCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGaa                                         >  2:1074744/1‑150 (MQ=255)
                                                                                                        gCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCg                                                                                             >  2:346076/1‑98 (MQ=255)
                                                                                                        gCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCg                                                                                             <  1:346076/98‑1 (MQ=255)
                                                                                                          aCTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCAGGGCGGGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCt                                      >  1:1152284/1‑150 (MQ=255)
                                                                                                              gCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCggg                                                                                                                               <  2:456317/58‑1 (MQ=255)
                                                                                                              gCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCggg                                                                                                                               >  1:456317/1‑58 (MQ=255)
                                                                                                                cccAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTgc                                >  2:1142490/1‑150 (MQ=255)
                                                                                                                  cAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGga                                            <  1:1238078/136‑1 (MQ=255)
                                                                                                                  cAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGga                                            >  2:1238078/1‑136 (MQ=255)
                                                                                                                       ccGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGcc                         <  1:1228228/150‑1 (MQ=255)
                                                                                                                       ccGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGcc                         >  2:1224564/1‑150 (MQ=255)
                                                                                                                         gcgcAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGa                       <  2:229005/150‑1 (MQ=255)
                                                                                                                         gcgcAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACACGGGCTACCGCTATCGTCGGGCGCGCCCACCGTCTGCTGAACCACGGCCAGTATGGGCCTAGCCGTGCACAGGTGCATGCGATGGAGCAACCTTCTGCGCACGCCGa                       >  2:688730/1‑150 (MQ=255)
                                                                                                                             aaCAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATg                                              <  1:1175879/123‑1 (MQ=255)
                                                                                                                             aaCAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATg                                              >  2:1175879/1‑123 (MQ=255)
                                                                                                                             aaCAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGCGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCt                   >  1:1670861/1‑150 (MQ=255)
                                                                                                                                      cTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATc          >  1:1528593/1‑150 (MQ=255)
                                                                                                                                        gTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCAt        <  2:1158045/150‑1 (MQ=255)
                                                                                                                                        gTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCAt        <  2:516376/150‑1 (MQ=255)
                                                                                                                                          tctcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTACGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATgg      >  1:499341/1‑150 (MQ=255)
                                                                                                                                           ccccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAgtgt                                                                                                <  2:706281/58‑1 (MQ=38)
                                                                                                                                           ccccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAgtgt                                                                                                >  1:706281/3‑60 (MQ=38)
                                                                                                                                           ccccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGa     >  2:423753/3‑150 (MQ=255)
                                                                                                                                            tcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAAAAAGGCTGCGCACGCCGATGCTGATGAAATCATGGaa    <  1:1096374/150‑1 (MQ=255)
                                                                                                                                            cccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGa                                                <  1:1303075/105‑1 (MQ=255)
                                                                                                                                            cccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGa                                                >  2:1303075/2‑106 (MQ=255)
                                                                                                                                            cccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGaa    >  1:1407713/2‑150 (MQ=255)
                                                                                                                                             ccGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAAt   <  1:213538/150‑1 (MQ=255)
                                                                                                                                             ccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAAt   >  2:1556792/1‑150 (MQ=255)
                                                                                                                                              cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATg  >  1:461819/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGACGGTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATG  >  NC_002947/3726310‑3726601

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