Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 206,064 T→A 6.4% L3857Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947206,0640TA 6.4% 222.0 / 6.1 78L3857Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base A (2/3);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

TCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTC  >  NC_002947/205917‑206213
                                                                                                                                                   |                                                                                                                                                     
tCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                     <  2:1170116/150‑1 (MQ=255)
 tggCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGt                                                                                                                                                    >  1:1733767/3‑150 (MQ=255)
  aGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   >  2:1362959/1‑150 (MQ=255)
  aGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  2:692751/150‑1 (MQ=255)
    ctcGTCGGTTACCGACGTGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 <  1:740782/148‑1 (MQ=255)
    cACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 <  2:685497/150‑1 (MQ=39)
      ggtctggtACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                               >  2:378105/8‑150 (MQ=255)
      cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTTGTCGTGa                                                                                                                                               <  2:842065/150‑1 (MQ=255)
      cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                               >  2:411507/1‑150 (MQ=40)
      cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                               <  2:1733767/150‑1 (MQ=255)
      agtaggatACCAACGAGCCAGGTACGCCGAGAAACGAAGGCGACCTGGTCAAAGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                               <  2:898625/143‑1 (MQ=14)
       gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGGCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACCACCTGGTCGTGAc                                                                                                                                              >  1:856829/1‑150 (MQ=255)
            ttACCGACGAGCCGGGTACGCGGGGCAACGATGGCGACCTGGTCAAGGTCACCATCACAGCCGACCATCTTCCGGGGCCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                         <  1:118980/150‑1 (MQ=14)
                 gacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    <  1:1726980/150‑1 (MQ=255)
                       ccAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   >  1:1000618/1‑129 (MQ=255)
                       ccAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                              >  1:1733227/1‑150 (MQ=255)
                       ccAGGTACGCCGGGCAACGAAGGCGACCTCGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  2:1000618/129‑1 (MQ=255)
                              cGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAGCCCCACTCTTGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCa                                                                                                                       >  2:990346/1‑150 (MQ=14)
                                       aCGAAGGCGACCTGGTCAAGGTCACAAGCTCTTCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGCTCTTCCCCGTGCACGTAAACCAGCCGCTGCCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCACgg                                                                                                              >  1:528149/1‑150 (MQ=11)
                                        cGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                             >  1:1775585/1‑150 (MQ=255)
                                        cGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                             >  2:1027520/1‑150 (MQ=255)
                                        cGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGCGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGc                                                                                                             >  1:1779361/1‑150 (MQ=255)
                                         gAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGcc                                                                                                            >  1:511197/1‑150 (MQ=40)
                                            ggCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGc                                                                                                         <  2:1335344/150‑1 (MQ=37)
                                             gCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACAAAGGTGGCCGAGACCGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCg                                                                                                        >  2:544389/1‑150 (MQ=255)
                                                aCCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                               >  1:62818/1‑140 (MQ=255)
                                                aCCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                               <  2:62818/140‑1 (MQ=255)
                                                 ccTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGAc                                                                                                    <  2:1733227/150‑1 (MQ=255)
                                                  cTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGAcc                                                                                                   >  1:484910/1‑150 (MQ=255)
                                                   tGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCa                                                                                                  <  2:1775585/150‑1 (MQ=255)
                                                     gTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgc                                                                                                >  1:452312/1‑150 (MQ=255)
                                                      tCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcg                                                                                               >  2:6453/1‑150 (MQ=255)
                                                      tCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcg                                                                                               >  1:1399563/1‑150 (MQ=255)
                                                        aaGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcgcg                                                                                             >  2:422888/1‑150 (MQ=255)
                                                         aGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAgcgcgc                                                                                            <  2:952242/150‑1 (MQ=255)
                                                              aCCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGcc                                                                                                            >  1:972272/1‑129 (MQ=255)
                                                              aCCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGcc                                                                                                            <  2:972272/129‑1 (MQ=255)
                                                               ccATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATaca                                                                                      <  2:1399563/150‑1 (MQ=255)
                                                                cATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATacac                                                                                     >  2:1289792/1‑150 (MQ=255)
                                                                        ccGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGGACGTAAACCAGCCGCTGGCGCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTg                                                                             >  2:1468260/1‑150 (MQ=255)
                                                                           aCCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAAGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTgcgc                                                                          >  2:976809/1‑150 (MQ=255)
                                                                                aCTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGc                                                                     >  2:594997/1‑150 (MQ=255)
                                                                                   tCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                                  <  2:603927/150‑1 (MQ=255)
                                                                                     ggTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGa                                                                >  2:1229087/1‑150 (MQ=255)
                                                                                         gCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCGCACGACCTGGTCGTGACCCTGCGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                            <  1:602944/150‑1 (MQ=255)
                                                                                           cGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                          <  1:6453/150‑1 (MQ=255)
                                                                                                      aaCCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTg                                               <  2:511197/150‑1 (MQ=255)
                                                                                                      aaCCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTg                                               >  2:1903796/1‑150 (MQ=255)
                                                                                                       aCCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTgg                                              >  2:378199/1‑150 (MQ=255)
                                                                                                        ccGATCTTCACCGCGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGc                                             <  2:1779361/150‑1 (MQ=255)
                                                                                                         cGATCTTCACCGGGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGcc                                            >  2:292353/1‑150 (MQ=255)
                                                                                                          gATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCa                                           <  1:1229087/150‑1 (MQ=255)
                                                                                                                   cGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGTTCACCATCAACGCCGCCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACCTCTACAACGACGCTGGCCAGATCAGCCt                                  <  2:259527/150‑1 (MQ=255)
                                                                                                                    gTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTg                                 <  2:499983/150‑1 (MQ=38)
                                                                                                                     tGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTgg                                >  1:590059/1‑150 (MQ=38)
                                                                                                                       cACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGc                              >  2:811257/1‑150 (MQ=39)
                                                                                                                          gTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                               >  2:443749/1‑66 (MQ=35)
                                                                                                                          gTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAg                                                                                                               <  1:443749/66‑1 (MQ=35)
                                                                                                                          gTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATc                           >  2:724488/1‑150 (MQ=38)
                                                                                                                            aaaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCAc                         <  2:1429371/150‑1 (MQ=37)
                                                                                                                             aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGGGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACg                        >  2:1093113/1‑150 (MQ=25)
                                                                                                                                   ccGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTggg                               >  1:979221/1‑137 (MQ=37)
                                                                                                                                   ccGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCAGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTggg                               <  2:979221/137‑1 (MQ=25)
                                                                                                                                     gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCtgt                <  1:378199/150‑1 (MQ=38)
                                                                                                                                      cTGGCTCACGACCTGGTCGTGCCCCTGAACAACATCGCCCACCAGTCGATTTTCCTCTGCGAGACCAGCGCGCCATACCCCCTCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCtgtg               <  2:528149/150‑1 (MQ=2)
                                                                                                                                        ggCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCCGCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGa             <  1:1153799/150‑1 (MQ=255)
                                                                                                                                        ggCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGa             <  1:724488/150‑1 (MQ=38)
                                                                                                                                         ggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTg                                               >  1:522831/3‑115 (MQ=255)
                                                                                                                                         ggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTg                                               <  2:522831/113‑1 (MQ=255)
                                                                                                                                          cTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGCCGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATg           >  2:340788/1‑150 (MQ=25)
                                                                                                                                          cTCACGACCTGGTCGGGACCCTGAGCACCACCGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATg           <  2:856829/150‑1 (MQ=255)
                                                                                                                                           tCACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCtgtg               <  1:332222/145‑1 (MQ=37)
                                                                                                                                           tCACGACCAGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCtgtg               >  2:332222/1‑145 (MQ=25)
                                                                                                                                           ccACGACCTGGTCGTGACCCGGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGt          <  1:322033/149‑1 (MQ=25)
                                                                                                                                               gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGAc      >  1:656401/1‑150 (MQ=34)
                                                                                                                                                 ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTgg              <  2:1555789/140‑1 (MQ=37)
                                                                                                                                                 ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTgg              >  1:1555789/1‑140 (MQ=37)
                                                                                                                                                 ccTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATAATTACGCGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACgg    <  1:594997/150‑1 (MQ=11)
                                                                                                                                                   tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTc  >  2:1765231/1‑150 (MQ=35)
                                                                                                                                                   tGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTc  >  1:1430045/1‑150 (MQ=35)
                                                                                                                                                   |                                                                                                                                                     
TCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACAGCCGACCAGACTTCGGTGGCCGAGAACGTCAAACCGATCTTCACCGTGCACGTAAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAGACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGACGCTGGCCAGATCAGCCTGGGCATCACGTCTGCTGTGGATGTAGACGGTC  >  NC_002947/205917‑206213

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