Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,581,316 G→A G60E (GGA→GAA)  yhhY → N‑acetyltransferase YhhY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,581,3160GA97.3% 250.1 / NA 75G60E (GGA→GAA) yhhYN‑acetyltransferase YhhY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (36/37);  minor base C (2/0);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGAC  >  NC_000913/3581073‑3581559
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
cATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAcgtc                                                                                                                                                                                                                                              >  1:200587/1‑251 (MQ=255)
    tGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGCCGTgg                                                                                                                                                                                                                                          >  2:295382/1‑251 (MQ=255)
     gACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTAATGATGAAGACGTCGTggg                                                                                                                                                                                                                                         <  1:537609/251‑1 (MQ=255)
       cTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGcc                                                                                                                                                                                                                                       >  1:113119/1‑251 (MQ=255)
       cTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGcc                                                                                                                                                                                                                                       >  1:134161/1‑251 (MQ=255)
                gAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCAcc                                                                                                                                                                                                                               >  1:559361/1‑250 (MQ=255)
                gAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCa                                                                                                                                                                                                                              <  1:492629/251‑1 (MQ=255)
                         tgatgTTGCAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGCTCGTCCCGGCATCAAGCACCTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTg                                                                                                                                                                                                                      >  2:363138/1‑250 (MQ=255)
                         tgatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGAAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGAAAAACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATAGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGc                                                                                                                                                                                                                     >  1:235668/1‑251 (MQ=255)
                          gatgTTGAAAGGGACATTTAACCCAAAGAGGACAAAGGAATGAGCGAGATAGTAATACGCCACGCAGAAACACGGGATTAAGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTACGATCATCTGTGGCAGGAGCGCCTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAAGTCGTGGGCCATCTCACCCTTGACGTGc                                                                                                                                                                                                                     >  2:558477/1‑250 (MQ=255)
                            tgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAAc                                                                                                                                                                                                                  >  1:39439/1‑251 (MQ=255)
                                  aaGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGcc                                                                                                                                                                                                            >  2:370046/1‑251 (MQ=255)
                                  aaGGGACACTTAACCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGcc                                                                                                                                                                                                            <  2:198121/251‑1 (MQ=255)
                                         aTTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgccg                                                                                                                                                                                                     >  2:503795/1‑251 (MQ=255)
                                            tACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGa                                                                                                                                                                                                                          <  1:285908/227‑1 (MQ=255)
                                            tACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGACACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGa                                                                                                                                                                                                                          >  2:285908/1‑227 (MQ=255)
                                             aCCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGt                                                                                                                                                                                                 <  1:197932/251‑1 (MQ=255)
                                               cccAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGAGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTc                                                                                                                                                                                                >  2:446765/1‑250 (MQ=255)
                                                ccAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCATAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGTAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACTGCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCAt                                                                                                                                                                                              <  1:266432/251‑1 (MQ=255)
                                                    agagGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTg                                                                                                                                                                                          <  2:272509/251‑1 (MQ=255)
                                                     gagGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCACGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGc                                                                                                                                                                                         >  2:513124/1‑251 (MQ=255)
                                                                    agtgagATAGTAATACGCCACGCAGAAACCCGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACCGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGCAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCCGGCATGTTGCCGATTTTGGTATCtg                                                                                                                                                                          >  2:587629/1‑251 (MQ=255)
                                                                                aTACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACGCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCg                                                                                                                                                              <  2:123464/251‑1 (MQ=255)
                                                                                 tACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGCCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATGTTGGTATCTGTGTCGActctagt                                                                                                                                                             >  2:158981/1‑248 (MQ=255)
                                                                                   cGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTg                                                                                                                                                           <  2:113119/251‑1 (MQ=255)
                                                                                    gCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTgg                                                                                                                                                          >  1:415166/1‑251 (MQ=255)
                                                                                      cACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGCCTCTCGTTGga                                                                                                                                                         >  1:410477/1‑250 (MQ=255)
                                                                                       ccGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaag                                                                                                                                                       <  2:410785/250‑1 (MQ=255)
                                                                                        cGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGGGTCGACTCTCGTTGgaaga                                                                                                                                                      <  2:385509/251‑1 (MQ=255)
                                                                                          cAGAAACACGGGATTACGAGGCCCTCAGGAAGATTCACGACCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGACCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCAATGACGTGCAACAGCGACCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAAc                                                                                                                                                    >  2:481424/1‑251 (MQ=255)
                                                                                                        tACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATCTGTGGCAGGGGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGAATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAAAAGCGCCCACGCCGAAGGCATGTTGCCGATTTTGGTATCTGTGTCGACTCCCGATGGCAGAACCGCGGAGTCGcc                                                                                                                                        >  2:313540/1‑249 (MQ=255)
                                                                                                         aCGAGGCCATCAGGCAGAGTCACGCCCAGCCGGAGGGGTATTGCAACACACTACAGGTGCCTCATCCTTCCGAGCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGTCGTCGCCAgc                                                                                                                                     <  2:431281/251‑1 (MQ=255)
                                                                                                            aGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGcc                                                                                                                                  >  2:279932/1‑251 (MQ=255)
                                                                                                              gCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCGTCCGATCATATTTGGCTGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGGCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAATAACCGCGGCGTCGCCAGCGACCt                                                                                                                                <  2:123468/251‑1 (MQ=255)
                                                                                                                   cAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                                                                                           <  2:415166/251‑1 (MQ=255)
                                                                                                                            tCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatg                                                                                                                   >  1:411199/1‑250 (MQ=255)
                                                                                                                                   cAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                                                                                           <  1:416746/251‑1 (MQ=255)
                                                                                                                                    aGCCGGAGGTGGATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtg                                                                                                          <  2:487839/251‑1 (MQ=255)
                                                                                                                                       cGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgcc                                                                                                                                                                                                      >  1:100062/1‑156 (MQ=255)
                                                                                                                                       cGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgcc                                                                                                                                                                                                      <  2:100062/156‑1 (MQ=255)
                                                                                                                                        ggAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGga                                                                                                                                                         >  1:321748/1‑200 (MQ=255)
                                                                                                                                        ggAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGga                                                                                                                                                         <  2:321748/200‑1 (MQ=255)
                                                                                                                                         gAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCg                                                                                                      >  2:174970/1‑250 (MQ=255)
                                                                                                                                                    cAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGCCTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACCACTGGGTGc                                                                                          >  2:160743/1‑251 (MQ=255)
                                                                                                                                                          actcccggTGCCTCATCCTTCCAATCATATGTGGCAGGAGCGGCTCGCCGATCGTCCCGGCCTCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACATCGCCCACGCCGCCGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGCTTCGGGAGCTGATTGAACTGTGCACCAACTGTTTGCggttcg                                                                                    >  2:497903/7‑247 (MQ=255)
                                                                                                                                                          acTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAg                                                                                    <  1:201248/251‑1 (MQ=255)
                                                                                                                                                           cTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGCCATCTCCCCCTTGACGTGCAACCGCGCCCAAGCCGCAGTCATGGTGCCGGTTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAAATGGTTGCGGGt                                                                                      >  1:223557/1‑247 (MQ=255)
                                                                                                                                                               aGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGc                                                                               <  2:39439/251‑1 (MQ=255)
                                                                                                                                                                    ccTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTg                                                                           >  1:271412/1‑250 (MQ=255)
                                                                                                                                                                       cATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACt                                                                       <  1:415533/251‑1 (MQ=255)
                                                                                                                                                                           cTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAAcc                                                                   >  2:341372/1‑251 (MQ=255)
                                                                                                                                                                              ccGATCATATGTGGCAGGAGCGACGCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCgtg                                                                <  1:308584/251‑1 (MQ=255)
                                                                                                                                                                                  ccATATGGGGCAGGAGCGATTCGTCGATCGTCCCGGCATCAAGCAACGCGTCGCCGGTATTGATGAAGACGTCGTGGGCCATCTAACCATTGACGTGCAACAGCGTCTACGCCGCAGTCATGTTGCCGATTTTGGGATCGGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTATCCGTGTTTg                                                            <  2:586466/250‑1 (MQ=255)
                                                                                                                                                                                    atatGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTc                                                          >  2:419895/1‑251 (MQ=255)
                                                                                                                                                                                          ggCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAAc                                                    <  1:279932/251‑1 (MQ=255)
                                                                                                                                                                                           gCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACg                                                   >  2:578971/1‑251 (MQ=255)
                                                                                                                                                                                               gAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGcccc                                               <  1:174970/251‑1 (MQ=255)
                                                                                                                                                                                                aGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCg                                              <  1:560068/251‑1 (MQ=255)
                                                                                                                                                                                                   gACTCGCCGATCGTACCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCATGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGc                                            >  2:208854/1‑250 (MQ=255)
                                                                                                                                                                                                    aCTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCaa                                          <  1:208854/251‑1 (MQ=255)
                                                                                                                                                                                                      ccGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAAtt                                        <  1:370046/250‑1 (MQ=255)
                                                                                                                                                                                                       cGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAAt                                         <  2:401847/249‑1 (MQ=255)
                                                                                                                                                                                                        gCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTaa                                      <  1:513124/251‑1 (MQ=255)
                                                                                                                                                                                                        gCCGAGCGTCCCGGCATCAAGCAACTCCTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGAATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTaa                                      <  2:159568/251‑1 (MQ=255)
                                                                                                                                                                                                               ggcccGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGAGTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCtat                               <  2:122873/249‑1 (MQ=255)
                                                                                                                                                                                                                   cGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATaaaa                           <  1:295382/251‑1 (MQ=255)
                                                                                                                                                                                                                          aaGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGCTTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTACGGTCTATAACAAATACgg                    >  2:565478/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cgtcgCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGCATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGaa             <  2:200587/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cgCCTGTATTGCTGCAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCtata                              >  2:533110/1‑230 (MQ=255)
                                                                                                                                                                                                                                     cgCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCtata                              <  1:533110/230‑1 (MQ=255)
                                                                                                                                                                                                                                      gCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTg         >  2:427730/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAggg    >  1:92358/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAggg    >  1:13217/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGa   <  2:125883/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGAc  <  1:587629/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGAC  >  NC_000913/3581073‑3581559

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