Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,616,073 T→C E331E (GAA→GAG yjjI ← DUF3029 domain‑containing protein YjjI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,616,0730TC94.0% 196.5 / ‑2.2 69E331E (GAA→GAGyjjIDUF3029 domain‑containing protein YjjI
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (31/32);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTA  >  NC_000913/4615829‑4616310
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
cATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTcacac                                                                                                                                                                                                                                          >  1:280213/1‑251 (MQ=255)
    atgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGc                                                                                                                                                                                                                                       >  1:494443/1‑250 (MQ=255)
      gCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAg                                                                                                                                                                                                                                    <  1:454002/251‑1 (MQ=255)
        ggggCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGtt                                                                                                                                                                                                                                  >  1:633029/1‑251 (MQ=255)
          ggCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTaa                                                                                                                                                                                                                                <  1:155359/251‑1 (MQ=255)
             gACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCg                                                                                                                                                                                                                             >  1:64576/1‑251 (MQ=255)
                 gTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                                                                                                                                                                                                         <  1:544542/251‑1 (MQ=255)
                      cAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCa                                                                                                                                                                                                                    <  1:12685/251‑1 (MQ=255)
                         aTGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGAAATCGGAACTGATCCCCGACTGTGCGTGTAACATGCCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGGTAACCCACTTCATTTGCGGCGGCTTCTTTACCGTATCGCGGGGCAATCCCTTTTTTCTCACACAGAAAGTTAAACGCATCCGCCAgac                                                                                                                                                                                                                 >  2:312075/1‑249 (MQ=255)
                              tGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                                                                                                                                                                                                           >  1:373890/1‑251 (MQ=255)
                               gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                                                                                                                                                                                                           <  1:116855/250‑1 (MQ=255)
                               gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACa                                                                                                                                                                                                          <  1:393804/251‑1 (MQ=255)
                                  cGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGCACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCCCGGGGGGATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGc                                                                                                                                                                                                       >  1:26758/1‑251 (MQ=255)
                                               ccGCCATACGGCAAACGGGGTCCCGGGGTGGTGCCGATATCGGACTTGTTCCCGGACTGGGGGGGTAACTTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGcac                                                                                                                                                                                          <  2:117295/250‑1 (MQ=255)
                                                 gCCATACGGCAAACGCGCGCCCGGCGTGGAGCCGATATCGGAACTGAACCCCGACTGTGCGTGTAACCTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTAGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCCTAGCGCGGGGCAATACATTCTTTCTCAAAAAGCAATTTAACCCCTTCCGCCCTCCCATCACTTGCAAACATTGGcaac                                                                                                                                                                                         >  2:591484/1‑248 (MQ=255)
                                                   cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATCCATGCCAAACATTGGCACaaa                                                                                                                                                                                       >  1:393597/1‑250 (MQ=255)
                                                   cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCCAACATTGGCACaaa                                                                                                                                                                                       >  1:35307/1‑250 (MQ=255)
                                                         gCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACACACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGACAATCCCTTCTTTCTCAAACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAGTGGCCCAACACgatc                                                                                                                                                                                 >  1:399293/1‑247 (MQ=255)
                                                                   gcCCGGCGGGGTGCCGATATCGGACTTGACCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTTGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCa                                                                                                                                                                      <  2:483957/251‑1 (MQ=255)
                                                                    cccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATc                                                                                                                                                                       >  1:259790/1‑249 (MQ=255)
                                                                     ccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGc                                                                                                                                                                    <  1:503573/251‑1 (MQ=255)
                                                                        gCGTGGGGGCGATATCGGACCTGAGCCCCGACTGGGCGTGTAACATGGCGAGTGTGTGCCATCCATATTGCGCGGGGGTATTGGCGACAAACGCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCGTCTTTACCGTAGCGCGCGGCAATCCCTTCATTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAGACATGCCAAACATTGGCACAAAGCGTTCAGGGTTAATCAGCCct                                                                                                                                                                 <  2:107184/251‑1 (MQ=255)
                                                                         cGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTAACGGGGGAATTGGCGCCAGACCCCGCCAGTTGCGCGCTGCTGCGATAAACCGCTTCATTTGCGGCGGCTTATTTAACGGAGAGCGCGGCAATCCCTTAATTCGCCCACAGCCAGTTCACCGCTTCCGGCCGCCCATACATGCCAAACATTGGAACAAAACATTCAAGGTTAATCAAcccct                                                                                                                                                                >  2:520898/1‑249 (MQ=255)
                                                                          ggggtgccgagagtgaaaCTGACCCCCGACTGGGCGTGGAACATGGCGCGTTGTTGCCGGCCATAGTTCACGGGGGTATTGGCGTCAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCAATTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttc                                                                                                                                                               <  2:44725/235‑1 (MQ=255)
                                                                           gggtgCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttct                                                                                                                                                              <  1:166918/250‑1 (MQ=255)
                                                                              tgCCGATATCGGAACTGATACCCGACTGTGCGTGTAACATGGCGCGTTTTTACCAGCCATATTTCACGGGGGTATTGGCGACAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGACTTCTTTACCGAAGCGCGCGGCAAACCCTTATTTATCACACAGAAAGTTAACCGCATCCGCCCGACCATACCTGCCAAACATTGGCACAAAACTGTCAGGGGTAATCAGcactacttac                                                                                                                                                           >  2:568712/1‑242 (MQ=255)
                                                                               cccGATATCGGAACCGATCCCCGTCTTGGCGTGTAACATGGCGCGTTTTTTCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGTCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCa                                                                                                                                                          <  1:40706/250‑1 (MQ=255)
                                                                                                            cGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGCTAACCCACTTCATTTGCGGCGGCTTCTTTACCGAAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACCTGCCAAACATTGGCACAAAACGTTACGGGTTCATCCGCCCTTCTTTCACCAGGAAGATATTCTCAAAGAAGGGTGAt                                                                                                                             >  2:284419/1‑251 (MQ=255)
                                                                                                             gtgtAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTTGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCATGAAGCTATTCTCAAAGAAGTTTGAtt                                                                                                                            >  2:651683/1‑251 (MQ=255)
                                                                                                                 accagggcccgTTTTTCCCAGCCATATTTCACGGGGGTATTGGCGACAATCTCCGCCAGTTGCGCGCTGATGCGATAACCCACGTCATGTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAttgttg                                                                                                                        <  2:360620/242‑1 (MQ=255)
                                                                                                                    aTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATa                                                                                                                     <  1:549201/251‑1 (MQ=255)
                                                                                                                                        atatTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCAt                                                                                                 <  1:137667/251‑1 (MQ=255)
                                                                                                                                          atTTCCCGGGGGTATTGGCGACAAACTCCGCCAGTGGCGCGCGGATGCGATAACCCACTTCATTGGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCg                                                                                               <  2:590987/251‑1 (MQ=255)
                                                                                                                                             tCACGTGGGTATTGGCGCAAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAAGCCATGCTTTCACACACAGCAAGGTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCGTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCgat                                                                                             <  2:232451/250‑1 (MQ=255)
                                                                                                                                                 gggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATgg                                                                                        >  1:301759/1‑251 (MQ=255)
                                                                                                                                                  ggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGc                                                                                       <  1:651683/251‑1 (MQ=255)
                                                                                                                                                         tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGatgatg                                                                                         <  1:440842/242‑1 (MQ=255)
                                                                                                                                                         tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGAATCGatgatg                                                                                         >  2:440842/1‑242 (MQ=255)
                                                                                                                                                             cacaaacatcGCCAGTTGCTCTCTGATGCTTTAACCCCTTTCATTTGCGGCGCCTTCTTTACCGTACATCGCGGCAAACTCTTCTTTTTCTCTCGTCAAGGTTACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGTTTAATCAGCCCTGCTGTTACCAGGAAGCTATTCTTAAAGAAGTGTGATTGTTGCTAGTGGAATTCTCACCGCGCAGCGATGATGGCGAGctgttgct                                                                            <  2:595524/242‑1 (MQ=255)
                                                                                                                                                             atcAAACTCCCCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTGCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATGCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                                                                            <  2:637630/249‑1 (MQ=255)
                                                                                                                                                              aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                           >  1:300515/1‑251 (MQ=255)
                                                                                                                                                              aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                           >  1:91505/1‑251 (MQ=255)
                                                                                                                                                                 acctacGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGACGGCTTCTTTCCCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACACTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCCCCAGGAAGCTATTCTCACAGAAGTGTGATTGTTGATAGATGGATTAACACCGCGCATCGATGATGGCAAGCTGCTGCAGGcc                                                                        >  2:526363/6‑250 (MQ=255)
                                                                                                                                                                 aaCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                                                                            >  1:637630/1‑247 (MQ=255)
                                                                                                                                                                  aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGc                                                                         >  1:589550/1‑249 (MQ=255)
                                                                                                                                                                  aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCag                                                                       <  1:532155/251‑1 (MQ=255)
                                                                                                                                                                   cTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGAGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGAACAACACGTTCAAGGTTAATCAACCCTTCTTTCACCCGGAAGCTATTCTCACAGCAGTGTGATTGTTGCTAGGGGAATTCACACCGCGCATCGATGATGGCGTTCGTCTGCTGGccgg                                                                      >  2:352538/1‑248 (MQ=255)
                                                                                                                                                                    tCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                                                                      >  1:577185/1‑250 (MQ=255)
                                                                                                                                                                       gCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTg                                                                  <  2:333292/251‑1 (MQ=255)
                                                                                                                                                                                         cgcaaacccccTTCATTTTCGGCGGTTGCTGTACCGTATCGAGCGGCCAGCCCTTCTTTCTTACACAGCTAGGTAACCGCTGTCGCCAGTCCATACATTCCAAACATTGGCACAAAACGTTGAGGGTGAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCATTAGTGCGGTATAGTGCGCGTaaa                                                <  2:26758/241‑1 (MQ=255)
                                                                                                                                                                                          gATAACCCACTTCATTTGCGGCTGTTTCTTTACCGTAGCGCGGGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATGTGCACAAAACGTTCAGGGTTAAGCAGCCCTGCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGTCAGTAGTGCGGTAGAGTGCGCGTAaag                                               <  2:393597/251‑1 (MQ=255)
                                                                                                                                                                                                 cacgttaTTTGCGGCGGCCTATTTACCGTTGCGCGCGTCAATACCTTCTTCCTCACACAGCAAGTTAACCGCTTCCGCCAGTCCATACGTGCCAAACAGTGTCACAAAACGGTCAGGGTTAATTAGCCCTTCGGTCTCCAGGAAGCTATTCTCAAAGAAGTTTGATTGGTGATAGAGGAATTCACACCGCGGGTCGTTGTGGGGGATCTGCTGCTTTCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCAt                                        <  1:282048/245‑1 (MQ=255)
                                                                                                                                                                                                     tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCa                                     <  1:224248/250‑1 (MQ=255)
                                                                                                                                                                                                      cATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGc                                   <  2:167390/251‑1 (MQ=255)
                                                                                                                                                                                                       atgtGCGGCGGCTTCTTTACCGGAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCg                                  <  2:494443/248‑1 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGtgtg                                                                                                                               <  1:492482/153‑1 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGtgtg                                                                                                                               >  2:492482/1‑153 (MQ=255)
                                                                                                                                                                                                                ggcTTCTTTACCTGAGCCGGCGTCAATCCCTTCTTTCTCTCACGGCAAGTTAGCCGCGTCCGCCAGCCCATACATGCCAAACATTGTCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATGGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTgc                         <  2:548194/251‑1 (MQ=255)
                                                                                                                                                                                                                ggcTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGCAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCGGTAGTGCGGTAGAGTGCTCGTAAAGAAGTCATCCAGCGATTCGCTgc                         >  2:69787/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcTTCTTTACCGTAGCGCCCGGCAATCCCTTCTTTTTCACACATCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAATGCTCGTAAATAAGTCATCCAGCTATTCGCTg                          <  2:300515/249‑1 (MQ=255)
                                                                                                                                                                                                                       ttACCGTAGCGCGCGGCAATCCCTTCTGTCTCCCTCAGCAAGTTAACCGCGTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGTCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCg                    <  2:586103/249‑1 (MQ=255)
                                                                                                                                                                                                                          ccGTAGCGCGCGGCAATCCCCTCTATACCACACAGCAAGTTAACCGCTGCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCGCAAAGAAGTGTGATTGGTGATAGAGGAAGTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCGGCGCTCGGCAAt               <  2:129894/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATgg             >  1:359829/1‑251 (MQ=255)
                                                                                                                                                                                                                             tAGCGCACAGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAATATTGGCACAAAACGTTCAGGGTTAAACAGCCCTTCTTTAAACAGGAAGCAATTCACAAAGCAGTGCGATTGTTGATAGAGGTATTCACACCGCACCTAGATGATGGCGATCTGCTGCTGGCAGTTGTGCGAGAGAGGACGCGGAAAGAAGTCATCCAGCGGGTCGCGGCGCTTGGCatgggc            >  2:150847/1‑246 (MQ=255)
                                                                                                                                                                                                                                   gcgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCATTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAtt                                                                                                             >  2:273815/1‑148 (MQ=255)
                                                                                                                                                                                                                                   gcgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCATTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAtt                                                                                                             <  1:273815/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGGGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAgg     >  2:420296/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAgg     >  1:102972/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaTCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTa  >  1:224590/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
CATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTA  >  NC_000913/4615829‑4616310

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