Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 937,766 A→T 10.6% L2620Q (CTG→CAG)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947937,7660AT10.6% 190.6 / 14.8 67L2620Q (CTG→CAG) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base T (2/5);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GCAGAACCGCTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGA  >  NC_002947/937619‑937910
                                                                                                                                                   |                                                                                                                                                
gCAGAACCGCTGATGCTGCTGCCATCGGCAGTGACGACCAGTTCGCTCGGGCTGGTCGGCGCGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGgcagc                                                                                                                                                <  2:1483211/150‑1 (MQ=255)
          tGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggt                                                                                                                                      <  2:1541592/150‑1 (MQ=255)
           gATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtg                                                                                                                                     <  2:284781/150‑2 (MQ=255)
              gctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggtgaac                                                                                                                                  <  2:397624/150‑5 (MQ=255)
              gctgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgag                                                                                                                                  <  1:1443510/150‑1 (MQ=255)
                tgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                <  2:916861/150‑1 (MQ=255)
                tgctgcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggtcgaggt                                                                                                                                >  1:803346/1‑150 (MQ=255)
                    gcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGAt                                                                                                                            >  1:344559/1‑150 (MQ=255)
                    gcCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGAt                                                                                                                            <  1:1556983/150‑1 (MQ=255)
                     ccATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATc                                                                                                                           >  2:97829/1‑150 (MQ=255)
                        tCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGTCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAgg                                                                                                                        >  1:609626/1‑150 (MQ=255)
                        tCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGCCGTACGTTAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAgg                                                                                                                        >  2:682650/1‑150 (MQ=255)
                         cGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGa                                                                                                                       <  2:609626/150‑1 (MQ=255)
                                  cgacCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGt                                                                                                              >  1:1050466/1‑150 (MQ=255)
                                     ccAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTcg                                                                                                           <  2:803346/150‑1 (MQ=255)
                                      cAGTTCGCTCGGGCTGGTCGGCGTGACCCCCAGTGGCACGTCGTAGGTCAGCGATTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGACATTCGCAACGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTcgc                                                                                                          <  1:1518980/150‑1 (MQ=255)
                                            gCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCAc                                                                                                    <  1:1470457/150‑1 (MQ=255)
                                                 ggCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGc                                                                                               <  1:1111133/150‑1 (MQ=255)
                                                 ggCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGc                                                                                               >  1:318796/1‑150 (MQ=255)
                                                   cTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCg                                                                                             <  1:97829/150‑1 (MQ=255)
                                                       tCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTAc                                                                                         >  1:1099841/1‑150 (MQ=255)
                                                         ggCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggt                                                                                                                                      <  1:128636/103‑1 (MQ=255)
                                                         ggCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggt                                                                                                                                      >  2:128636/1‑103 (MQ=255)
                                                           cGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGAc                                                                                     >  1:1109508/1‑150 (MQ=255)
                                                                ccaccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCCGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACGGACCAGCg                                                                                >  2:728371/1‑150 (MQ=255)
                                                                 caccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGt                                                                               >  1:1592830/1‑150 (MQ=255)
                                                                  accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTTGAGCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTg                                                                              >  1:1177543/1‑150 (MQ=255)
                                                                  accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGaaa                                                                                                                     <  2:730167/111‑1 (MQ=255)
                                                                  accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGaaa                                                                                                                     >  1:730167/1‑111 (MQ=255)
                                                                  accaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGCCTCAGGCGTTAACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCCACCAGCGTg                                                                              >  2:683288/1‑150 (MQ=255)
                                                                   ccaGTGGCACGTCGTATGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGcc                                                                                                                                          <  1:1018610/89‑1 (MQ=255)
                                                                   ccaGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGcc                                                                                                                                          >  2:1018610/1‑89 (MQ=255)
                                                                         gCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTgg                                                                       <  1:655394/150‑1 (MQ=255)
                                                                           aCGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCa                                                                     >  2:494955/1‑150 (MQ=255)
                                                                           aCGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCa                                                                     >  2:1452648/1‑150 (MQ=255)
                                                                              tcgtAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGcc                                                                            <  1:1472262/140‑1 (MQ=255)
                                                                              tcgtAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGcc                                                                            >  2:1472262/1‑140 (MQ=255)
                                                                                gtAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATcgcg                                                                <  1:1065404/150‑1 (MQ=255)
                                                                                gtAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATcgcg                                                                >  2:1065404/1‑150 (MQ=255)
                                                                                     tCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCt                                                           >  2:1646800/1‑150 (MQ=255)
                                                                                     ccAGCGCTTCAGAGCCGTTGCCGCCCGGGTCGTTCTGGGCCGGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCt                                                           <  2:14086/149‑1 (MQ=255)
                                                                                       aGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGGAATTACCGACCAGCGTGCCATTGGCATCGCGCACCTCg                                                         >  1:1319897/1‑150 (MQ=255)
                                                                                                 agGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAg                                                                                  <  2:1175719/115‑1 (MQ=255)
                                                                                                 agGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAg                                                                                  >  1:1175719/1‑115 (MQ=255)
                                                                                                   gCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCg                                             <  1:1646800/150‑1 (MQ=255)
                                                                                                     gTTGCCGCTCGGGTCGGTCTTTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggt                                                                                                                                      <  2:1006135/59‑1 (MQ=37)
                                                                                                     gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggt                                                                                                                                      >  1:1006135/1‑59 (MQ=38)
                                                                                                     gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGt                                                    <  2:481929/141‑1 (MQ=255)
                                                                                                     gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGt                                                    >  1:481929/1‑141 (MQ=255)
                                                                                                        gCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTc                                        <  1:1527915/150‑1 (MQ=255)
                                                                                                           gCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCAGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTcgcc                                     >  1:1138218/1‑150 (MQ=255)
                                                                                                             tctggtcggtcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTcgccgc                                   <  1:682650/147‑1 (MQ=255)
                                                                                                                ggtcggtcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGAc                                <  1:494955/150‑1 (MQ=255)
                                                                                                                 gtcggtcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGAcc                               >  2:717391/1‑150 (MQ=255)
                                                                                                                      tcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGc                                           >  2:65877/1‑133 (MQ=255)
                                                                                                                      tcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGc                                           <  1:65877/133‑1 (MQ=255)
                                                                                                                            ccAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGAc                                                       >  1:568462/1‑115 (MQ=255)
                                                                                                                            ccAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGAc                                                       <  2:568462/115‑1 (MQ=255)
                                                                                                                             cAGGCTCAGGCGTTCACCCGGCCGCGCAACAAGGCCCTGGCCGTTCACGACAGTATCGTGCTCTTTCACTGTGCCGTTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTACACCCGGCTTCGCCGCGACCGCCCAAGGTCAg                   <  2:1177543/150‑1 (MQ=255)
                                                                                                                              aGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAgg                  <  2:318796/150‑1 (MQ=255)
                                                                                                                               ggCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGc                 <  1:717391/150‑1 (MQ=255)
                                                                                                                                      gCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGgc          <  1:1452648/150‑1 (MQ=255)
                                                                                                                                        gTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGgcgc        >  2:1448152/1‑150 (MQ=255)
                                                                                                                                           gaCCCGGCTGCGCTGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTa                                                   <  2:656278/103‑1 (MQ=255)
                                                                                                                                           gaCCCGGCTGCGCTGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTa                                                   >  1:656278/2‑104 (MQ=255)
                                                                                                                                             cccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacg   >  1:922616/1‑150 (MQ=255)
                                                                                                                                              ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacga  <  1:728371/150‑1 (MQ=255)
                                                                                                                                              ccGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCcgacga  >  2:849431/1‑150 (MQ=255)
                                                                                                                                                 gcagcgcagcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTc                                                          <  1:1101051/91‑1 (MQ=255)
                                                                                                                                                 gcagcgcagcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTc                                                          >  2:1101051/1‑91 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GCAGAACCGCTGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGA  >  NC_002947/937619‑937910

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: ‑