Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 198,777 T→A 6.3% L1428Q (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_002947198,7770TA 6.3% 261.5 / 5.5 90L1428Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (38/46);  new base A (3/3);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTG  >  NC_002947/198628‑198923
                                                                                                                                                     |                                                                                                                                                  
ggTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCt                                                                                                                                                    <  2:706795/150‑1 (MQ=255)
  tCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTgg                                                                                                                                                  <  1:857991/150‑1 (MQ=255)
    aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGTCTTCACTTCGGTAGCCGTTATCGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                                <  1:870930/150‑1 (MQ=18)
    aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                                <  2:318941/150‑1 (MQ=255)
    aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                                <  1:1726508/150‑1 (MQ=255)
    aGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTc                                                                                                                                                <  1:960447/150‑1 (MQ=255)
     gCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCg                                                                                                                                               <  2:2254570/150‑1 (MQ=255)
      cACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGt                                                                                                                                              <  2:1089996/150‑1 (MQ=255)
      cACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGt                                                                                                                                              <  2:1800021/150‑1 (MQ=255)
        ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGa                                                                                                                                            >  1:48178/1‑150 (MQ=255)
        ccGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATACCGGAAGAAACGCTATCTGTAGCCGACAACTTCAGACTGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGa                                                                                                                                            <  1:1451952/150‑1 (MQ=18)
         cGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGAc                                                                                                                                           <  2:48178/150‑1 (MQ=255)
         cGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTCGtcag                                                                                                                                           <  1:2335054/150‑4 (MQ=255)
         cGAGGTCACCGACGAGCCAGGTACGCCAAGCAACGAAGGCGATATCGTCAAGGTCACCATCACCGCCGACCAGACTTCGGTCGCCGAGAACGGCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCATGACCTAGTCATGAc                                                                                                                                           >  1:703300/1‑150 (MQ=17)
          gAGGTCACCTACGAGCCAGGTACGCCAGGCAACGAGGGCGATAGGGTCAAGGTCACGATCACGCCGCAGCCGACTTCGGGATCCGAGAACGTCAACACGACCTTCACTGTGCACATCACCACCGCCGTGCCTCACGACTTGGTCgccacc                                                                                                                                          >  2:295842/1‑145 (MQ=255)
           aGGTCACCGACGAGCCAGGTACGCCAGGAAACGAAGACGAACTCGTCTAGCACACCAGCACGGCCGACCCGACTTGGGGAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGAccc                                                                                                                                         >  2:1433957/1‑150 (MQ=17)
             gTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTg                                                                                                                                       >  2:411558/1‑150 (MQ=255)
              tCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                      <  1:800539/150‑1 (MQ=255)
                      gaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                      <  2:1245271/142‑1 (MQ=255)
                      gaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                      >  1:1245271/1‑142 (MQ=255)
                      gaGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACACCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                              <  1:1949515/150‑1 (MQ=255)
                            ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                        <  1:790181/150‑1 (MQ=255)
                            ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                        <  1:685089/150‑1 (MQ=255)
                            ggTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                        <  2:2227897/150‑1 (MQ=255)
                                  ccAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAcc                                                                                                                  <  2:1648325/150‑1 (MQ=255)
                                   cAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCa                                                                                                                 <  1:777354/150‑1 (MQ=255)
                                     ggCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATc                                                                                                               >  2:78878/1‑150 (MQ=255)
                                        aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                            <  1:2076017/150‑1 (MQ=255)
                                        aaCGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                            >  1:1879312/1‑150 (MQ=255)
                                            aaGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCg                                                                                                        >  2:1211355/1‑150 (MQ=255)
                                                 gATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGTCCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                              <  1:1895043/123‑1 (MQ=255)
                                                 gATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                              >  2:1895043/1‑123 (MQ=255)
                                                  atatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACg                                                                                                                             <  2:1733634/123‑1 (MQ=255)
                                                  atatCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACg                                                                                                                             >  1:1733634/1‑123 (MQ=255)
                                                     tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACAGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCa                                                                                               <  1:1495989/150‑1 (MQ=255)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGGTCACGACCAGGTCGTGAgcca                                                                                                                                        <  2:1105892/107‑2 (MQ=17)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGGTCACGACCAGGTCGTGAgcca                                                                                                                                        >  1:1105892/1‑106 (MQ=17)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                   >  2:1970123/1‑128 (MQ=255)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                   <  1:1970123/128‑1 (MQ=255)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                             <  2:411309/150‑1 (MQ=255)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                             <  1:411558/150‑1 (MQ=255)
                                                              tCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGt                                                                                      >  2:2097492/1‑150 (MQ=255)
                                                                 ccATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTaca                                                                                   >  1:1986884/1‑150 (MQ=255)
                                                                 ccATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTaca                                                                                   >  1:242686/1‑150 (MQ=255)
                                                                     cACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCa                                                                               <  2:929435/150‑1 (MQ=255)
                                                                        ggCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCacga                                                                            <  1:314117/150‑1 (MQ=255)
                                                                         gCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCacgac                                                                           <  1:2097492/150‑1 (MQ=255)
                                                                              ccAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCa                                                                      <  2:1879312/150‑1 (MQ=255)
                                                                               cAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACAGCCCTGGCTAACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAg                                                                     >  2:451930/1‑150 (MQ=255)
                                                                                 gACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                   >  2:846654/1‑150 (MQ=255)
                                                                                 gACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                   >  1:1071871/1‑150 (MQ=255)
                                                                                    ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                            <  2:961683/106‑1 (MQ=255)
                                                                                    ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                            >  1:961683/1‑106 (MQ=255)
                                                                                              gagaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCt                                                                                                                                                    <  2:296537/56‑1 (MQ=255)
                                                                                              gagaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCt                                                                                                                                                    >  1:296537/1‑56 (MQ=255)
                                                                                               agaACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATc                                                     >  1:1706389/1‑150 (MQ=255)
                                                                                                    gTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAg                                                                                              >  2:388065/1‑104 (MQ=255)
                                                                                                    gTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAg                                                                                              <  1:388065/104‑1 (MQ=255)
                                                                                                     tCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                 <  2:884545/68‑1 (MQ=255)
                                                                                                     tCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                 >  1:884545/1‑68 (MQ=255)
                                                                                                     tCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATg                                               >  1:168473/1‑150 (MQ=255)
                                                                                                      cAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTc                                                         >  1:1907300/1‑139 (MQ=255)
                                                                                                      cAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTc                                                         <  2:1907300/139‑1 (MQ=255)
                                                                                                        aaCCGACCTTCCCTCCGCACATCAACACCCCCCTTGCTCACGACCTGGTCTTGCCACTGAGCAACAACGCTCAGGTCACCAGCAAGGACGGCGAAACCAGCACACCGTACACTCACGACGCGCAGGGCGCTGACGTCTATCAGGATGCgg                                            <  2:242686/150‑1 (MQ=17)
                                                                                                        aaCCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                            >  1:767788/1‑150 (MQ=255)
                                                                                                         aCCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCggg                                           >  2:1005048/1‑150 (MQ=255)
                                                                                                          ccGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGc                                          <  1:451930/150‑1 (MQ=255)
                                                                                                          ccGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGc                                          <  1:132269/150‑1 (MQ=255)
                                                                                                            gACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTc                                                                                                                                                >  1:1490865/1‑46 (MQ=255)
                                                                                                            gACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCAGGTc                                                                                                                                                <  2:1490865/46‑1 (MQ=255)
                                                                                                               cTTCACTGTGCACATCAACACCGCCCTGGCTCACGACGTGGTCGTGACCCTGAGAAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGAt                                     >  1:179342/1‑150 (MQ=255)
                                                                                                               cTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGAt                                     >  2:623251/1‑150 (MQ=255)
                                                                                                                  cACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGAcgc                                                                         <  1:961721/111‑1 (MQ=255)
                                                                                                                  cACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGAcgc                                                                         >  2:961721/1‑111 (MQ=255)
                                                                                                                     tgtgCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACCACGCGCAGGGCCATGACGTCGATCAGTATTCGGGCCAGATCAGCCt                               >  2:1191388/1‑150 (MQ=18)
                                                                                                                      gtgCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGTTCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTg                              <  1:1191388/150‑1 (MQ=255)
                                                                                                                      gtgCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTg                              <  1:623251/150‑1 (MQ=255)
                                                                                                                            aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                            <  1:2059522/130‑1 (MQ=255)
                                                                                                                            aTCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCgg                                            >  2:2059522/1‑130 (MQ=255)
                                                                                                                               aaCACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAAc                     >  2:886168/1‑150 (MQ=255)
                                                                                                                                 cacCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAcc                                                                                                                  <  1:720935/55‑1 (MQ=255)
                                                                                                                                 cacCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAcc                                                                                                                  >  2:720935/1‑55 (MQ=255)
                                                                                                                                 cacCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                   >  1:2273488/1‑102 (MQ=255)
                                                                                                                                 cacCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAggg                                                                   <  2:2273488/102‑1 (MQ=255)
                                                                                                                                    cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCg                                             >  2:480912/1‑121 (MQ=255)
                                                                                                                                    cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCg                                             <  1:480912/121‑1 (MQ=255)
                                                                                                                                    cGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTcggc                >  2:1639672/1‑150 (MQ=255)
                                                                                                                                       tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGt             >  2:1616267/7‑150 (MQ=255)
                                                                                                                                        cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACCCTCAGGTGACCATCAAGCGCGGCGAAACCAGCGCACCGTACACCCACGACGCGCAGGGCGATGACTTCTATCAGGATGCGTGCCAGATCAGCCTGGGCATCAACTCGGCGGTg            <  2:846668/150‑1 (MQ=14)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                             >  1:2254708/1‑60 (MQ=38)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                             <  2:2254708/60‑1 (MQ=38)
                                                                                                                                                  aCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTg  <  2:1986884/150‑1 (MQ=255)
                                                                                                                                                  aCCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTg  >  1:458969/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTCAGCACCGAGGTCACCGACGAGCCAGGTACGCCAGGCAACGAAGGCGATATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACATCAACACCGCCCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGCGAAACCAGCGCACCGTACACTCACGACGCGCAGGGCGATGACGTCTATCAGGATGCGGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTG  >  NC_002947/198628‑198923

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