Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 937,771 A→T 10.2% A2618A (GCT→GCA PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947937,7710AT10.2% 197.8 / 12.1 67A2618A (GCT→GCAPP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base A (31/29);  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 = 7.81e-01

TGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATGCTGCT  >  NC_002947/937629‑937920
                                                                                                                                              |                                                                                                                                                     
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                                                                                                                                                >  2:128636/1‑103 (MQ=255)
                                               ggCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGggt                                                                                                                                                <  1:128636/103‑1 (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                                                                                            >  1:1175719/1‑115 (MQ=255)
                                                                                       agGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAg                                                                                            <  2:1175719/115‑1 (MQ=255)
                                                                                         gCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCg                                                       <  1:1646800/150‑1 (MQ=255)
                                                                                           gTTGCCGCTCGGGTCGGTCTTTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggt                                                                                                                                                <  2:1006135/59‑1 (MQ=37)
                                                                                           gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCTGCGCTGCCGggt                                                                                                                                                >  1:1006135/1‑59 (MQ=38)
                                                                                           gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGt                                                              >  1:481929/1‑141 (MQ=255)
                                                                                           gTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGt                                                              <  2:481929/141‑1 (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                                                     <  1:65877/133‑1 (MQ=255)
                                                                                                            tcTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGc                                                     >  2:65877/1‑133 (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                                                             >  1:656278/2‑104 (MQ=255)
                                                                                                                                 gaCCCGGCTGCGCTGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTa                                                             <  2:656278/103‑1 (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)
                                                                                                                                              agcCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATgctgct  >  2:832288/1‑150 (MQ=255)
                                                                                                                                              agcCGAGTCGAGGTCGATCAGGAAAGGACCGTTCGCGCTCACTGTGCCGGTACCGACCAGCGTCCCATTCGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGTGACCGCCCAAGTTCAGGCCGTCGGCGCCGACGACAATgatgct  <  1:830313/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TGATGCTGCTGCCATCGGCAGCGACGACCAGTTCGCTCGGGCTGGTCGGCGTGACCACCAGTGGCACGTCGTAGGTCAGCGCTTCAGAGGCGTTGCCGCTCGGGTCGGTCTGTACCAGGCTCAGGCGTTCACCCGGCAGCGCAGCCGGGTCGAGGTCGATCAGGAAAGTACCGTTCGCGGTCACTGTGCCGGTACCGACCAGCGTGCCATTGGCATCGCGCACCTCGACGGTAGCACCGGCTTCGCCGCGACCGCCCAAGGTCAGGCCGTCGGCGCCGACGACAATGCTGCT  >  NC_002947/937629‑937920

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: