Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 96,134:1 +G 8.7% coding (551/1005 nt) PP_0088 ← alkanal monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294796,1341.G 8.7% 400.1 / 32.0 104coding (551/1005 nt)PP_0088alkanal monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (44/51);  new base G (4/5);  total (48/56)
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.54e-01

Only 100 of 106 total aligned reads displayed.

CGGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGAA  >  NC_002947/95988‑96279
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                                                        <  1:2286463/142‑4 (MQ=255)
                caTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCGCGGCATGccctac                                                                                                                                        >  2:2286463/1‑139 (MQ=255)
                 aTGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCa                                                                                                                               >  1:615564/1‑150 (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:1906619/1‑150 (MQ=255)
                           ccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGc                                                                                                                     >  1:1663821/1‑150 (MQ=255)
                           ccAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGc                                                                                                                     <  2:615564/150‑1 (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)
                                    cgtACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCAGCCAGCTGGGCACTGAACAGGCTGGATCCCa                                                                                                            <  2:373103/148‑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                                                                                                       <  2:909024/147‑1 (MQ=255)
                                            ttGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAAc                                                                                                       >  1:909024/1‑147 (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)
                                                                             cGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCGGGGTACgaca                                                                   <  1:1535125/150‑4 (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                                                                                      >  2:1163954/1‑125 (MQ=255)
                                                                                   acGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGgt                                                                                      <  1:1163954/125‑1 (MQ=255)
                                                                                     gcgAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           >  2:1804826/1‑150 (MQ=255)
                                                                                       gAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttttt                                                         <  2:112928/150‑1 (MQ=255)
                                                                                          tGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGat                                                      <  2:1034044/150‑1 (MQ=255)
                                                                                               tCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGAtcatca                                                 >  1:453456/1‑150 (MQ=255)
                                                                                                cGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAg                                                >  2:1303483/1‑150 (MQ=255)
                                                                                                   gCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGTGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                                             <  1:2038592/150‑1 (MQ=255)
                                                                                                   gCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCTCATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGc                                             <  2:1559193/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:1576715/1‑150 (MQ=255)
                                                                                                                cgcgAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTa                                >  2:1578089/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)
                                                                                                                                aaGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAgg                                                                                                                      >  1:2332531/1‑48 (MQ=255)
                                                                                                                                aaGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCg                >  1:911449/1‑150 (MQ=255)
                                                                                                                                  ggCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACttt                                                           <  1:2226939/105‑1 (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)
                                                                                                                                        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)
                                                                                                                                                   |                                                                                                                                                 
CGGTCTGCGCCACCACCATGGGAATGCCCAGCATCACATACGGCTTGTCCAGCGTGGTCGAGGGCTTGAAATGGTTGCGGTACACGCGAATGGCCTCGTGCATGTAGCGTGGCGCGAAATGCGAGGCGAAGGCGTAGGGCATGCCGCGCATGCCCGCCAGCTGGGCACTGAACAGGCTGGATCCCAGCAACCACATCGGCACGTCGGTGTCGTACCCCGGTACGGCGATCACTTTTTGATCATCAGTGCGTGGCCCGAGGTAGCGCGAAAGCTCCTCGACATCGTCTGGGAA  >  NC_002947/95988‑96279

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