Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 96,131 Δ1 bp 8.5% coding (554/1005 nt) PP_0088 ← alkanal monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294796,1310C. 8.5% 327.9 / 25.8 106coding (554/1005 nt)PP_0088alkanal monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (45/52);  new base . (4/5);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

Only 100 of 108 total aligned reads displayed.

TCTCGTCGGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGAA  >  NC_002947/95982‑96279
                                                                                                                                                     |                                                                                                                                                     
tctcGTCGGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGc                                                                                                                                                       >  2:2038592/1‑150 (MQ=255)
    atcgGTCTGCGCCACCACCATGGGAATACCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACACGAATGGCCTCGTGAATATAGCGTGGCGCGAAATGCAAGGCGAAGGCGTAAGGCATTcccgg                                                                                                                                                   >  1:1893408/2‑147 (MQ=255)
      cgGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCa                                                                                                                                                <  2:1261731/150‑1 (MQ=255)
       ggTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCAt                                                                                                                                               <  1:1115840/150‑1 (MQ=255)
        gTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATg                                                                                                                                              >  1:937169/1‑150 (MQ=255)
         tCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGc                                                                                                                                             <  2:1111862/150‑1 (MQ=255)
             cgcCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGc                                                                                                                                         <  1:1226647/150‑1 (MQ=255)
              gcCACCACCATGGGAATACCCCGCATCACATACGACCTGTCCAGCTGTGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctac                                                                                                                                        <  2:2071607/150‑4 (MQ=255)
               ccaccaccaTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCa                                                                                                                                       <  1:282317/150‑1 (MQ=255)
                   caccaTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCCAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTg                                                                                                                                   <  1:2292631/150‑1 (MQ=255)
                      caTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctac                                                                                                                                        >  2:2286463/1‑139 (MQ=255)
                      caTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctac                                                                                                                                        <  1:2286463/142‑4 (MQ=255)
                           gAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGa                                                                                                                           <  2:937169/150‑1 (MQ=255)
                              tGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAGTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACa                                                                                                                        <  1:1343716/150‑1 (MQ=255)
                                cccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAgg                                                                                                                      >  2:1902796/1‑150 (MQ=255)
                                 ccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGc                                                                                                                     <  2:615564/150‑1 (MQ=255)
                                 ccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGc                                                                                                                     >  1:1663821/1‑150 (MQ=255)
                                   aGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCCCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTg                                                                                                                   <  1:1280039/150‑1 (MQ=255)
                                    gaaTCACATACGGCTTGTCAAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAACTGGGCACTGAAc                                                                                                                         >  1:940811/3‑143 (MQ=255)
                                    gaaTCACATACGGCTTGTCAAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAACTGGGCACTGAAc                                                                                                                         <  2:940811/141‑1 (MQ=255)
                                      aTCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGAt                                                                                                                >  2:1535125/1‑150 (MQ=255)
                                         acaTACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATccc                                                                                                             <  2:1663821/150‑1 (MQ=255)
                                          caTACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCCGCGGTGCACTGAACAGGCTGGATCCCa                                                                                                            <  1:1530960/150‑1 (MQ=255)
                                          caTACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCa                                                                                                            <  1:142380/150‑1 (MQ=255)
                                          caTACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCa                                                                                                            <  1:1902796/150‑1 (MQ=255)
                                                gCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTAGACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCgcgc                                                                                                                                                 >  1:424874/1‑107 (MQ=255)
                                                gCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCgcgc                                                                                                                                                 <  2:424874/107‑1 (MQ=255)
                                                  ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTTGGGCATGCGCGGCATGccctac                                                                                                                                        >  1:1861783/1‑111 (MQ=255)
                                                  ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTTGGGCATGCGCGGCATGccctac                                                                                                                                        <  2:1861783/114‑4 (MQ=255)
                                                  ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAAc                                                                                                       >  1:909024/1‑147 (MQ=255)
                                                  ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAAc                                                                                                       <  2:909024/147‑1 (MQ=255)
                                                     tCCAGCGTGGTCGGGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTGGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATc                                                                                                 <  1:921328/150‑1 (MQ=255)
                                                      ccAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGCGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCg                                                                                                <  2:285415/150‑1 (MQ=255)
                                                       cAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCgg                                                                                               >  2:1442826/1‑150 (MQ=255)
                                                          cGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCAc                                                                                            <  2:2267380/150‑1 (MQ=255)
                                                               tCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACcac                                                                                                    >  2:2288606/1‑137 (MQ=255)
                                                               tCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACcac                                                                                                    <  1:2288606/137‑1 (MQ=255)
                                                                cGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACt                                                                                                                             >  2:795571/1‑111 (MQ=255)
                                                                cGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACt                                                                                                                             <  1:795571/111‑1 (MQ=255)
                                                                cGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGgt                                                                                      <  1:1476610/150‑1 (MQ=255)
                                                                  aGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGgtgt                                                                                    >  2:1800733/1‑150 (MQ=255)
                                                                     gCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTggg                                                                                                                                 >  2:806504/1‑102 (MQ=255)
                                                                     gCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTggg                                                                                                                                 <  1:806504/102‑1 (MQ=255)
                                                                        tGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGACGTAGGGCATGCACGTCTGGCCCGCCAGCTGGGCACTGAACAGGCTGGATTCCAGAAACCCAATGGGATCGTCTGTGTCGTAcc                                                                              >  1:761379/1‑150 (MQ=255)
                                                                         gAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTAccc                                                                             <  1:1800733/150‑1 (MQ=255)
                                                                          aaaTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAATGGTCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTAcccc                                                                            >  1:1467206/1‑150 (MQ=255)
                                                                          aaaTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTAcccc                                                                            <  2:2097825/150‑1 (MQ=255)
                                                                           aaTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCg                                                                           <  1:1906619/150‑1 (MQ=255)
                                                                            aTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctaa                                                                                                                                        >  2:2132141/1‑85 (MQ=255)
                                                                            aTGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctac                                                                                                                                        <  1:2132141/88‑4 (MQ=255)
                                                                               gTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCGCATCGGCACGt                                                                                          >  1:1235596/1‑131 (MQ=255)
                                                                               gTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCACCAACCGCATCGGCACGt                                                                                          <  2:1235596/131‑1 (MQ=255)
                                                                                  gCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGTCCGCCAGCTGGGCACTGAGCAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGc                                                                    <  1:433387/150‑1 (MQ=255)
                                                                                    ggTACACGCGAATGGCTTCGTGCATGTAGCGTGGAGTGAATTGCGAGGCGAAAGCGTACGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATTCCACCAACCCAATCCATACGTCGGTGTTTTATACCGGTACGGCGa                                                                  <  2:760032/150‑1 (MQ=255)
                                                                                    ggTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGCTAAGGCGGTATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGa                                                                  <  2:1791134/150‑1 (MQ=255)
                                                                                     gTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGACCGCGTAGGGCATGCCGCGCAt                                                                                                                                               >  1:1978784/1‑72 (MQ=255)
                                                                                     gTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCAt                                                                                                                                               <  2:1978784/72‑1 (MQ=255)
                                                                                         acGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAGATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACt                                                             >  2:578005/1‑150 (MQ=255)
                                                                                         acGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGgt                                                                                      <  1:1163954/125‑1 (MQ=255)
                                                                                         acGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGgt                                                                                      >  2:1163954/1‑125 (MQ=255)
                                                                                           gcgAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           >  2:1804826/1‑150 (MQ=255)
                                                                                             gAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttttt                                                         <  2:112928/150‑1 (MQ=255)
                                                                                                tGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGat                                                      <  2:1034044/150‑1 (MQ=255)
                                                                                                    cTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGatcatc                                                  >  2:1621225/1‑150 (MQ=255)
                                                                                                     tCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGAtcatca                                                 >  1:453456/1‑150 (MQ=255)
                                                                                                      cGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAg                                                >  2:1303483/1‑150 (MQ=255)
                                                                                                         gCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGTGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                                             <  1:2038592/150‑1 (MQ=255)
                                                                                                                gCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGccc                                      >  1:893266/1‑150 (MQ=255)
                                                                                                                  gtcgcgcGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGa                                    >  1:2400394/4‑150 (MQ=255)
                                                                                                                   tGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAg                                   <  2:1570289/150‑1 (MQ=255)
                                                                                                                      cgcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTa                                >  2:1578089/1‑150 (MQ=255)
                                                                                                                      cgcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTa                                >  2:1576715/1‑150 (MQ=255)
                                                                                                                      cgcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTa                                >  1:279765/1‑150 (MQ=255)
                                                                                                                       gcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCCAc                                                                                                       <  1:2405748/78‑1 (MQ=255)
                                                                                                                       gcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAAc                                                                                                       >  2:2405748/1‑78 (MQ=255)
                                                                                                                          aaaTGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAgcgc                            >  2:1264967/1‑150 (MQ=255)
                                                                                                                                gAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTg                                          <  1:1094420/130‑1 (MQ=255)
                                                                                                                                gAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTg                                          >  2:1094420/1‑130 (MQ=255)
                                                                                                                                gAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGc                      >  1:168236/1‑150 (MQ=255)
                                                                                                                                gAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATAAGATTGGGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGc                      >  1:24155/1‑150 (MQ=255)
                                                                                                                                 aGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGct                     >  1:1627549/1‑150 (MQ=255)
                                                                                                                                     gAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTTGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAgcgc                            <  2:315280/139‑1 (MQ=255)
                                                                                                                                     gAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAgcgc                            >  1:315280/1‑139 (MQ=255)
                                                                                                                                      aaGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAgg                                                                                                                      >  1:2332531/1‑48 (MQ=255)
                                                                                                                                      aaGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAgg                                                                                                                      <  2:2332531/48‑1 (MQ=255)
                                                                                                                                      aaGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCg                >  1:911449/1‑150 (MQ=255)
                                                                                                                                        ggCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           >  2:2226939/1‑105 (MQ=255)
                                                                                                                                           gTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTGGAACCCCGGTACGTCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATc           >  1:898477/1‑150 (MQ=255)
                                                                                                                                              gggCATGCGCGGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCACAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           >  1:62211/1‑99 (MQ=255)
                                                                                                                                              gggCATGCGCGGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           <  2:62211/99‑1 (MQ=255)
                                                                                                                                              gggCATGCGCGGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGaa                         >  2:1690093/1‑133 (MQ=255)
                                                                                                                                              gggCATGCGCGGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGaa                         <  1:1690093/133‑1 (MQ=255)
                                                                                                                                              gggCATGCGCGGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAAAAACATAGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTc        >  2:1190508/1‑150 (MQ=255)
                                                                                                                                               ggCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                                             <  2:1441233/112‑1 (MQ=255)
                                                                                                                                               ggCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                                             >  1:1441233/1‑112 (MQ=255)
                                                                                                                                                gCATGCCGGTCTTGTCCGCCTGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCATGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGTGTGGCCCGAGGTAGCGCGAATGCTCCTCGACATCGTCTg      <  2:1893408/150‑1 (MQ=255)
                                                                                                                                                 cATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACCTCGGCACGTCGGTGTCGTACCCCGGTACgg                                                                     <  1:65903/86‑1 (MQ=255)
                                                                                                                                                 cATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACgg                                                                     >  2:65903/1‑86 (MQ=255)
                                                                                                                                                   tGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGa   <  2:911449/150‑1 (MQ=255)
                                                                                                                                                    gCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGaa  <  1:1264967/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTCGTCGGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGAA  >  NC_002947/95982‑96279

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: