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,8770GA89.8% 136.1 / 7.7 49P237P (CCG→CCArpnArecombination‑promoting nuclease RpnA
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (20/24);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGA  >  NC_000913/3543638‑3544122
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gatGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCAACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGga                                                                                                                                                                                                                                            <  1:148705/251‑1 (MQ=255)
       ccGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGttt                                                                                                                                                                                                                                      >  1:242396/1‑250 (MQ=255)
           cATGCCTTGTGATGAAATCATGGAGCACCGACGTATGGCGCTGCTGTAGTTAATTCAAAAACATATTCGGCAACGCGACCTGAATGGGATGGTAGAGCAAATTTCCCGCTTATTAAGTAGTGGGTACGCTAATGACAGACAAATCAAAGGTCTGTTTAATTACATACTGCAAACTGGCGACGCTCTACGTTATAACTATTTTATCGACGGCGTTGCCGAACGTTCACCAAATTACAAGGAGAGTTTAATgt                                                                                                                                                                                                                                 <  2:50604/251‑2 (MQ=255)
                  gatgatgaAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACAGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGTCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTg                                                                                                                                                                                                                           <  1:373142/250‑1 (MQ=255)
                   atgatgaAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCg                                                                                                                                                                                                                         >  1:196276/1‑251 (MQ=255)
                             cATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGc                                                                                                                                                                                                               <  2:447153/251‑1 (MQ=255)
                              aTGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCg                                                                                                                                                                                                              >  2:365993/1‑251 (MQ=255)
                                gcagcaCCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTTTTAAGTAGTGGATACGCTAATAACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGc                                                                                                                                                                                                            >  1:28379/1‑251 (MQ=255)
                                  agcaCCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAAcaca                                                                                                                                                                                                                                                <  1:499318/213‑1 (MQ=255)
                                  agcaCCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAAcaca                                                                                                                                                                                                                                                >  2:499318/1‑213 (MQ=255)
                                  agcaCCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACCCAAGGCGAGTTTAATGACTATTGCGGAAAGAtttccgca                                                                                                                                                                                                           >  2:511195/1‑244 (MQ=255)
                                  agcaCCGACGGATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCag                                                                                                                                                                                                          <  2:41942/251‑1 (MQ=255)
                                         aCGTATGGCGTTTCTGGATTTAATTCTAAAACATATTCGTCAACGCGTCCTGATGGGTCTGGTACAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTTCGTTTTAACTATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAATGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAggggg                                                                                                                                                                                                   <  1:211167/251‑1 (MQ=255)
                                                    tgctgGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGGGGGGGAACAATCCaa                                                                                                                                                                                         >  2:532761/1‑250 (MQ=255)
                                                         gAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCt                                                                                                                                                                                   <  1:149419/251‑1 (MQ=255)
                                                            ttAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGc                                                                                                                                                                                 >  1:332744/1‑250 (MQ=255)
                                                                                    cgcgACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATc                                                                                                                                                        <  2:502751/251‑1 (MQ=255)
                                                                                     gcgACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATcc                                                                                                                                                       <  1:55755/251‑1 (MQ=255)
                                                                                        aCCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGa                                                                                                                                                    <  1:58733/251‑1 (MQ=255)
                                                                                         ccTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCgg                                                                                                                                                     >  1:11819/1‑249 (MQ=255)
                                                                                         ccTGATGGGGCTGGTAGAGCAAATGGACTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAg                                                                                                                                                   <  2:563360/251‑1 (MQ=255)
                                                                                               ggggCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCctc                                                                                                                                             <  1:511022/251‑1 (MQ=255)
                                                                                                     ggTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACTCTAATGACAGACAAAGCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCCCTTGCAg                                                                                                                                       <  2:41797/251‑1 (MQ=255)
                                                                                                        agagCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTGAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGAc                                                                                                                                     <  1:26376/250‑1 (MQ=255)
                                                                                                            cAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCTCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGAcatcat                                                                                                                                >  1:487287/1‑251 (MQ=255)
                                                                                                            cAAATGGCCTGATTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGAcatcat                                                                                                                                <  2:337072/251‑1 (MQ=255)
                                                                                                                    cTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGGACAATCCAAAGCCCAGCATAAAACCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGAtttt                                                                                                                        >  2:536289/1‑250 (MQ=255)
                                                                                                                          attaAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAAAAATCCAAAGCCCTGCATATAGCCAAAATAATTCTTAAATCCGGAGTTCCTCTTGCAAACATCATGCGCTTTACCGGGc                                                                                                                  >  2:217641/1‑251 (MQ=255)
                                                                                                                            taAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTg                                                                                                                >  1:559534/1‑251 (MQ=255)
                                                                                                                             aaGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAgggg                                                                                                                                                                                                    >  1:387618/1‑166 (MQ=255)
                                                                                                                             aaGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAgggg                                                                                                                                                                                                    <  2:387618/166‑1 (MQ=255)
                                                                                                                               gtagtGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGCACAATCCAAAGCCCTGCATATTGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCCGACCACATTCGCTTTACCGGGCTGTc                                                                                                              >  2:460007/1‑250 (MQ=255)
                                                                                                                                         cGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGt                                                                                                   <  1:245472/251‑1 (MQ=255)
                                                                                                                                          gCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGGCAGAAGAAGAGtt                                                                                                  >  1:284541/1‑251 (MQ=255)
                                                                                                                                              aTGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTggtt                                                                                              <  2:335882/251‑3 (MQ=255)
                                                                                                                                               tGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTg                                                                                             >  1:338144/1‑251 (MQ=255)
                                                                                                                                                      caAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTgcggcg                                                                                        >  2:162119/1‑249 (MQ=255)
                                                                                                                                                         aTCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCa                                                                                   <  1:359277/251‑1 (MQ=255)
                                                                                                                                                                     tttAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGCTTTTATCGACGGCGTTGCCGAACGTTCACCAAACCACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAAAAATCCAAAGCCCTGCATATAGACAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCACGCGCTTTACCGGGCTGTCACACGAAGAGTTGGCTGCGGCGAGGCAGTGAAgatcggt                                                                       >  1:342739/1‑245 (MQ=255)
                                                                                                                                                                             cATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCAt                                                               <  1:309839/251‑1 (MQ=255)
                                                                                                                                                                                          ggCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCAc                                                   <  2:284541/250‑1 (MQ=255)
                                                                                                                                                                                                   gTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATTCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACAccc                                           <  2:28379/249‑1 (MQ=255)
                                                                                                                                                                                                   gTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCTCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACAccc                                           <  2:487287/249‑1 (MQ=255)
                                                                                                                                                                                                   gTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCct                                          <  1:532761/250‑1 (MQ=255)
                                                                                                                                                                                                                    tATCGACGGCGTTGCCGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGccacc                          >  1:151014/1‑249 (MQ=255)
                                                                                                                                                                                                                                  ccGAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACt          <  2:338144/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gAACGTTCACCAACACACAAGGAGAGTTTAATGAATATTGCGGAAAGATTGCGGCAGGAGGGGGAAAAATCCAAAGCCCTGCATATAACCAACATAATTCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGTGTTGGCTGCGGAGAGGAAATAAAGTTCTGTCTCCCCATATCAAAAGACACCTACACCCTCTGATTCAACACCACCCGCAGGTGACAAAACTcc        >  2:380272/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gAACGTTCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCg        <  1:33829/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tCACCAAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCCCTGCTTCAACGCCACCAGCAGGGGGCAAAACTCGGCCggg  >  2:161684/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GATGTCACCGTCATGCCTGATGATGAAATCATGCAGCACCGACGTATGGCGCTGCTGGAGTTAATTCAAAAACATATTCGTCAACGCGACCTGATGGGGCTGGTAGAGCAAATGGCCTGCTTATTAAGTAGTGGATACGCTAATGACAGACAAATCAAAGGGCTGTTTAATTACATACTGCAAACTGGCGACGCTGTACGTTTTAACGATTTTATCGACGGCGTTGCCGAACGTTCACCGAAACACAAGGAGAGTTTAATGACTATTGCGGAAAGATTGCGGCAGGAGGGGGAACAATCCAAAGCCCTGCATATAGCCAAAATAATGCTTGAATCCGGAGTTCCTCTTGCAGACATCATGCGCTTTACCGGGCTGTCAGAAGAAGAGTTGGCTGCGGCGAGTCAGTAAAGTTCTGTCTCGCCATTTCAAAAGCCACCTACACCCTCTGCTTCAACGCCACCAGCAGGTGACAAAACTCGGCCGGA  >  NC_000913/3543638‑3544122

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