Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,543,877 G→A P237P (CCG→CCA rpnA → recombination‑promoting nuclease RpnA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,543,8770GA97.7% 149.4 / ‑3.9 44P237P (CCG→CCArpnArecombination‑promoting nuclease RpnA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (20/23);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTC  >  NC_000913/3543631‑3544115
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
gCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAac                                                                                                                                                                                                                                             <  1:144609/251‑1 (MQ=255)
        atGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGgag                                                                                                                                                                                                                                     <  1:416490/251‑1 (MQ=255)
                     gCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACt                                                                                                                                                                                                                         <  2:66719/250‑1 (MQ=255)
                              agaAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGaaa                                                                                                                                                                                                               <  2:360240/250‑1 (MQ=255)
                                  aTCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGAtt                                                                                                                                                                                                           <  2:315840/251‑1 (MQ=255)
                                            aCCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTAACCAAAaca                                                                                                                                                                                                                                            >  1:429141/1‑208 (MQ=255)
                                             ccGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTTGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCaggagg                                                                                                                                                                                                <  2:195760/251‑1 (MQ=255)
                                                    aTGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAAc                                                                                                                                                                                          >  1:246700/1‑250 (MQ=255)
                                                                gAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCt                                                                                                                                                                             <  1:374311/251‑1 (MQ=255)
                                                                  gTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGc                                                                                                                                                                           >  1:83408/1‑251 (MQ=255)
                                                                       ttCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCatat                                                                                                                                                                       >  1:287592/1‑250 (MQ=255)
                                                                        tCAAAAACATCTTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAGGTAGAGGATACGCTAATGACACACAGATCAAAAGGCTGTTTCATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGTTTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGGGGGGGAACAATCCAAAGCCCTGCatat                                                                                                                                                                       >  1:62585/1‑249 (MQ=255)
                                                                                      gTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGTCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGGGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCtt                                                                                                                                                       >  1:405967/1‑251 (MQ=255)
                                                                                      gTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCCGGAGGGGGAACAATCCCAAGACCTGCATATAGCCAAAATAATGCtt                                                                                                                                                       >  1:71573/1‑251 (MQ=255)
                                                                                             cgACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCAGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATcct                                                                                                                                                <  2:145757/251‑2 (MQ=255)
                                                                                               aCCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAGCAATACAAAACCCTGCCACTAGCCAAAATAATGATTGaaacccgg                                                                                                                                              >  2:362875/1‑245 (MQ=255)
                                                                                               aCCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGa                                                                                                                                              <  1:193002/251‑1 (MQ=255)
                                                                                                ccTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGTATACTCTAATGACAGACAAATCAAAGGTCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAg                                                                                                                                             <  1:432946/251‑1 (MQ=255)
                                                                                                  tGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCata                                                                                                                                                                        <  1:348577/222‑1 (MQ=255)
                                                                                                               agagCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGTAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGAc                                                                                                                               <  1:38761/250‑1 (MQ=255)
                                                                                                                       tGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATgcgc                                                                                                                      >  1:452563/1‑251 (MQ=255)
                                                                                                                           cTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGAAATCATGCGCttt                                                                                                                   >  1:286231/1‑250 (MQ=255)
                                                                                                                                tattaAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCgg                                                                                                              <  1:213005/250‑1 (MQ=255)
                                                                                                                                 attaAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTAACGGGc                                                                                                            >  1:249042/1‑251 (MQ=255)
                                                                                                                                    acgtagtGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGACGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCTGCAGGAGGGGGAACAATCCAAAGCCCGGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTGGCAGACATCATGCGCTTTACCGGGCTGt                                                                                                         <  2:237911/249‑1 (MQ=255)
                                                                                                                                       tagtGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTc                                                                                                        >  1:100056/1‑249 (MQ=255)
                                                                                                                                           ggATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCtt                                                                                                                                     <  1:235876/216‑1 (MQ=255)
                                                                                                                                                      tGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTg                                                                                       >  1:434742/1‑251 (MQ=255)
                                                                                                                                                       gacagacaAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTgc                                                                                      <  1:267619/251‑1 (MQ=255)
                                                                                                                                                              aaaTCAAAGGGCTGTTTAATTACATACTGCAAACTAGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGCGGGGGAACAATCCAAAACCCCGCATATAGACAAAATAATACTTGAATCCGGAATTCCTCTTGCAGAAATCATGCGCCTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAgg                                                                               >  2:11886/1‑250 (MQ=255)
                                                                                                                                                                 tCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAgggg                                                                                                                                                                                              >  2:288944/1‑137 (MQ=255)
                                                                                                                                                                 tCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAgggg                                                                                                                                                                                              <  1:288944/137‑1 (MQ=255)
                                                                                                                                                                    aaGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGTTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTaa                                                                         <  1:355769/251‑1 (MQ=255)
                                                                                                                                                                     aGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTGAACGATGTTATCGACGGCGTTGCCGAACGTGCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATGGCGGCAGGAGGGGGAACAAGCCAAAGCCCTGCATATAGCCCAAATAATGCTTg                                                                                                                                                      <  2:340092/173‑1 (MQ=255)
                                                                                                                                                                     aGGGCTGTTTAATTACATACTGCAAACTAGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAACCCTGCATCTAGCCAAAATAATGCTTg                                                                                                                                                      >  1:340092/1‑173 (MQ=255)
                                                                                                                                                                                   aCATTCTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCTGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGcc                                                           <  2:338720/250‑1 (MQ=255)
                                                                                                                                                                                   aCATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGcc                                                           <  1:421885/250‑1 (MQ=255)
                                                                                                                                                                                    cATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCAt                                                         >  1:400099/1‑251 (MQ=255)
                                                                                                                                                                                           cAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAACCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAAACATCATGTGCTTTACCGGGCTGTCAGCATAAAAGTTGGCTTCGGCGAATCCGTATAATTCTGTCTCCCCATTTCaac                                                   >  2:350619/1‑249 (MQ=255)
                                                                                                                                                                                               cTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGGCAGTAAAGTTCTGTCTCGCCATTTCAAAAGcc                                               >  2:222163/1‑250 (MQ=255)
                                                                                                                                                                                                                      gATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGcc                      <  1:273416/251‑1 (MQ=255)
                                                                                                                                                                                                                            aTCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACcag                 <  1:125133/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTGGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGACCCCTACACCCTACTCTTCAACGCCACCAGCCGGTGACa        >  2:322933/1‑250 (MQ=255)
                                                                                                                                                                                                                                          cGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGAGACAAAACTc  >  1:54988/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
GCTGATTGATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTC  >  NC_000913/3543631‑3544115

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: