Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,129,326 A→C 6.3% F693C (TTC→TGC)  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,3260AC 6.3% 296.2 / 7.8 96F693C (TTC→TGC) recQATP‑dependent DNA 3'‑5' helicase
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base C (2/4);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GTTACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGCTAT  >  NC_002947/5129180‑5129473
                                                                                                                                                    |                                                                                                                                                   
gTTACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGaa                                                                                                                                                    <  7:374558/148‑1 (MQ=255)
 ttACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACt                                                                                                                                                  >  3:26883/1‑149 (MQ=255)
   aCGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTAAGAAGATGGCGCTTCTGCCTTCGATCATGCATGCGTTTACAATTCGAATTCCGCAGCCATTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGc                                                                                                                                                >  7:171436/1‑149 (MQ=255)
    cgcgTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTCCGCCCTCCGGCACCCGCGCCCCGAACTGCg                                                                                                                                               <  4:67235/149‑1 (MQ=255)
           gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgc                                                                                                                                         <  4:9097/148‑3 (MQ=255)
           gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgc                                                                                                                                         <  2:277297/148‑3 (MQ=255)
           gATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAAt                                                                                                                                        >  7:362817/1‑149 (MQ=255)
            aTAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATg                                                                                                                                       <  8:381233/149‑1 (MQ=255)
            aTAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATg                                                                                                                                       <  3:42396/149‑1 (MQ=255)
            aTAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGAACCCGCGACCCGAACTGCGGGGCAATg                                                                                                                                       >  3:127054/1‑149 (MQ=255)
               aTCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgggtg                                                                                                                                    <  5:298447/149‑8 (MQ=255)
               aTCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGcc                                                                                                                                    <  7:83750/149‑1 (MQ=255)
                tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCTGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGcc                                                                                                                                    >  1:278547/1‑149 (MQ=255)
                tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgg                                                                                                                                       >  2:179715/1‑141 (MQ=255)
                tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGCAGTTCGGGGcgcgg                                                                                                                                       <  1:179715/145‑5 (MQ=255)
                tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGcc                                                                                                                                    >  3:211962/1‑148 (MQ=255)
                tCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCg                                                                                                                                   >  2:307436/1‑149 (MQ=255)
                  gTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTCCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGa                                                                                                                                  >  5:201100/1‑148 (MQ=255)
                   tCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGaca                                                                                                                                <  6:312456/149‑1 (MQ=255)
                    cAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTGCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGaca                                                                                                                                >  1:278739/1‑149 (MQ=255)
                    aaaaaaGAATGAGAAAGGTAAGCAGATGGCGCTTCTGCCTTCGATACTGCATGAGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGAGCACGAAATAGAGTACGACCTCCGGAACCCGCGACCCGAAATGCGGGGCAATGGCCGaca                                                                                                                                <  4:127054/147‑1 (MQ=255)
                      aaaaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGGTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCg                                                                                                                             <  8:362104/149‑1 (MQ=255)
                        aaGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGt                                                                                                                           <  2:278547/149‑1 (MQ=255)
                          gCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTgg                                                                                                                         <  8:171436/149‑1 (MQ=255)
                            aTGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCa                                                                                                                       <  6:20618/149‑1 (MQ=255)
                              gACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGc                                                                                                                     >  8:300790/1‑149 (MQ=255)
                               aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGc                                                                                                                     <  7:296409/148‑1 (MQ=255)
                               aCAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCt                                                                                                                    >  7:123807/1‑149 (MQ=255)
                                 aaaGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCa                                                                                                                  >  3:154348/1‑149 (MQ=255)
                                   aGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCAcc                                                                                                                >  5:122380/1‑149 (MQ=255)
                                     gTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCt                                                                                                              >  6:52913/1‑149 (MQ=255)
                                       cagcagATCGCGCTTCTGCCTTCGATCCTGAATGCCTTTACAAATCGAATTCCGCCGCCAGTTCTGCGAGAACGAAATAGAGTACGCCCTCCGGCACCCGCTCCCCGAACTGCGGGGCAATGGCCGACACCAGTGGCAGCTCACCCtag                                                                                                            >  4:258485/1‑147 (MQ=255)
                                          cagATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCa                                                                                                         <  5:23802/149‑1 (MQ=255)
                                               ggCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATc                                                                                                       >  5:277293/1‑146 (MQ=255)
                                                 cgcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaa                                                                                                  <  2:278739/149‑1 (MQ=255)
                                                   gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                                 <  1:68674/149‑1 (MQ=255)
                                                   gcTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                                 <  3:118089/149‑1 (MQ=255)
                                                    cTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAg                                                                                                <  8:49023/149‑1 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                                 <  6:156295/142‑1 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGaaa                                                                                                 >  5:156295/1‑142 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGCAAGGCGTcc                                                                                          >  4:247362/1‑149 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTc                                                                                           <  1:420621/148‑1 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTc                                                                                           <  5:52913/148‑1 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAATACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTc                                                                                                                   >  7:125172/1‑124 (MQ=255)
                                                           ccTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAATACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTc                                                                                                                   <  8:125172/124‑1 (MQ=255)
                                                               cGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTg                                                                                        <  6:357189/147‑1 (MQ=255)
                                                                  tCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCt                                                                                   >  8:202808/1‑149 (MQ=255)
                                                                    cTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTcc                                                                                 <  4:90991/149‑1 (MQ=255)
                                                                    cTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTcc                                                                                 <  7:8756/149‑1 (MQ=255)
                                                                      gCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCcat                                                                               >  3:167130/1‑149 (MQ=255)
                                                                       cATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCcatc                                                                              <  4:211962/149‑1 (MQ=255)
                                                                       cATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCcat                                                                               >  7:198631/1‑148 (MQ=255)
                                                                         tGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCcatca                                                                             >  8:322867/1‑148 (MQ=255)
                                                                           cGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCAtcatca                                                                          <  2:218613/149‑1 (MQ=255)
                                                                           cGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCAtcatca                                                                          <  1:307436/149‑1 (MQ=255)
                                                                            gTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAg                                                                         <  5:427669/149‑1 (MQ=255)
                                                                             tttACAATTAGAATTCCGCCGGCTGTTCTGCGCGCACGCATTTGAGTACGCCCTCCGGCACCAGCGCGCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAgg                                                                        <  3:258485/149‑1 (MQ=255)
                                                                               tACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAg                                                                         >  4:335233/1‑146 (MQ=255)
                                                                               tACAATTCGAATTCCCCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAg                                                                         <  3:335233/146‑1 (MQ=255)
                                                                                aCAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCt                                                                     >  2:163649/1‑149 (MQ=255)
                                                                                aCAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCt                                                                     >  2:80500/1‑149 (MQ=255)
                                                                                  aaTTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTc                                                                   >  7:9326/1‑149 (MQ=255)
                                                                                   aTTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTcc                                                                  <  1:174785/149‑1 (MQ=255)
                                                                                   aTTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCACCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTcc                                                                  >  2:280189/1‑149 (MQ=255)
                                                                                    ttCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCg                                                                 >  2:369937/1‑149 (MQ=255)
                                                                                      cGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCgg                                                                <  8:9326/148‑1 (MQ=255)
                                                                                       gAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCgg                                                                >  6:3105/1‑147 (MQ=255)
                                                                                              gccgccAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGAt                                                       >  4:258969/1‑149 (MQ=255)
                                                                                                    aGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTg                                                 <  8:362817/149‑1 (MQ=255)
                                                                                                      tgctgcGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGt                                               <  1:44505/149‑1 (MQ=255)
                                                                                                         tgcGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTc                                              >  7:52584/1‑147 (MQ=255)
                                                                                                          gcgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTc                                              >  7:197373/1‑146 (MQ=255)
                                                                                                           cgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCa                                            <  8:123678/147‑1 (MQ=255)
                                                                                                           cgcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCa                                            <  7:202808/147‑1 (MQ=255)
                                                                                                            gcgcACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGa                                          >  4:69802/1‑148 (MQ=255)
                                                                                                             cgcacgcaATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACt                                        >  8:258692/1‑149 (MQ=255)
                                                                                                                acgcaATAGAGTACGCCCTCTGGCACCCGCGCCCCGAACTGCCGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAg                                     >  4:409589/1‑149 (MQ=255)
                                                                                                                acgcaATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAg                                     <  1:163649/149‑1 (MQ=255)
                                                                                                                  gcaATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCt                                   >  4:194318/1‑149 (MQ=255)
                                                                                                                  gcaATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCt                                   >  6:339879/1‑149 (MQ=255)
                                                                                                                    aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTcc                                 <  8:197373/149‑1 (MQ=255)
                                                                                                                    aaTAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGCCTCCTGCACCTACATGATGAGGTATTCCCGCAGGTCGATGGCCTGGTCCAGACTCCGCTCCt                                >  3:252280/1‑149 (MQ=255)
                                                                                                                          gTACGCCCTCCGGCACCCGCACCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggcgg                           >  2:5476/1‑149 (MQ=255)
                                                                                                                           tACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggcggc                          >  7:403741/1‑149 (MQ=255)
                                                                                                                           tACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggcggc                          >  7:118861/1‑149 (MQ=255)
                                                                                                                           tACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTggcgg                           <  6:277293/148‑1 (MQ=255)
                                                                                                                                 cTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACg                    >  1:428877/1‑149 (MQ=255)
                                                                                                                                 cTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACg                    >  1:392049/1‑149 (MQ=255)
                                                                                                                                      gcccccGCTCCCCCTTCTGCGGGGCAATGGCCGTCCCCGGGGGCATCTCTCCCTTGCCTTCTTTTTCTGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg               <  2:392049/146‑1 (MQ=255)
                                                                                                                                      gCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCg               >  3:180861/1‑149 (MQ=255)
                                                                                                                                          ccGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCGCGGc             >  6:425056/1‑147 (MQ=255)
                                                                                                                                           attgcCCCGCAGTTCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCGTGCTCCGCCATCCTCAGGTCTTCCGGCAGGTCGATGGCCTGGTCGAGGCTCAGCTCCTGGCGGCGCGACGCCGGGGCcag          >  8:377486/4‑149 (MQ=255)
                                                                                                                                            gcgcCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCcagc         <  7:300790/149‑1 (MQ=255)
                                                                                                                                                cccGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGc     >  7:356008/1‑149 (MQ=255)
                                                                                                                                                  cGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGCta   <  3:409589/149‑1 (MQ=255)
                                                                                                                                                  cGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGCta   >  7:293942/1‑149 (MQ=255)
                                                                                                                                                   gAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGCtat  <  8:52584/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTACGCGTTCGATAATCGTCAAAAAGCATGACAAAGGTCAGCAGATGGCGCTTCTGCCTTCGATCCTGCATGCGTTTACAATTCGAATTCCGCCGCCAGTGCTGCGCGCACGCAATAGAGTACGCCCTCCGGCACCCGCGCCCCGAACTGCGGGGCAATGGCCGACACCGGTGGCAGCTCACCCTCGCCATCAAGAAAGGCGTCCTGCACCTCCATCATCAGGTCTTCCGGCAGGTCGATGGCCTGGTCCAGACTCAGCTCCTGGCGGCCCAACGCCTCGGCCAGCAGGCTAT  >  NC_002947/5129180‑5129473

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