Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,016,880 A→G 7.4% intergenic (‑196/+31) PP_5728 ← / ← PP_5267 rhs repeat‑associated core domain‑containing protein/cytochrome c5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,016,8800AG 7.4% 266.8 / 9.1 95intergenic (‑196/+31)PP_5728/PP_5267rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base A (42/45);  new base G (3/4);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

AGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGCCATC  >  NC_002947/6016739‑6017027
                                                                                                                                              |                                                                                                                                                    
aGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAAcggcgg                                                                                                                                               >  6:21543/1‑149 (MQ=255)
   aCTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCGTTGATGGCGGGct                                                                                                                                            <  4:158106/149‑2 (MQ=255)
      ggTAAAAATGCTAGTTCTCATGGAATCAGGCTCACCGCGATTGGAGCTTCTGGAGACATCCATGGACATGCCTCTAGCGGGTTGGTTGCCTGTGCTGGACCATTCTCGGGCACCACAAAAACAAAGCCCGCCATCAGCGGCGGCCCTCa                                                                                                                                         >  4:132596/1‑149 (MQ=255)
       gTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTc                                                                                                                                          <  5:96778/147‑1 (MQ=255)
       gTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCa                                                                                                                                         <  2:171837/148‑1 (MQ=255)
        tAAAAATGCTTGTTCTCTTGGAATCAGGCTCAGCCCGATTGGTGCTTTTGGAGACATCGATGGTCATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACGACAAAAACAAAGCCCGCCGTTGATGGCGGGctttgt                                                                                                                                        <  6:101467/148‑6 (MQ=255)
          aaaaTGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCGTTGATGGCGGGc                                                                                                                                             >  4:125913/1‑141 (MQ=255)
          aaaaTGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCGTTGATGGCGGGc                                                                                                                                             <  3:125913/141‑1 (MQ=255)
          aaaaTGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGc                                                                                                                                      >  3:83076/1‑148 (MQ=255)
           aaaTGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCaa                                                                                                                                    <  8:143495/149‑1 (MQ=255)
           aaaTGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCaa                                                                                                                                    <  2:373750/149‑1 (MQ=255)
             aTGCTAGTTGTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAAtt                                                                                                                                  <  3:228131/149‑1 (MQ=255)
                cTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAACAGCCCGCCATCAACGGCGGGCCTCATGCAATTCa                                                                                                                                <  1:252031/149‑1 (MQ=255)
                cTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAg                                                                                                                               >  3:90328/1‑149 (MQ=255)
                 tAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGc                                                                                                                              >  2:379605/1‑149 (MQ=255)
                   gTTCTCATGGAATCAGGCTCAGCCCTATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCt                                                                                                                             <  7:156926/148‑1 (MQ=255)
                        cATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAAt                                                                                                                       >  5:38962/1‑149 (MQ=255)
                          tccaaTCGGGCTGAGCCTGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATca                                                                                                                     >  5:225465/4‑149 (MQ=255)
                          tccaaTCGGGCTGAGCCTGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCGGTAAAtaa                                                                                                                     >  5:111283/4‑147 (MQ=255)
                            gAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATcac                                                                                                                    >  2:263845/1‑148 (MQ=255)
                            gAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATcac                                                                                                                    <  8:102289/148‑1 (MQ=255)
                             aaTCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACaa                                                                                                                  >  8:197459/1‑149 (MQ=255)
                             aaTCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGCCGGGCCTCATGCAATTCAGCTGTAAATCACaa                                                                                                                  >  5:56464/1‑149 (MQ=255)
                                  ggCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCGTCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCg                                                                                                              >  3:323656/1‑148 (MQ=255)
                                  ggCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGCGCCTCATGCAATTCAGCTGTAAATCACAACCCgg                                                                                                             >  5:118101/1‑149 (MQ=255)
                                   gCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGa                                                                                                            <  4:83076/149‑1 (MQ=255)
                                   gCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGa                                                                                                            <  3:352134/149‑1 (MQ=255)
                                          ccagcaggAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTc                                                                                                     <  1:333556/143‑1 (MQ=255)
                                          ccGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTc                                                                                                     <  3:132596/149‑1 (MQ=255)
                                          ccGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTc                                                                                                     <  2:277842/149‑1 (MQ=255)
                                                 gAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATc                                                                                               <  2:276314/148‑1 (MQ=255)
                                                 gAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCCTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATc                                                                                               <  3:165631/148‑1 (MQ=255)
                                                    cTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAACCAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGccc                                                                                           <  2:373553/149‑1 (MQ=255)
                                                          gagaCATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCAAGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCttt                                                                                      >  7:340560/1‑147 (MQ=255)
                                                              cATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTc                                                                                <  1:379605/149‑1 (MQ=255)
                                                                  cATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTc                                                                             <  6:56464/148‑1 (MQ=255)
                                                                    tGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATc                                                                          <  5:38933/149‑1 (MQ=255)
                                                                    tGGACAGGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTATCGGGCAACACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGGAATTCAGCTGTAAATCACAACCCGGTCATTTTCTTGATCGCCCCTTTCAGGTCTTCATc                                                                          <  4:151469/149‑1 (MQ=255)
                                                                      gACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTc                                                                             >  5:261797/1‑144 (MQ=255)
                                                                      gACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTc                                                                             <  6:261797/144‑1 (MQ=255)
                                                                       aCATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTg                                                                                                  <  2:386623/123‑1 (MQ=255)
                                                                       aCATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTg                                                                                                  >  1:386623/1‑123 (MQ=255)
                                                                        cATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGa                                                                       <  5:238438/148‑1 (MQ=255)
                                                                          tGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCa                                                                    <  7:318470/149‑1 (MQ=255)
                                                                              tctAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCa                                                                            >  3:293118/1‑137 (MQ=255)
                                                                              tctAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCa                                                                            <  4:293118/137‑1 (MQ=255)
                                                                              tctAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCt                                                                >  3:80985/1‑149 (MQ=255)
                                                                                tAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTg                                                               >  8:71634/1‑148 (MQ=255)
                                                                                  gCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTg                                                               <  1:376223/146‑1 (MQ=255)
                                                                                  gCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGca                                                             >  1:179355/1‑148 (MQ=255)
                                                                                    gggTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAg                                                          <  6:118101/149‑1 (MQ=255)
                                                                                    gggTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAg                                                          >  6:362193/1‑149 (MQ=255)
                                                                                      gttggttgCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGt                                                        >  5:30931/1‑149 (MQ=255)
                                                                                           ttgCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCt                                                   >  2:135600/1‑149 (MQ=255)
                                                                                             gCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCttt                                                 >  5:221634/1‑149 (MQ=255)
                                                                                             gCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCtt                                                  <  7:197459/148‑1 (MQ=255)
                                                                                             cccTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATAAACGGCGGGCCTCATGCAATAAAGCTGTAAATCACAAACCGGACAATTTCTTGATCGCCCCTTTAAGGTCATCATCCGAGAAATCTGCACAGCTGCCCttt                                                 >  2:305106/2‑149 (MQ=255)
                                                                                               cTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTgg                                               >  7:213651/1‑149 (MQ=255)
                                                                                               cTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTgg                                               <  8:213651/149‑1 (MQ=255)
                                                                                               cTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTgg                                               <  5:362193/149‑1 (MQ=255)
                                                                                                tgtgCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTggc                                              >  1:241912/1‑149 (MQ=255)
                                                                                                tgtgCTGGACCTATCTCGGGCACAACAAAAACAAAGCCCGCCAACAACGGCGGGCCACATGCAATTCAGCTGTAAACCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTATGCACAGGTGCCCTTTggc                                              >  2:301868/1‑149 (MQ=255)
                                                                                                       gACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGGACAg                                                          <  5:289868/130‑1 (MQ=255)
                                                                                                       gACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAg                                                          >  6:289868/1‑130 (MQ=255)
                                                                                                          cTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGt                                    <  6:30931/149‑1 (MQ=255)
                                                                                                              ctcGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGa                                 <  1:263845/148‑1 (MQ=255)
                                                                                                              ctcGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGa                                 >  4:382416/1‑148 (MQ=255)
                                                                                                                 gggCACCACAAAAACAAAGCCCGCCGTCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGAt                                <  4:323656/146‑1 (MQ=255)
                                                                                                                 gggCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATa                               <  4:90328/147‑1 (MQ=255)
                                                                                                                 gggCAACAAAAAAACAAATCCCGGCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCATTCCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATAc                              <  1:305106/149‑1 (MQ=255)
                                                                                                                   gCACCACAAAAACAAAGCCCGCCATCAACGTCGGGCCTCATGCAATTCAGCTGTAAACCACAACCCGGACATTCTCTTGATCGACCCTTTCAGGTCTTCATCCGAGCAGTCTGCTCAGGTGCCCTTTGGCGGCATGGCGTTGATACCAg                           <  1:301868/149‑1 (MQ=255)
                                                                                                                       cacaAAAACAAAGCCCGCCATCCACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAAt                       >  6:357357/1‑149 (MQ=255)
                                                                                                                        acaAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGa                                 >  8:182646/1‑138 (MQ=255)
                                                                                                                        acaAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGa                                 <  7:182646/138‑1 (MQ=255)
                                                                                                                        acaAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATg                      <  3:361587/149‑1 (MQ=255)
                                                                                                                           aaaaCAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGcc                   >  8:79416/1‑149 (MQ=255)
                                                                                                                           aaaaCAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGcc                   >  2:126778/1‑149 (MQ=255)
                                                                                                                             aaCAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCACTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCtt                 >  7:234670/1‑149 (MQ=255)
                                                                                                                              catgaggCCCGCCGTTGATGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCtt                 >  4:202365/7‑149 (MQ=255)
                                                                                                                               cAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTgg               >  4:93123/1‑149 (MQ=255)
                                                                                                                                tgaggCCCGCCGTTGATGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGCGCCCTTTGGCGGCATGGCGTTGATAc                              >  1:377226/5‑134 (MQ=255)
                                                                                                                                tgaggCCCGCCGTTGATGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGCGCCCTTTGGCGGCATGGCGTTGATAc                              <  2:377226/130‑1 (MQ=255)
                                                                                                                                gaaGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGc              >  6:316417/2‑149 (MQ=255)
                                                                                                                                aaaGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCAt                                                                           <  4:90853/88‑1 (MQ=255)
                                                                                                                                aaaGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCAt                                                                           >  3:90853/1‑88 (MQ=255)
                                                                                                                                aaaGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGc              <  7:71634/149‑1 (MQ=255)
                                                                                                                                   aGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGAAACCCCCTTTAAGGTATTGATACGATCAGTCAGCACAAGTGCCCATTGGCGGCATGGCGTTGGTAACAGTAATGGCCTTGGc              >  1:314181/1‑147 (MQ=255)
                                                                                                                                      ccGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGa         >  3:343686/1‑149 (MQ=255)
                                                                                                                                       cGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGAt        <  4:343686/149‑1 (MQ=255)
                                                                                                                                       cGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGAt        <  4:80985/149‑1 (MQ=255)
                                                                                                                                         ccATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCGTTTCAGGTCTTCATCCGAGCAGGCTGCACAGGTGCCCTTTGGCGGCTTGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGc      <  8:340560/149‑1 (MQ=255)
                                                                                                                                          cATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGc      >  3:174849/1‑148 (MQ=255)
                                                                                                                                           aTCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGCCa    >  3:222877/1‑149 (MQ=255)
                                                                                                                                             cAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGGTACCAGTAATGGCCt                  >  1:233690/1‑133 (MQ=255)
                                                                                                                                             cAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGGTACCAGTAATGGCCt                  <  2:233690/133‑1 (MQ=255)
                                                                                                                                              aaCGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGCCATc  <  5:21543/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGCCATCAACGGCGGGCCTCATGCAATTCAGCTGTAAATCACAACCCGGACATTTTCTTGATCGCCCCTTTCAGGTCTTCATCCGAGCAGTCTGCACAGGTGCCCTTTGGCGGCATGGCGTTGATACCAGTAATGGCCTTGGCCAGGATGCCATC  >  NC_002947/6016739‑6017027

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: