Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,982,882 T→A 10.5% H119Q (CAT→CAA PP_3510 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,982,8820TA10.5% 256.8 / 15.4 89H119Q (CAT→CAAPP_3510hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/39);  new base A (4/5);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAG  >  NC_002947/3982733‑3983024
                                                                                                                                                     |                                                                                                                                              
gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                >  2:988570/1‑150 (MQ=255)
gCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAt                                                                                                                                                >  1:1632074/1‑150 (MQ=255)
 cAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAtt                                                                                                                                               <  1:1608110/150‑1 (MQ=255)
 cAAAGGCACGTATGCTGGAACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCAtt                                                                                                                                               <  2:591408/150‑1 (MQ=255)
  aaaGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTg                                                                                                                                              >  1:58012/1‑150 (MQ=255)
   aaGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGa                                                                                                                                             <  2:2293231/150‑1 (MQ=255)
    aGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGAc                                                                                                                                            <  1:1638685/150‑1 (MQ=255)
     ggCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCCCGCCGGCATGCTCCTCCGCCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGAcc                                                                                                                                           <  2:1269504/150‑1 (MQ=255)
          gTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTCCACCGCTGGGACTGGGGAAAACGCCCGCAGGGTCAATGCCCCTGca                                                                                                                                      <  2:467934/150‑2 (MQ=255)
            aTGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCggg                                                                                                                                    <  1:988570/150‑2 (MQ=255)
            aTGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                    <  2:644707/150‑1 (MQ=255)
                 ggCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGATCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGAc                                                                                                                               <  2:2245031/150‑1 (MQ=255)
                          cTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATc                                                                                                                      >  2:502609/1‑150 (MQ=255)
                               ttgCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTc                                                                                                                 <  1:610782/148‑1 (MQ=255)
                               aGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACTTTAGTATCGCTTc                                                                                                                 >  1:497868/1‑150 (MQ=255)
                                      ggACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCt                                                                                                          >  1:699799/1‑150 (MQ=255)
                                        aCCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTAt                                                                                                        <  1:2011884/150‑1 (MQ=255)
                                         ccACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATc                                                                                                       >  2:1535954/1‑150 (MQ=255)
                                          cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCt                                                                                                                                         >  1:789032/1‑115 (MQ=255)
                                          cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCt                                                                                                                                         <  2:789032/115‑1 (MQ=255)
                                          cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAGTGCCCCTGcaggcgt                                                                                                                                 >  2:1940828/1‑117 (MQ=255)
                                          cACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCCGCAGGGTCAGTGCCCCTGcaggcgt                                                                                                                                 <  1:1940828/123‑7 (MQ=255)
                                             gCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAtt                                                                                                   >  2:1485022/1‑150 (MQ=255)
                                             gCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAAtt                                                                                                   >  1:344225/1‑150 (MQ=255)
                                              ccGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTc                                                                                                  <  1:2318716/150‑1 (MQ=255)
                                               cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATa                                                                                                                           >  2:2344360/1‑124 (MQ=255)
                                               cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATa                                                                                                                           <  1:2344360/124‑1 (MQ=255)
                                               cGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCg                                                                                                 >  2:1466185/1‑150 (MQ=255)
                                                  aaTACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                              <  1:1485022/150‑1 (MQ=255)
                                                    tACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                   <  1:737067/111‑5 (MQ=255)
                                                    tACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGcaggc                                                                                                                                   >  2:737067/1‑107 (MQ=255)
                                                       aGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGc                                                                                         >  2:1275188/1‑150 (MQ=255)
                                                        gCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGcc                                                                                        <  2:1632074/150‑1 (MQ=255)
                                                         cTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCg                                                                                       >  1:1440795/1‑150 (MQ=255)
                                                               cATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCg                                                                                 >  2:228013/1‑150 (MQ=255)
                                                                  cccGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                    >  1:2305838/1‑96 (MQ=255)
                                                                  cccGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                    <  2:2305838/96‑1 (MQ=255)
                                                                  cccGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                              >  1:241268/1‑134 (MQ=255)
                                                                  cccGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                              <  2:241268/134‑1 (MQ=255)
                                                                     gCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                  <  2:1103853/95‑1 (MQ=255)
                                                                     gCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGttt                                                                                                                                  >  1:1103853/1‑95 (MQ=255)
                                                                           cATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGGGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAgca                                                                     >  2:57026/1‑148 (MQ=255)
                                                                                  gCTCTGACAGAGATGACGCGTGTTCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTa                                                              >  2:1306408/1‑150 (MQ=255)
                                                                                  gCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGca                                                                                                                                      <  1:2315161/78‑2 (MQ=255)
                                                                                  gCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGTCAATGCCCCTGca                                                                                                                                      >  2:2315161/1‑77 (MQ=255)
                                                                                  gCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTa                                                              >  2:2093333/1‑150 (MQ=255)
                                                                                    tATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAg                                                            <  2:2288253/150‑1 (MQ=255)
                                                                                        aCAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTgg                                                        <  1:502609/150‑1 (MQ=255)
                                                                                         cAGCGAAGAGTCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                              <  1:2106249/111‑1 (MQ=255)
                                                                                         cAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGa                                                                                              >  2:2106249/1‑111 (MQ=255)
                                                                                           gCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                    >  2:1257643/1‑71 (MQ=255)
                                                                                           gCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGt                                                                                                                                    <  1:1257643/71‑1 (MQ=255)
                                                                                               agagGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTg                                                 <  1:1466185/150‑1 (MQ=255)
                                                                                                   gCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCa                                             <  2:1440795/150‑1 (MQ=255)
                                                                                                     cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTc                                                                                                                 >  1:2058015/1‑80 (MQ=255)
                                                                                                     cTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTc                                                                                                                 <  2:2058015/80‑1 (MQ=255)
                                                                                                        tgCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATaa                                        <  2:58012/150‑1 (MQ=255)
                                                                                                            cGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCCCACCCTATCAATTCGTGACTGGGCTTTCCCGAGCCGCCAAGGATGTATTGAGGTGGCCTGGTGAGCATATAAcagc                                    <  2:645210/150‑1 (MQ=255)
                                                                                                             gggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGt                                                                  <  1:435148/119‑1 (MQ=255)
                                                                                                             gggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGt                                                                  >  2:435148/1‑119 (MQ=255)
                                                                                                             gggTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagca                                   <  1:228013/150‑1 (MQ=255)
                                                                                                                   aTGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGa                             >  2:456260/1‑150 (MQ=255)
                                                                                                                     ggCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACt                                                                                            >  1:1461257/1‑85 (MQ=255)
                                                                                                                     ggCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACt                                                                                            <  2:1461257/85‑1 (MQ=255)
                                                                                                                     ggCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGc                           <  1:1275188/150‑1 (MQ=255)
                                                                                                                        ggggACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTTGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                      <  1:974983/136‑1 (MQ=255)
                                                                                                                        ggggACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAca                                      >  2:974983/1‑136 (MQ=255)
                                                                                                                            aCTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCggg                    <  1:932014/150‑1 (MQ=255)
                                                                                                                              tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAg                            >  1:1402953/1‑140 (MQ=255)
                                                                                                                              tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAg                            <  2:1402953/140‑1 (MQ=255)
                                                                                                                              tGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCa                  >  1:19410/1‑150 (MQ=255)
                                                                                                                                 ggAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAATATCGCTTCGCACCCTATCAATTCGTGACTGACCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACtt               >  2:210141/1‑150 (MQ=255)
                                                                                                                                    aaaCGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGATCGCTCGGGctgctgtt            <  1:1788420/150‑8 (MQ=255)
                                                                                                                                        gCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGAc        >  2:792161/1‑150 (MQ=255)
                                                                                                                                        aaccgCAGGGTCAATGCCCCTGCAGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTg                                                         >  2:655233/5‑101 (MQ=255)
                                                                                                                                        aaccgCAGGGTCAATGCCCCCGCAGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTg                                                         <  1:655233/97‑1 (MQ=255)
                                                                                                                                         ccTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACa       >  2:69024/1‑150 (MQ=255)
                                                                                                                                          ccgCAGGGTCAATGCCCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGa                             <  2:1304066/125‑1 (MQ=255)
                                                                                                                                          ccgCAGGGTCAATGCCCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGa                             >  1:1304066/3‑127 (MQ=255)
                                                                                                                                          cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTAt                                                                                                        >  2:1365372/1‑52 (MQ=255)
                                                                                                                                          cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTAt                                                                                                        <  1:1365372/52‑1 (MQ=255)
                                                                                                                                          cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACaa      >  1:1555748/1‑150 (MQ=255)
                                                                                                                                          cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACaa      >  1:1555688/1‑150 (MQ=255)
                                                                                                                                          cTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACaa      >  2:1812831/1‑150 (MQ=255)
                                                                                                                                            gCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagc                                    >  1:512118/1‑118 (MQ=255)
                                                                                                                                            gCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcagc                                    <  2:512118/118‑1 (MQ=255)
                                                                                                                                              aGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAg  <  1:133007/150‑1 (MQ=255)
                                                                                                                                                gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcag                                     <  2:2063687/113‑1 (MQ=255)
                                                                                                                                                gggCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAAcag                                     >  1:2063687/1‑113 (MQ=255)
                                                                                                                                                  gCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa           <  2:1894251/137‑1 (MQ=255)
                                                                                                                                                  gCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCaa           >  1:1894251/1‑137 (MQ=255)
                                                                                                                                                   cATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAtt                                                               <  1:404831/84‑1 (MQ=255)
                                                                                                                                                   cATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTAtt                                                               >  2:404831/1‑84 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCAAAGGCACGTATGCTGGCACTGACCTGCCAGGCATTGGACCACGCCGGAATACAGCTGGACCATCCCGCTTGGCATGTACGCTATGACAGCGAAGAGGCCTGTGCACGGGTGCATGGCGGGGACTGGGGAAAACGCCTGCAGGGGCATTGACCCTGCGGTTTGACGATAGTATCGCTTCGCACCCTATCAATTCGTGACTGGCCGTTCCCGATCCGCCAAGGATGTATTAAGGTGGCATGGTGCGCATATAACAGCAGCCCGAGCGATCGGGCACTTGCAAGACAATGAG  >  NC_002947/3982733‑3983024

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: