Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,859,142 G→A 100% A684T (GCC→ACC)  mutS → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,859,1420GA97.8% 2660.8 / 23.3 736A684T (GCC→ACC) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (4/12);  new base A (382/335);  total (386/350)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 743 total aligned reads displayed.

TACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTG  >  NC_000913/2858995‑2859288
                                                                                                                                                   |                                                                                                                                                  
 aCATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCa                                                                                                                                                 >  1:9478939/1‑150 (MQ=255)
 aCATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCa                                                                                                                                                 >  1:9302136/1‑150 (MQ=255)
   aTCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACaccacc                                                                                                                                               <  2:6110709/150‑1 (MQ=255)
      ggCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGaa                                                                                                                                            <  2:10531085/150‑1 (MQ=255)
      ggCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGaa                                                                                                                                            <  1:11898245/150‑1 (MQ=255)
      ggCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGaa                                                                                                                                            >  1:2164913/1‑150 (MQ=255)
       gCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAAt                                                                                                                                           >  2:3539664/1‑150 (MQ=255)
         aGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATAc                                                                                                                                         >  2:7976079/1‑150 (MQ=255)
          gCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACa                                                                                                                                        <  2:10307209/150‑1 (MQ=255)
             aTGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTc                                                                                                                                     >  2:4459673/1‑150 (MQ=255)
               gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTg                                                                                                                                   <  2:2510381/150‑1 (MQ=255)
               gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTg                                                                                                                                   >  2:4102951/1‑150 (MQ=255)
                  ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGgtg                                                                                                                                >  1:4753883/1‑150 (MQ=255)
                   cGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGgtgt                                                                                                                               >  1:3114668/1‑150 (MQ=255)
                   cGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGgtgt                                                                                                                               >  1:6442273/1‑150 (MQ=255)
                     gCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTa                                                                                                                             >  2:12040700/1‑150 (MQ=255)
                         aaaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATgg                                                                                                                         <  1:1383081/150‑1 (MQ=255)
                         aaaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATgg                                                                                                                         >  1:9040680/1‑150 (MQ=255)
                          aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGa                                                                                                                        <  2:11883857/150‑1 (MQ=255)
                          aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGa                                                                                                                        >  2:5234789/1‑150 (MQ=255)
                                 gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                 <  1:7283832/150‑1 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                          >  1:11162335/1‑150 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                          <  1:10004129/150‑1 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                          <  1:5288384/150‑1 (MQ=255)
                                         aCCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGAGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTCATGGATGAGATCGGGCGTgg                                                                                                         >  1:5550808/1‑150 (MQ=255)
                                           cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa                                                                                                       >  2:5877329/1‑150 (MQ=255)
                                           cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa                                                                                                       >  2:8397645/1‑150 (MQ=255)
                                             atcgatcgCATCTTTACCCGCGGAGGCGCGGCAGAGGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGCATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACg                                                                                                     <  2:8557032/150‑1 (MQ=255)
                                                gatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                                  <  2:8272256/150‑1 (MQ=255)
                                                gatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                                  >  1:10105441/1‑150 (MQ=255)
                                                gatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACACAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                                  >  1:11398435/1‑150 (MQ=255)
                                                  tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                                >  1:10622957/1‑150 (MQ=255)
                                                      aTCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTAc                                                                                            >  1:8358896/1‑150 (MQ=255)
                                                       tCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACg                                                                                           >  1:9768906/1‑150 (MQ=255)
                                                         tttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGAt                                                                                         >  2:11225666/1‑150 (MQ=255)
                                                         tttACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGGTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGAt                                                                                         >  1:8885148/1‑150 (MQ=255)
                                                           tACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATgg                                                                                       <  2:1682142/150‑1 (MQ=255)
                                                              ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt                                                                                    >  2:3663735/1‑150 (MQ=255)
                                                              ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt                                                                                    >  2:2105677/1‑150 (MQ=255)
                                                              ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCGTTATGGTGGAGATGACTGAAACCGCCAAGATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGAGGGTcg                                                                                    >  2:6849891/1‑149 (MQ=255)
                                                                gcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGt                                                                                  >  2:1422891/1‑150 (MQ=255)
                                                                 cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                 <  2:1214789/150‑1 (MQ=255)
                                                                 cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                 <  1:5950128/150‑1 (MQ=255)
                                                                  gTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCg                                                                                <  2:7238766/150‑1 (MQ=255)
                                                                    aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCt                                                                              <  2:5544404/150‑1 (MQ=255)
                                                                    aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCt                                                                              <  2:4818881/150‑1 (MQ=255)
                                                                    aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCt                                                                              <  1:6072058/150‑1 (MQ=255)
                                                                     ggCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTg                                                                             <  2:2809680/150‑1 (MQ=255)
                                                                      gcgcGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGTCTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg                                                                            <  2:9687958/150‑1 (MQ=255)
                                                                      gcgcGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg                                                                            >  2:6486396/1‑150 (MQ=255)
                                                                      gcgcGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg                                                                            <  1:6214286/150‑1 (MQ=255)
                                                                      gcgcGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAAGACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTgg                                                                            <  1:10349101/150‑1 (MQ=255)
                                                                         cgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgt                                                                         <  1:4143569/150‑1 (MQ=255)
                                                                          ggCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtg                                                                        >  1:3064419/1‑150 (MQ=255)
                                                                          ggCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtg                                                                        <  1:1302276/150‑1 (MQ=255)
                                                                             aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggc                                                                     >  2:10142016/1‑150 (MQ=255)
                                                                             aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggc                                                                     >  1:12141198/1‑150 (MQ=255)
                                                                             aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggc                                                                     >  1:657953/1‑150 (MQ=255)
                                                                             aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgggggc                                                                     >  1:5161902/1‑150 (MQ=255)
                                                                               aTGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcgt                                                                   <  1:1744786/150‑1 (MQ=255)
                                                                               aTGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcgt                                                                   <  1:710782/150‑1 (MQ=255)
                                                                                   ccTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTgcgc                                                               <  2:11376469/150‑1 (MQ=255)
                                                                                    cTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGcgcg                                                              <  1:8716380/150‑1 (MQ=255)
                                                                                     tGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCgg                                                             >  2:2330665/1‑150 (MQ=255)
                                                                                      ggCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACAACACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGa                                                            >  1:3772007/1‑150 (MQ=255)
                                                                                       gCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaa                                                           >  1:8189554/1‑150 (MQ=255)
                                                                                         gTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGaaaa                                                         <  2:7450578/150‑1 (MQ=255)
                                                                                            cGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCt                                                      <  2:328747/150‑1 (MQ=255)
                                                                                             gggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTg                                                     <  2:6230359/150‑1 (MQ=255)
                                                                                              ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                    >  1:8022551/1‑150 (MQ=255)
                                                                                               gcgcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGc                                                   <  2:3983954/150‑1 (MQ=255)
                                                                                               gcgcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGc                                                   <  1:2470056/150‑1 (MQ=255)
                                                                                                   tCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAAt                                               >  2:11476460/1‑150 (MQ=255)
                                                                                                   tCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAAt                                               >  1:788702/1‑150 (MQ=255)
                                                                                                   tCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAAt                                               >  1:10810072/1‑150 (MQ=255)
                                                                                                     aaCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATaa                                             >  1:6877732/1‑150 (MQ=255)
                                                                                                      aCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAg                                            <  1:5730800/150‑1 (MQ=255)
                                                                                                      aCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAg                                            >  2:5748497/1‑150 (MQ=255)
                                                                                                       ccTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGa                                           <  1:7397383/150‑1 (MQ=255)
                                                                                                            aTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAGTCTGGCGAATAAGATTAAg                                      <  2:9472600/150‑1 (MQ=255)
                                                                                                             tggtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAgg                                     <  2:5385096/150‑1 (MQ=255)
                                                                                                               gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                   >  2:3019733/1‑150 (MQ=255)
                                                                                                                          tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                        >  1:991642/1‑150 (MQ=255)
                                                                                                                          tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                        >  1:2699602/1‑150 (MQ=255)
                                                                                                                          tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                        >  2:10777566/1‑150 (MQ=255)
                                                                                                                            aaaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCTCCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTg                      <  1:11707239/150‑1 (MQ=255)
                                                                                                                            aaaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTg                      <  2:920182/150‑1 (MQ=255)
                                                                                                                            aaaCCGCCAATATTTTACATAACACCACCGAAAACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCAACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCAATGACGTTATTTg                      >  2:8791675/1‑150 (MQ=255)
                                                                                                                               ccGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTa                   <  1:436501/150‑1 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                 >  2:6046378/1‑150 (MQ=255)
                                                                                                                                    aaTATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc              >  2:11663764/1‑150 (MQ=255)
                                                                                                                                     atatTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACt             <  1:11553193/150‑1 (MQ=255)
                                                                                                                                         tttACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTAttt         <  1:4460669/150‑1 (MQ=255)
                                                                                                                                          ttACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTc        <  1:8474167/150‑1 (MQ=255)
                                                                                                                                            aCATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGa      <  2:7742726/150‑1 (MQ=255)
                                                                                                                                            aCATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGa      >  2:8971726/1‑150 (MQ=255)
                                                                                                                                             cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAg     <  2:788969/150‑1 (MQ=255)
                                                                                                                                             cATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAg     <  1:7673666/150‑1 (MQ=255)
                                                                                                                                               tAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCt   <  2:11031616/150‑1 (MQ=255)
                                                                                                                                               tAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCt   >  1:3466275/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCTG  >  NC_000913/2858995‑2859288

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