Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,388,285 G→A 8.0% A281V (GCT→GTT) ‡ aaeA ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,2850GA 8.0% 376.9 / 13.9 124A281V (GCT→GTT) ‡aaeAp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base G (55/59);  new base A (4/6);  total (59/65)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 127 total aligned reads displayed.

GCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCT  >  NC_000913/3388136‑3388434
                                                                                                                                                     |                                                                                                                                                     
gCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAg                                                                                                                                                       >  1:1384638/1‑150 (MQ=255)
 ccAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGc                                                                                                                                                      <  1:1126467/150‑1 (MQ=255)
  cAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCa                                                                                                                                                     >  1:164633/1‑150 (MQ=255)
    gTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAgt                                                                                                                                                   <  2:767977/150‑1 (MQ=255)
       ttACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGTCCATTTTACGGATGAACGAATCCAGGCTATCGTCGCGATCTTGTTTGCAAGTGACCACCAATGAAGCAgtggt                                                                                                                                                >  2:214760/1‑150 (MQ=255)
        tACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAgtggtg                                                                                                                                               <  1:22814/150‑1 (MQ=255)
         aCCGCTAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtgg                                                                                                                                              <  2:1459244/150‑2 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGc                                                                                                                                              >  1:704401/1‑150 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGc                                                                                                                                              <  2:650287/150‑1 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtgg                                                                                                                                              <  1:1018411/150‑2 (MQ=255)
         aCCGCAAAGCGAATATGTTGGTTAGCAATCGAGAAAATACCCATCGAGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCCCCACTGCTACAgtggtgg                                                                                                                                              <  1:1125900/150‑2 (MQ=255)
          ccGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCTCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtggt                                                                                                                                             <  2:959865/150‑3 (MQ=255)
           cGCAAAGCGAAGATGTTGGTTAGCATTGGAGAATATACCCATCGTGATTAACCTACCTCTCGCAGGCGAGGGGCCATTTTACGGAAGAACGAGTCCTGGCTTTCGTCGCGATCTTGGGTGCCAGTGAACACCACTGGAGCAGTGGTGccc                                                                                                                                            >  1:290103/1‑150 (MQ=255)
            gCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCg                                                                                                                                           <  1:823986/150‑1 (MQ=255)
            gCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGCTACAgtggtggtca                                                                                                                                           <  1:283655/150‑5 (MQ=255)
             cAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGc                                                                                                                                          >  1:429474/1‑150 (MQ=255)
             cAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGc                                                                                                                                          <  1:1279379/150‑1 (MQ=255)
                aGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCCCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGGCGCGATCTTGTTTGACAGTGACCACCACTGTAGCAGTGGTGCCCGCAgg                                                                                                                                       <  1:1200703/150‑1 (MQ=255)
                aGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAgg                                                                                                                                       >  2:867855/1‑150 (MQ=255)
                aGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAgg                                                                                                                                       <  1:1377361/150‑1 (MQ=255)
                 gCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGc                                                                                                                                      >  1:853308/1‑150 (MQ=255)
                 gCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAATTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGc                                                                                                                                      >  2:1018277/1‑150 (MQ=255)
                   gAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCa                                                                                                                                    <  1:827089/150‑1 (MQ=255)
                       atGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGATGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATg                                                                                                                                >  2:1009077/1‑150 (MQ=255)
                        tGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGt                                                                                                                               >  1:754746/1‑150 (MQ=255)
                          ttGGTTAGCAATGGAGAAAATACACATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCAAGTGCCCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtc                                                                                                                             >  1:1080307/1‑150 (MQ=255)
                            ggTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtctc                                                                                                                           >  1:5042/1‑150 (MQ=255)
                            ggTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTtctc                                                                                                                           >  1:1213894/1‑150 (MQ=255)
                             gTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTctct                                                                                                                          <  2:208853/150‑1 (MQ=255)
                                 gCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCt                                                                                                                      <  1:211294/150‑1 (MQ=255)
                                  cAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCATGCTg                                                                                                                     >  2:49504/1‑150 (MQ=255)
                                  cAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAATTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTg                                                                                                                     <  1:1018277/150‑1 (MQ=255)
                                   aaTGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTgg                                                                                                                    <  1:1485642/150‑1 (MQ=255)
                                      ggAGAAAATACCCATCGTGATTAACCAAACTCCCGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTTGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTg                                                                                                                 <  1:1412713/150‑1 (MQ=255)
                                       gagaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGt                                                                                                                >  2:303101/1‑150 (MQ=255)
                                         gaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCg                                                                                                              <  2:754746/150‑1 (MQ=255)
                                         gaAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCg                                                                                                              <  2:164633/150‑1 (MQ=255)
                                            aaTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCgaga                                                                                                           <  2:5042/150‑1 (MQ=255)
                                             aTACCCATCGTGATTAACCAAACTCACGCAGGCTATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGAc                                                                                                          <  2:318936/150‑1 (MQ=255)
                                             aTACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGAGGTTCTCTTGCTGGTTGTCGAGAc                                                                                                          <  1:113138/150‑1 (MQ=255)
                                                cccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                       >  1:66576/1‑150 (MQ=255)
                                                 ccATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAAt                                                                                                      >  2:422627/1‑150 (MQ=255)
                                                  cATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTTGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                     <  2:580947/150‑1 (MQ=255)
                                                  cATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATa                                                                                                     <  2:1585045/150‑1 (MQ=255)
                                                          ttAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCgg                                                                                             <  2:847149/150‑1 (MQ=255)
                                                             aCCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                          >  1:1114189/1‑150 (MQ=255)
                                                             aCCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAac                                                                                          <  1:130935/150‑1 (MQ=255)
                                                               cAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAacac                                                                                        >  2:554409/1‑150 (MQ=255)
                                                                aaaCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACg                                                                                       <  1:303101/150‑1 (MQ=255)
                                                                   cTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTg                                                                                    >  1:587109/1‑150 (MQ=255)
                                                                   cTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTg                                                                                    >  1:580479/1‑150 (MQ=255)
                                                                     cACGCAGGCGATGGGCCATTTTACGCAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGAAGGCCAGATGTACTCATGCTGGTAGTCGAGACGAATACGAACCGGTACACGTTgcg                                                                                  >  1:885235/1‑150 (MQ=255)
                                                                        gCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACGGTAGCAGTGGTGCCCGCAGGCCAGATGTTATCTTGCTGGTTGTCGAGACGAATTCGACCCGGACCCCGTTACGCaa                                                                               >  2:79842/1‑150 (MQ=255)
                                                                         cAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAg                                                                              <  1:1656510/150‑1 (MQ=255)
                                                                         cAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGGGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAg                                                                              <  2:486948/150‑1 (MQ=255)
                                                                                tGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCTATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCAcc                                                                       <  2:847876/150‑1 (MQ=255)
                                                                                tGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCAcc                                                                       <  2:1513335/150‑1 (MQ=255)
                                                                                   gCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCaa                                                                    >  1:776359/1‑149 (MQ=255)
                                                                                     cATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTAGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGCCGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTc                                                                  >  2:142402/1‑150 (MQ=255)
                                                                                       ttttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCaa                                                                <  1:981613/150‑1 (MQ=255)
                                                                                         ttACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAgg                                                              >  1:1483443/1‑150 (MQ=255)
                                                                                            tggAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACCCGTTGCGCAAGACGCACCCAGTCACGGTTa                                                           >  1:1549695/2‑150 (MQ=255)
                                                                                              gaagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTaga                                                         <  2:915654/150‑1 (MQ=255)
                                                                                              gaagaaCGAATCCTGGCTTTCGTCGCCATCTTGTTTGCCAGTGACCACCAGTGTAGCAGAGGTGCCCGCAGGCCAGATGTTATCTTGCTGGTTGTCGAGCCGAATACGAACCGGAACACGTTGCGCAAGACCCCCCCATTCAAGGTTtga                                                         >  2:1228110/1‑150 (MQ=255)
                                                                                               aagaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTagag                                                        <  2:806249/150‑1 (MQ=255)
                                                                                                agaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGt                                                       <  1:418319/150‑1 (MQ=255)
                                                                                                 gaaCGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTc                                                      <  1:1450158/150‑1 (MQ=255)
                                                                                                    cGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCtat                                                   <  2:66576/150‑1 (MQ=255)
                                                                                                           tGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGcc                                            <  2:853308/150‑1 (MQ=255)
                                                                                                             gCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCAt                                          >  2:556694/1‑150 (MQ=255)
                                                                                                               tttCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGCCGAATACGAACCGGAACACGTTGCGCCAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATcc                                        >  2:1665402/1‑150 (MQ=255)
                                                                                                                  cgtcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCtt                                     <  2:100871/150‑1 (MQ=255)
                                                                                                                   gtcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCttt                                    >  1:1666161/1‑150 (MQ=255)
                                                                                                                    tcgCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGGTAGAGTCTATAGTCGCCATCCCTTTg                                   <  1:867855/150‑1 (MQ=255)
                                                                                                                      gcgATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCAGTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                                 >  1:1681190/1‑150 (MQ=255)
                                                                                                                       cgATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTAAGAGTCTATAGTCGCCATCCCTTTgtcg                                >  1:983510/1‑150 (MQ=255)
                                                                                                                            ttGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAAAACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcg                           <  1:857549/150‑1 (MQ=255)
                                                                                                                            ttGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCATGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcg                           <  1:49504/150‑1 (MQ=255)
                                                                                                                               tttGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGt                        >  2:616958/1‑150 (MQ=255)
                                                                                                                               tttGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGt                        >  2:1697266/1‑150 (MQ=255)
                                                                                                                                     aGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGCACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTgctgct                  >  1:217866/1‑150 (MQ=255)
                                                                                                                                     aGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTgctgct                  >  1:479216/1‑150 (MQ=255)
                                                                                                                                      gTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCCTTCAAGGTTAGAGTCAATAGACGCCATCCCTTTGTCGTCCAGCGTGctgctg                 >  2:238844/1‑150 (MQ=255)
                                                                                                                                      gTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGctgctg                 >  2:337035/1‑150 (MQ=255)
                                                                                                                                       ggACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGGTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTgg                <  2:704401/149‑1 (MQ=255)
                                                                                                                                        gggcaccacTGCTACAGTGGTGCCCGCAGGCCAGATGTGCGCGTGCTGGTTGTCGAGACGAATACGAACCGGACCACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGc               >  2:25427/4‑150 (MQ=255)
                                                                                                                                         ggcaccacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCg              >  2:1311446/3‑150 (MQ=255)
                                                                                                                                         ggcaccacTGCTACAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCg              >  2:684890/3‑150 (MQ=255)
                                                                                                                                         accaccACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCg              >  1:957709/1‑150 (MQ=255)
                                                                                                                                           caccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt            <  1:554409/150‑1 (MQ=255)
                                                                                                                                           caccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt            >  2:997110/1‑150 (MQ=255)
                                                                                                                                           caccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGGGCTGCTGGCGtt            <  1:694763/150‑1 (MQ=255)
                                                                                                                                           caccacTGTAGCAGTGGTGCCCGCAGGCCACATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt            >  2:1339825/1‑150 (MQ=255)
                                                                                                                                            accacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTg           <  2:217866/150‑1 (MQ=255)
                                                                                                                                             ccacTGTAGCAGTGGTGCCCGCAGTCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg          >  1:58713/1‑150 (MQ=255)
                                                                                                                                             ccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg          <  2:45511/150‑1 (MQ=255)
                                                                                                                                             ccacTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTgg          <  2:1213894/150‑1 (MQ=255)
                                                                                                                                                 tGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGAc      <  2:1080307/150‑1 (MQ=255)
                                                                                                                                                   tAGCAGGGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGAccc    <  1:1100861/150‑1 (MQ=255)
                                                                                                                                                     gCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCt  >  2:740404/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATACGAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCT  >  NC_000913/3388136‑3388434

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: