Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,859,142 G→A 68.9% 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,1420GA68.9% 266.8 / 182.1 200A684T (GCC→ACC) mutSmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (75/61);  total (103/97)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 200 total aligned reads displayed.

TACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCT  >  NC_000913/2858995‑2859287
                                                                                                                                                   |                                                                                                                                                 
 aCATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCa                                                                                                                                                >  2:2223978/1‑150 (MQ=255)
      ggCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGaa                                                                                                                                           >  2:966943/1‑150 (MQ=255)
      ggCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCGTTATGGTGGAGATGACTGAAACCGCCCATATTTTACATAACACCACCGaa                                                                                                                                           >  2:1362792/1‑150 (MQ=255)
       gCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAAt                                                                                                                                          >  1:1331338/1‑150 (MQ=255)
            tATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGt                                                                                                                                     >  2:501765/1‑150 (MQ=255)
             aTGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTc                                                                                                                                    <  2:238541/150‑1 (MQ=255)
              tGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCt                                                                                                                                   >  1:820976/1‑150 (MQ=255)
              tGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCt                                                                                                                                   >  1:1652943/1‑150 (MQ=255)
               gTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTg                                                                                                                                  >  1:2560981/1‑150 (MQ=255)
                  ccGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGgtg                                                                                                                               >  1:2621535/1‑150 (MQ=255)
                     gCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTa                                                                                                                            <  1:1388874/150‑1 (MQ=255)
                      cacaAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa                                                                                                                           <  2:1525411/150‑1 (MQ=255)
                      cacaAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTaa                                                                                                                           >  1:839993/1‑150 (MQ=255)
                       acaAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAAt                                                                                                                          <  1:1570628/150‑1 (MQ=255)
                          aaaaGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGa                                                                                                                       <  1:2211911/150‑1 (MQ=255)
                              gTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATgag                                                                                                                   >  1:2811904/1‑150 (MQ=255)
                                 gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                <  1:889215/150‑1 (MQ=255)
                                 gagaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATc                                                                                                                <  1:1750074/150‑1 (MQ=255)
                                  agaTTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCg                                                                                                               >  1:510070/1‑150 (MQ=255)
                                       ggACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGt                                                                                                          <  2:2091060/150‑1 (MQ=255)
                                       ggACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGt                                                                                                          <  2:1152799/150‑1 (MQ=255)
                                       ggACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGt                                                                                                          <  1:703153/150‑1 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                         <  1:1439158/150‑1 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                         <  2:1294353/150‑1 (MQ=255)
                                        gACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTg                                                                                                         <  1:2423221/150‑1 (MQ=255)
                                           cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGac                                                                                                      >  2:1313059/1‑149 (MQ=255)
                                           cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa                                                                                                      >  2:1032366/1‑150 (MQ=255)
                                           cTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGaa                                                                                                      >  2:229437/1‑150 (MQ=255)
                                            tATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAAc                                                                                                     >  2:407928/1‑150 (MQ=255)
                                             atcgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCCATATTTTACATAACACCACCGAATCCAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACg                                                                                                    >  1:1869828/1‑150 (MQ=255)
                                               cgatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACCCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTc                                                                                                  <  1:1501448/150‑1 (MQ=255)
                                                gatcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTcc                                                                                                 <  2:63802/150‑1 (MQ=255)
                                                 atcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTCTGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCa                                                                                                <  1:3106560/150‑1 (MQ=255)
                                                  tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                               >  1:1573120/1‑150 (MQ=255)
                                                  tcgCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTAGGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCAc                                                                                               <  1:1377479/150‑1 (MQ=255)
                                                    gCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt                                                                                             >  1:957638/1‑150 (MQ=255)
                                                    gCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCt                                                                                             >  1:721670/1‑150 (MQ=255)
                                                     cATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTa                                                                                            >  1:1582320/1‑150 (MQ=255)
                                                     cATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTa                                                                                            >  1:1790292/1‑150 (MQ=255)
                                                      aTCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCCCCTAc                                                                                           >  1:876223/1‑150 (MQ=255)
                                                            aCCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGt                                                                                     <  1:1701835/150‑1 (MQ=255)
                                                             cccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCCCCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTc                                                                                    <  1:1365674/150‑1 (MQ=255)
                                                              ccGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCt                                                                                   >  1:423628/1‑150 (MQ=255)
                                                               cgcgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTg                                                                                  >  2:223514/1‑150 (MQ=255)
                                                                 cgTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTc                                                                                <  1:1226719/150‑1 (MQ=255)
                                                                    aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCt                                                                             <  2:2210354/150‑1 (MQ=255)
                                                                    aGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTGATGGTGGAGATGACTGAAACCGCCAATATTTCACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCt                                                                             <  2:1100405/150‑1 (MQ=255)
                                                                     ggCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTg                                                                            <  2:53503/150‑1 (MQ=255)
                                                                     ggCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATCCAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTg                                                                            >  2:1794148/1‑150 (MQ=255)
                                                                        gcgGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcg                                                                         >  1:2259418/1‑150 (MQ=255)
                                                                           gCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgg                                                                      <  2:148341/150‑1 (MQ=255)
                                                                             aGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggagtgggc                                                                    >  1:2152415/1‑150 (MQ=255)
                                                                              gATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcg                                                                   >  2:2823353/1‑150 (MQ=255)
                                                                              gATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTggcgtgggcg                                                                   >  1:550924/1‑150 (MQ=255)
                                                                                    cTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGcgcg                                                             <  1:1485914/150‑1 (MQ=255)
                                                                                           ccGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATc                                                      <  2:2710176/150‑1 (MQ=255)
                                                                                             gggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTg                                                    >  1:471417/1‑150 (MQ=255)
                                                                                             gggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTg                                                    <  1:2320904/150‑1 (MQ=255)
                                                                                              ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                   <  2:1863730/150‑1 (MQ=255)
                                                                                              ggCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTgg                                                   <  2:3100390/150‑1 (MQ=255)
                                                                                                 gcTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGa                                                <  1:323079/150‑1 (MQ=255)
                                                                                                  cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                               >  2:2681622/1‑150 (MQ=255)
                                                                                                  cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                               <  1:1049664/150‑1 (MQ=255)
                                                                                                  cTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGaa                                               >  1:1839582/1‑150 (MQ=255)
                                                                                                     aaCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATaa                                            <  2:862451/150‑1 (MQ=255)
                                                                                                     aaCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATaa                                            <  2:2921577/150‑1 (MQ=255)
                                                                                                     aaCCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATaa                                            <  2:2512346/150‑1 (MQ=255)
                                                                                                        cTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAt                                         >  1:2024225/1‑150 (MQ=255)
                                                                                                         tttATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGAtt                                        >  2:1038235/1‑150 (MQ=255)
                                                                                                               gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                  >  2:1925459/1‑150 (MQ=255)
                                                                                                               gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                  >  2:1925524/1‑150 (MQ=255)
                                                                                                               gtggAGATGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCa                                  >  1:2953150/1‑150 (MQ=255)
                                                                                                                   agaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGa                              <  2:731369/150‑1 (MQ=255)
                                                                                                                   agaTGACTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGa                              <  1:2823353/150‑1 (MQ=255)
                                                                                                                    gaTGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGAc                             <  2:1747339/150‑1 (MQ=255)
                                                                                                                        aCTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTa                         >  2:1898428/1‑150 (MQ=255)
                                                                                                                         cTGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAt                        <  1:736702/150‑1 (MQ=255)
                                                                                                                          tGAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAtt                       <  2:1818702/150‑1 (MQ=255)
                                                                                                                           gAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAttt                      <  2:2994018/150‑1 (MQ=255)
                                                                                                                           gAAACCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTAttt                      <  2:41004/150‑1 (MQ=255)
                                                                                                                            aaaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTg                     >  2:489499/1‑150 (MQ=255)
                                                                                                                            aaaCCGCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTg                     >  1:2566450/1‑150 (MQ=255)
                                                                                                                             aaCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                    <  1:2098594/150‑1 (MQ=255)
                                                                                                                             aaCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                    >  2:713075/1‑150 (MQ=255)
                                                                                                                             aaCCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGc                    <  1:1516246/150‑1 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                >  1:1287789/1‑150 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                <  1:2449480/150‑1 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                <  1:2119066/150‑1 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                >  2:2337794/1‑150 (MQ=255)
                                                                                                                                 gCCAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAcc                >  1:1104783/1‑150 (MQ=255)
                                                                                                                                  ccAATATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTAccc               >  2:1233295/1‑150 (MQ=255)
                                                                                                                                    aaTATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc             <  2:3032995/150‑1 (MQ=255)
                                                                                                                                    aaTATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc             >  2:117737/1‑150 (MQ=255)
                                                                                                                                    aaTATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc             <  2:1093772/150‑1 (MQ=255)
                                                                                                                                    aaTATTTTACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCAc             >  1:939542/1‑150 (MQ=255)
                                                                                                                                         tttACATAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTAttt        <  2:2621535/150‑1 (MQ=255)
                                                                                                                                             cATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAg    <  2:2590788/150‑1 (MQ=255)
                                                                                                                                              aTAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGc   >  2:1245851/1‑150 (MQ=255)
                                                                                                                                              aTAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGc   >  2:782253/1‑150 (MQ=255)
                                                                                                                                               tAACACCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCt  >  1:689801/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACATCGGCAGCTATGTACCGGCACAAAAAGTCGAGATTGGACCTATCGATCGCATCTTTACCCGCGTAGGCGCGGCAGATGACCTGGCGTCCGGGCGCTCAACCTTTATGGTGGAGATGACTGAAACCGCCAATATTTTACATAACGCCACCGAATACAGTCTGGTGTTAATGGATGAGATCGGGCGTGGAACGTCCACCTACGATGGTCTGTCGCTGGCGTGGGCGTGCGCGGAAAATCTGGCGAATAAGATTAAGGCATTGACGTTATTTGCTACCCACTATTTCGAGCT  >  NC_000913/2858995‑2859287

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: