Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,768,318 C→A 10.4% P97Q (CCG→CAG)  PP_5600 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,768,3180CA10.4% 226.0 / 5.4 73P97Q (CCG→CAG) PP_5600hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base A (3/5);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCGGGCGTGTGCTCACC  >  NC_002947/3768169‑3768467
                                                                                                                                                     |                                                                                                                                                      
cACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGcc                                                                                                                                                        >  1:201410/1‑150 (MQ=255)
  cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                                      <  2:565246/150‑1 (MQ=255)
  cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAgg                                                                                                                                                      <  2:1439033/150‑1 (MQ=255)
  cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAgg                                                                                                                                                      <  1:592375/150‑1 (MQ=255)
  cGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAgg                                                                                                                                                      <  1:1723318/150‑1 (MQ=255)
      cGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAg                                                                                                                                                 >  2:2416586/1‑150 (MQ=255)
      cGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAg                                                                                                                                                 >  1:1919415/1‑150 (MQ=255)
      cGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAg                                                                                                                                                 >  2:276543/1‑150 (MQ=255)
       gACAGCAGCGGCATGGTCCAAGGTGAGTAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGc                                                                                                                                                >  1:1815795/1‑150 (MQ=255)
        aCAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGcc                                                                                                                                               >  1:1714167/1‑150 (MQ=255)
           gcagcGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCt                                                                                                                                            >  2:2159084/1‑150 (MQ=255)
            cagcGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCtt                                                                                                                                           <  1:2159084/150‑1 (MQ=255)
               cGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTAcc                                                                                                                                        <  2:2149407/150‑1 (MQ=255)
                 gCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACcca                                                                                                                                      <  2:201410/150‑1 (MQ=255)
                 gCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACcca                                                                                                                                      >  1:1029134/1‑150 (MQ=255)
                  cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaa                                                                                                                                     <  1:1826458/150‑1 (MQ=255)
                       tCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGcc                                                                                                                                                        >  1:1408344/1‑127 (MQ=255)
                       tCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGcc                                                                                                                                                        <  2:1408344/127‑1 (MQ=255)
                           aGGTGAGCAAGACTTCAGCGATCCGACCGTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                                      <  2:1305546/125‑1 (MQ=255)
                           aGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                                      >  1:1305546/1‑125 (MQ=255)
                           aGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTcc                                                                                                                            >  2:777832/1‑150 (MQ=255)
                           aGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTcc                                                                                                                            <  2:1919415/150‑1 (MQ=255)
                           aGGCGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                                      <  2:1299399/125‑1 (MQ=255)
                           aGGCGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                                      >  1:1299399/1‑125 (MQ=255)
                             gTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGAACCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGc                                                                                                                          >  1:1959428/1‑150 (MQ=255)
                                      aCTTCAGCGCTCCGCCCCGGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACTCAACCAACCTCCGCCACGCCTGc                                                                                                                 <  2:1029134/150‑1 (MQ=255)
                                          cAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCaa                                                                                                             >  2:299664/1‑150 (MQ=255)
                                           aGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAag                                                                                                            <  1:193122/150‑1 (MQ=255)
                                                tCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAgcg                                                                                                       >  2:2044681/1‑150 (MQ=255)
                                                  cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaa                                                                                                                                     >  2:1911993/1‑118 (MQ=255)
                                                  cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaa                                                                                                                                     <  1:1911993/118‑1 (MQ=255)
                                                  cGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCt                                                                                                     >  1:2075386/1‑150 (MQ=255)
                                                    aCCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTgg                                                                                                   <  1:276543/150‑1 (MQ=255)
                                                          tCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTg                                                                                             <  1:1281327/150‑1 (MQ=255)
                                                              ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagca                                                                                                          >  2:266860/1‑133 (MQ=255)
                                                              ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagca                                                                                                          <  1:266860/133‑1 (MQ=255)
                                                                    gTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATcccc                                                                                   <  1:2416586/150‑1 (MQ=255)
                                                                      gACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgc                                                                                 >  2:1636293/1‑150 (MQ=255)
                                                                      gACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgc                                                                                 <  2:2075386/150‑1 (MQ=255)
                                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATcccc                                                                                   >  2:1927367/1‑145 (MQ=255)
                                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATcccc                                                                                   <  1:1927367/145‑1 (MQ=255)
                                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGcc                                                                              >  1:2196139/1‑150 (MQ=255)
                                                                         tGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGcc                                                                              >  2:1580598/1‑150 (MQ=255)
                                                                               aTCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCGGCGAGCAgg                                                                                                                                                      <  2:1851319/73‑1 (MQ=255)
                                                                                tagCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcg                                                                       >  1:1910890/3‑150 (MQ=255)
                                                                                   ccAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgctt                                                                    <  1:777832/150‑1 (MQ=255)
                                                                                    cAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttc                                                                   >  1:1977974/1‑150 (MQ=255)
                                                                                    cAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTCTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGACTCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttc                                                                   <  1:156163/150‑1 (MQ=255)
                                                                                       ccTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgct                                                                >  1:1678848/1‑150 (MQ=255)
                                                                                        cTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgc                                                                 >  1:2418759/1‑148 (MQ=255)
                                                                                        cTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgc                                                                 <  2:2418759/148‑1 (MQ=255)
                                                                                         tGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCgcttcgcttc                                                              >  1:1954644/1‑150 (MQ=255)
                                                                                          gTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTtct                                                             >  1:911869/1‑150 (MQ=255)
                                                                                              aGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTcgcg                                                         >  1:1312716/1‑150 (MQ=255)
                                                                                                  ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCt                                                     >  1:482328/1‑150 (MQ=255)
                                                                                                    acacTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGa                                                   <  2:1714167/150‑1 (MQ=255)
                                                                                                    acacTGCTGGCCGTCGTTTCCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGGCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTgc                                                   <  2:912968/150‑2 (MQ=255)
                                                                                                      acTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACa                                                 >  1:443646/1‑150 (MQ=255)
                                                                                                      acTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCAGGGTCAGCCCCTTGCCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGAc                                                  >  1:1526834/1‑150 (MQ=255)
                                                                                                       ctgctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACAt                                                <  2:1815795/150‑1 (MQ=255)
                                                                                                        tgctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCATCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATc                                               <  1:2260299/150‑1 (MQ=255)
                                                                                                          ctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGc                                             <  1:1058065/150‑1 (MQ=255)
                                                                                                                 tCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTg                                                                                                    <  2:2174340/88‑1 (MQ=255)
                                                                                                                 tCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTg                                                                                                    >  1:2174340/1‑88 (MQ=255)
                                                                                                                      ttCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTg                                 <  2:1275867/150‑1 (MQ=255)
                                                                                                                          gCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTgccgc                             >  2:37967/1‑150 (MQ=255)
                                                                                                                          gCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTgccgc                             <  1:1580598/150‑1 (MQ=255)
                                                                                                                             tgctgACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTg                          >  1:1808423/1‑150 (MQ=255)
                                                                                                                                tgACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCTCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATc                       <  1:1636293/150‑1 (MQ=255)
                                                                                                                                 tgcccggCTCGCCGGCGAGCAGGGTCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATc                       >  1:2363099/7‑150 (MQ=255)
                                                                                                                                 gACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCg                      <  2:1910890/150‑1 (MQ=255)
                                                                                                                                  aCCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCgg                     <  2:1071129/150‑1 (MQ=255)
                                                                                                                                  aCCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCgg                     >  2:2053714/1‑150 (MQ=255)
                                                                                                                                         tCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCggg              >  1:860376/1‑150 (MQ=255)
                                                                                                                                          cGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCGGGc             <  1:2053714/150‑1 (MQ=255)
                                                                                                                                                 gAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCGGGCGTGTGCt      >  1:204581/1‑150 (MQ=255)
                                                                                                                                                     cGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCGGGCGTGTGCTCAct  <  2:1678848/150‑2 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGGCAGCGGGCGTGTGCTCACC  >  NC_002947/3768169‑3768467

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