Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,129,376 C→G 6.8% Q692H (CAG→CAC) ‡ recQ ← ATP‑dependent DNA 3'‑5' helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,129,3760CG 6.8% 309.0 / 6.5 97Q692H (CAG→CAC) ‡recQATP‑dependent DNA 3'‑5' helicase
Reads supporting (aligned to +/- strand):  ref base C (44/46);  new base G (4/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGC  >  NC_002947/5129232‑5129518
                                                                                                                                                 |                                                                                                                                              
cgcgTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgg                                                                                                                                           <  2:1380178/150‑1 (MQ=255)
cgcgGTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGTCGTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAACTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgg                                                                                                                                           <  1:264229/150‑1 (MQ=255)
   gTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGACCTCCGGCACCCGCGCCCCGAACTGCGGGGc                                                                                                                                        >  1:1929425/1‑150 (MQ=255)
    ttCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGTAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCa                                                                                                                                       >  2:1740784/1‑150 (MQ=255)
       gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCCCACTGCGAGGCAATg                                                                                                                                    >  2:1376837/1‑150 (MQ=255)
       gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTACGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgg                                                                                                                                    <  2:1564945/150‑5 (MQ=255)
         tAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGc                                                                                                                                  <  1:436651/150‑1 (MQ=255)
         tAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGc                                                                                                                                  <  2:2047743/150‑1 (MQ=255)
                 aaaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGcc                                                                                                                                 <  1:1584198/143‑1 (MQ=255)
                 aaaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGcc                                                                                                                                 >  2:1584198/1‑143 (MQ=255)
                  aaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATGCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCgg                                                                                                                         <  2:1081637/150‑1 (MQ=255)
                       cATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCCGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACGGGTGGCa                                                                                                                    <  2:532629/150‑1 (MQ=255)
                         tGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGc                                                                                                                  >  2:846884/1‑150 (MQ=255)
                         tGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTAGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGc                                                                                                                  <  1:1602786/150‑1 (MQ=255)
                          gACAAAGGTTAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCt                                                                                                                 >  2:2284346/1‑150 (MQ=255)
                          gACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGAGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCt                                                                                                                 >  1:1523290/1‑150 (MQ=255)
                           aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgggt                                                                                                                                  <  1:1613440/132‑7 (MQ=255)
                           aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgggt                                                                                                                                  >  2:1613440/1‑126 (MQ=255)
                             aaaGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAc                                                                                                              >  1:956304/1‑150 (MQ=255)
                             aaaGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAc                                                                                                              >  1:129138/1‑150 (MQ=255)
                             aaaGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAc                                                                                                              >  2:1277008/1‑150 (MQ=255)
                             aaaGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAc                                                                                                              >  1:2127460/1‑150 (MQ=255)
                              aaGTTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAcc                                                                                                             <  2:418630/150‑1 (MQ=255)
                               aGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgg                                                                                                                                           >  1:1217353/1‑119 (MQ=255)
                               aGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCgg                                                                                                                                           <  2:1217353/119‑1 (MQ=255)
                                ggTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCt                                                                                                           >  2:134238/1‑150 (MQ=255)
                                  tCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCg                                                                                                         >  2:2350341/1‑150 (MQ=255)
                                   cagcagATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGc                                                                                                        <  2:2348342/150‑1 (MQ=255)
                                      cagATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCAt                                                                                                     >  2:2182747/1‑150 (MQ=255)
                                       agATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATc                                                                                                    <  2:956304/150‑1 (MQ=255)
                                                ttCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGc                                                                                           <  1:1715200/150‑1 (MQ=255)
                                                 tCTGCCTTCGATCCTGCATGCGTTTACAATTCGACTTCCGCCGCAAGGTGCTCGGGCACGCAATAGAGTACGCACTACGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCg                                                                                          <  1:1740784/150‑1 (MQ=255)
                                                    gCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              >  1:732412/1‑143 (MQ=255)
                                                    gCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              <  2:732412/143‑1 (MQ=255)
                                                     ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCt                                                                                      <  1:1176282/150‑1 (MQ=255)
                                                     ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCt                                                                                      <  1:1277008/150‑1 (MQ=255)
                                                      cTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              <  2:501851/141‑1 (MQ=255)
                                                      cTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                              >  1:501851/1‑141 (MQ=255)
                                                      cTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTg                                                                                     >  1:504283/1‑150 (MQ=255)
                                                      cTTCGATCCTGCATGCGTTCACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTg                                                                                     <  2:1584840/150‑1 (MQ=255)
                                                       ttCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTTCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGc                                                                                    >  1:183378/1‑150 (MQ=255)
                                                        tCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCa                                                                                   <  2:2127460/150‑1 (MQ=255)
                                                          gATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                          <  1:2278874/93‑1 (MQ=255)
                                                          gATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCggg                                                                                                                                          >  2:2278874/1‑93 (MQ=255)
                                                            tCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTc                                                                               <  1:2284346/150‑1 (MQ=255)
                                                              cTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCca                                                                             <  2:431602/150‑1 (MQ=255)
                                                               tGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACTCCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCcat                                                                            >  2:1921738/1‑150 (MQ=255)
                                                                   tGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCatcatc                                                                        <  1:846884/150‑1 (MQ=255)
                                                                    gCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATc                                                                                                    >  1:89690/1‑121 (MQ=255)
                                                                    gCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATc                                                                                                    <  2:89690/121‑1 (MQ=255)
                                                                     cGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAg                                                                      <  2:1051142/150‑1 (MQ=255)
                                                                      gTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAgg                                                                     >  2:1829492/1‑150 (MQ=255)
                                                                       tttACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGt                                                                    >  2:2239027/1‑150 (MQ=255)
                                                                       tttACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGt                                                                    <  1:425919/150‑1 (MQ=255)
                                                                       tttACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGt                                                                    >  1:289745/1‑150 (MQ=255)
                                                                             aTTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAcc                                                                                                             <  2:2043012/103‑1 (MQ=255)
                                                                             aTTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAcc                                                                                                             >  1:2043012/1‑103 (MQ=255)
                                                                              ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgc                                                                                                                                      >  2:2322630/1‑75 (MQ=255)
                                                                              ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCgg                                                             >  1:290020/1‑150 (MQ=255)
                                                                              ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCGCCGGCACCCGCGCCCCGCAGTTCGGGGcgc                                                                                                                                      <  1:2322630/77‑3 (MQ=255)
                                                                               tCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGc                                                            >  1:468220/1‑150 (MQ=255)
                                                                                 gAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAg                                                          >  1:688213/1‑150 (MQ=255)
                                                                                  aaTTCCGCCGCCCGTGCTGCGCGCACTCAAGAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCCTTGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAgg                                                         <  1:1731247/150‑1 (MQ=255)
                                                                                  aaTTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGCAAGGCGTCCTGCACCTCCATCAGCAGGTCTTccggccgg                                                         >  2:523151/1‑150 (MQ=255)
                                                                                    ttCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTc                                                       >  1:1795637/1‑150 (MQ=255)
                                                                                      ccgccgCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGa                                                     <  1:134238/150‑1 (MQ=255)
                                                                                       cgccgcCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGCCAGGTCGAt                                                    >  2:1979453/1‑150 (MQ=255)
                                                                                          cgccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGc                                                 <  2:290020/150‑1 (MQ=255)
                                                                                                gTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAg                                                                                             >  1:393359/1‑101 (MQ=255)
                                                                                                gTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAg                                                                                             <  2:393359/101‑1 (MQ=255)
                                                                                                  gctgcGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTcc                                          <  2:1523290/150‑1 (MQ=255)
                                                                                                   ctgcGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCa                                         <  2:289745/150‑1 (MQ=255)
                                                                                                     gcgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                       <  1:41229/150‑1 (MQ=255)
                                                                                                     gcgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                       <  2:129138/150‑1 (MQ=255)
                                                                                                     gcgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                       <  2:1499817/150‑1 (MQ=255)
                                                                                                               aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGcc                                                >  2:256029/1‑131 (MQ=255)
                                                                                                               aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGcc                                                <  1:256029/131‑1 (MQ=255)
                                                                                                               aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                       >  1:110086/1‑140 (MQ=255)
                                                                                                               aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                       <  2:110086/140‑1 (MQ=255)
                                                                                                                 tAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTgg                           >  1:2032681/1‑150 (MQ=255)
                                                                                                                 tAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTgg                           >  2:1142992/1‑150 (MQ=255)
                                                                                                                 tAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCCCCTGCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTgg                           <  1:1376837/150‑1 (MQ=255)
                                                                                                                  agagTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggc                          <  2:504283/150‑1 (MQ=255)
                                                                                                                   gagTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggcg                         >  1:1219963/1‑150 (MQ=255)
                                                                                                                         gCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCaa                   >  1:723224/1‑150 (MQ=255)
                                                                                                                           ccTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCg                                                                                                         >  1:269471/1‑62 (MQ=255)
                                                                                                                           ccTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCg                                                                                                         <  2:269471/62‑1 (MQ=255)
                                                                                                                               cGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTc             >  1:203421/1‑150 (MQ=255)
                                                                                                                               cGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTc             >  2:736785/1‑150 (MQ=255)
                                                                                                                                ggCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg            <  2:688213/150‑1 (MQ=255)
                                                                                                                                  ggccattgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGc          >  1:1079136/8‑150 (MQ=255)
                                                                                                                                     ccGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGc                                 <  1:705301/124‑1 (MQ=255)
                                                                                                                                     ccGAGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGc                                 >  2:705301/1‑124 (MQ=255)
                                                                                                                                      cttgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCAGGCAGGTCGATGGCTTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCcagc      >  1:266117/4‑150 (MQ=255)
                                                                                                                                      cgcgCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCcagc      >  1:306233/1‑150 (MQ=255)
                                                                                                                                        cgcCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCcagcag    <  2:306233/150‑1 (MQ=255)
                                                                                                                                          ccccGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGc  <  2:203421/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGC  >  NC_002947/5129232‑5129518

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: