Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,902,717 T→C 7.1% V196A (GTG→GCG) ‡ PP_4310 → hydantoin racemase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,902,7170TC 7.1% 433.7 / 14.6 151V196A (GTG→GCG) ‡PP_4310hydantoin racemase
Reads supporting (aligned to +/- strand):  ref base T (62/77);  new base C (6/5);  total (69/82)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 152 total aligned reads displayed.

TGGAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTAGAGATCCT  >  NC_002947/4902570‑4902866
                                                                                                                                                   |                                                                                                                                                     
tggaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGc                                                                                                                                                     >  2:1859395/1‑150 (MQ=255)
tggaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGc                                                                                                                                                     >  1:651879/1‑150 (MQ=255)
 ggaactggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGcc                                                                                                                                                    >  1:1367722/1‑150 (MQ=255)
    tctggacGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTGGTGTCGTGGTGGCCGGCCTGGGCGAGCAGAGTCGCCCACGCACCGAGGTGCgggg                                                                                                                                                 >  2:1494670/2‑149 (MQ=255)
    actggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCg                                                                                                                                                   <  2:1841710/148‑1 (MQ=255)
    actggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCg                                                                                                                                                   >  1:1841710/1‑148 (MQ=255)
    actggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggt                                                                                                                                                 <  2:214623/150‑1 (MQ=255)
      tggaaGCCGACCCGCTGCGGGCCGTGGAGGAGATTGTCGAGCAGGCCGAACGCGCCGAGCGCGACAACAAGGCCGAAATGGACGGCCTGGTTTGTGGCGGGTTGGTCGGGCTGGACGAGGAGATTCGCCAACGCACCGGCGTGCCggtgg                                                                                                                                               >  2:347401/1‑150 (MQ=255)
       ggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtggt                                                                                                                                              <  1:1859395/150‑1 (MQ=255)
       ggaaGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtggt                                                                                                                                              <  2:1420149/150‑1 (MQ=255)
           gCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGAt                                                                                                                                          <  2:1420241/150‑1 (MQ=255)
           gCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGAt                                                                                                                                          <  1:611935/150‑1 (MQ=255)
            ccGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATg                                                                                                                                         <  1:204703/150‑1 (MQ=255)
             cGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggtgcgttg                                                                                                                                           <  2:977096/147‑6 (MQ=255)
             cGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggtgcgttg                                                                                                                                           >  1:977096/1‑142 (MQ=255)
                 ccGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCCTGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGtgtg                                                                                                                                    >  1:1496873/1‑150 (MQ=255)
                      tcgtgCCGTGAAGACGATTGCTGCGCAGGCCGAGCTCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTTAGTGc                                                                                                                               <  2:27795/146‑1 (MQ=255)
                        gggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAg                                                                                                                             >  1:2332075/1‑150 (MQ=255)
                        gggCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAg                                                                                                                             <  1:887496/150‑1 (MQ=255)
                          gCCGTGGAGGCGATTGTCGAGCAGTCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGGTGTGGCGGCATGGCCGGGCTGGACGAGCAGCTTCGCCAACGCACCGGCGGGCAGGTGGGGGATGGTGTGAGAGCAGCg                                                                                                                           >  2:1483685/1‑150 (MQ=255)
                          gCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCg                                                                                                                           <  1:1945187/150‑1 (MQ=255)
                            catGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGt                                                                                                                         <  2:2641552/148‑1 (MQ=255)
                              tGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGa                                                                                                                       <  1:1314199/150‑1 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGgtgt                                                                                                                                     >  1:129993/1‑129 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGgtgt                                                                                                                                     <  2:129993/129‑1 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGgtgt                                                                                                                                     >  2:377394/1‑129 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGgtgt                                                                                                                                     <  1:377394/129‑1 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGAc                                                                                                                      >  1:1625112/1‑144 (MQ=255)
                                     gATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGAc                                                                                                                      <  2:1625112/144‑1 (MQ=255)
                                      aTTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGcc                                                                                                               <  1:1239912/150‑1 (MQ=255)
                                         gTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGaa                                                                                                            <  1:1785306/150‑1 (MQ=255)
                                          tCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAAt                                                                                                           >  2:229268/1‑150 (MQ=255)
                                           cGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggt                                                                                                                                                 <  1:383009/111‑1 (MQ=255)
                                           cGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCACGCCggt                                                                                                                                                 >  2:383009/1‑111 (MQ=255)
                                              gCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtgcgtt                                                                                                                                            <  1:1370186/113‑5 (MQ=255)
                                              gCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCggtgcgtt                                                                                                                                            >  2:1370186/1‑109 (MQ=255)
                                                   ccGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATCGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgc                                                                                                  <  2:1044550/150‑1 (MQ=255)
                                                    cGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgcg                                                                                                 >  1:1511366/1‑150 (MQ=255)
                                                    cGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgcg                                                                                                 >  2:2320852/1‑150 (MQ=255)
                                                       aCGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCAt                                                                                              <  1:1937324/150‑1 (MQ=255)
                                                        cgcgCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATg                                                                                             <  1:523246/150‑1 (MQ=255)
                                                         gcgcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATgg                                                                                            <  1:1691293/150‑1 (MQ=255)
                                                         gcgcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATgg                                                                                            >  2:2960866/1‑150 (MQ=255)
                                                         gcgcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATgg                                                                                            <  2:2376490/150‑1 (MQ=255)
                                                           gcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCAtt                                                                                                                  >  2:1178292/1‑126 (MQ=255)
                                                           gcCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCAtt                                                                                                                  <  1:1178292/126‑1 (MQ=255)
                                                              gTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTg                                                                                       >  1:2372339/1‑150 (MQ=255)
                                                                gcgcGACGACAAGGCGGAGGTCATCCGCCGGGGTTGTGGCGGCATGGGCGGGCTGGACGAGCAGGTTCGCCAACGCACCGGCGCGCCGGTGGTGGATGGTGTCGGTGCAGCGGTGACGATTGCGGAATCGCTGGTGCGCATGGGTTTGAg                                                                                     <  1:1494670/150‑1 (MQ=255)
                                                                   cgacgaCAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGAAAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCAc                                                                                  <  2:730547/150‑1 (MQ=255)
                                                                       gacAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGttt                                                                                        <  1:1436995/140‑1 (MQ=255)
                                                                        acAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCAc                                                                                  >  1:248996/1‑145 (MQ=255)
                                                                          aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgcgc                                                                                                <  1:2193691/129‑1 (MQ=255)
                                                                          aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTgcgc                                                                                                >  2:2193691/1‑129 (MQ=255)
                                                                          aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                           >  1:1016185/1‑150 (MQ=255)
                                                                          aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                           <  2:1733510/150‑1 (MQ=255)
                                                                          aaGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                           >  1:2542837/1‑150 (MQ=255)
                                                                           aGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATAGCGGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAg                                                                          <  2:2253596/150‑1 (MQ=255)
                                                                            ggCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGt                                                                         <  1:173989/150‑1 (MQ=255)
                                                                            ggCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGt                                                                         <  1:229268/150‑1 (MQ=255)
                                                                                  ggTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGc                                                                                                        <  1:2167132/113‑1 (MQ=255)
                                                                                    tGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCt                                                                 <  2:1367722/150‑1 (MQ=255)
                                                                                      aTCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGa                                                                                      >  1:1539771/1‑127 (MQ=255)
                                                                                      aTCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGa                                                                                      <  2:1539771/127‑1 (MQ=255)
                                                                                       tCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATg                                                              <  1:1604726/150‑1 (MQ=255)
                                                                                         tGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGa                                                                                                                                   <  1:987539/79‑1 (MQ=255)
                                                                                          gCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGa                                                           <  2:2542837/150‑1 (MQ=255)
                                                                                          gCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGa                                                           <  2:1016185/150‑1 (MQ=255)
                                                                                             tGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGc                                                                                                                <  1:720962/94‑1 (MQ=255)
                                                                                             tGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCa                                                                                   >  2:331316/1‑123 (MQ=255)
                                                                                             tGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCa                                                                                   <  1:331316/123‑1 (MQ=255)
                                                                                                      gcggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                           <  1:2681833/122‑1 (MQ=255)
                                                                                                      gcggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCaaa                                                                           >  2:2681833/1‑122 (MQ=255)
                                                                                                      gcggcATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCga                                               >  1:2023057/1‑150 (MQ=255)
                                                                                                                cGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgtggt                                     >  1:2688475/1‑150 (MQ=255)
                                                                                                                  ggCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGtggtgg                                   <  1:816367/150‑1 (MQ=255)
                                                                                                                   gCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTggg                                  >  2:2708423/1‑150 (MQ=255)
                                                                                                                    caggACGAGCAGATTCGCCAACGCATCGGCGTGCCGGAGGTGGAGGGTGTGAGTGCAGCGGTCACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACCCCACGGGGGAAGAAAgtgggtggg                                 >  2:832252/3‑145 (MQ=255)
                                                                                                                      ggACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTg                               <  2:2118464/150‑1 (MQ=255)
                                                                                                                          gAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCg                           >  1:1879213/1‑150 (MQ=255)
                                                                                                                           ggccgATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGa                          <  2:1902257/146‑1 (MQ=255)
                                                                                                                           aGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGa                          <  2:2372339/150‑1 (MQ=255)
                                                                                                                            gCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGAt                         >  2:1705696/1‑150 (MQ=255)
                                                                                                                             cAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGGGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATg                        <  1:1441741/150‑1 (MQ=255)
                                                                                                                                aTTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGc                             <  2:2918233/142‑1 (MQ=255)
                                                                                                                                aTTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGc                             >  1:2918233/1‑142 (MQ=255)
                                                                                                                                   cGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTc                  <  1:347401/150‑1 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 <  1:2708423/150‑1 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 >  2:789373/1‑150 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 >  1:2316860/1‑150 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 >  1:855982/1‑150 (MQ=255)
                                                                                                                                    gCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGCTTGAATACGTGCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 >  1:2065588/1‑150 (MQ=255)
                                                                                                                                    catccaccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCg                 >  1:1006519/8‑150 (MQ=255)
                                                                                                                                      cAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGc               >  2:687834/1‑150 (MQ=255)
                                                                                                                                        aCGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCg             >  2:606608/1‑150 (MQ=255)
                                                                                                                                         accACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc            >  1:2246728/3‑150 (MQ=255)
                                                                                                                                         accACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGc            >  2:1999082/3‑150 (MQ=255)
                                                                                                                                          gCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGAAGCGCATCGGcca             >  2:2541841/1‑147 (MQ=255)
                                                                                                                                          gCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGAAGCGCATCGGcca             <  1:2541841/148‑2 (MQ=255)
                                                                                                                                          ccACCGGCACGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAgtggt                                     <  2:726076/123‑1 (MQ=255)
                                                                                                                                           cACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTa          <  2:2316860/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGGAACTGGAAGCCGACCCGCTGCGGGCCGTGGAGGCGATTGTCGAGCAGGCCGAACGCGCCGTGCGCGACGACAAGGCCGAGGTGATCTGCCTGGGTTGTGGCGGCATGGCCGGGCTGGACGAGCAGATTCGCCAACGCACCGGCGTGCCGGTGGTGGATGGTGTGAGTGCAGCGGTGACCATTGCCGAATCGCTGGTGCGCATGGGTTTGAGCACGTCCAAAGTGCGCACCTATGCGACGCCACGGGCGAAGAAAGTGGTGGGCTGGCCGATGCGCTTCGGCCGCTAGAGATCCT  >  NC_002947/4902570‑4902866

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: