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→GAG) yjjIDUF3029 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

CATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTA  >  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)
                              tGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑Ac                                                                                                                                                                                                           >  1:373890/1‑251 (MQ=255)
                               gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑Ac                                                                                                                                                                                                           <  1:116855/250‑1 (MQ=255)
                               gATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACa                                                                                                                                                                                                          <  1:393804/251‑1 (MQ=255)
                                  cGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGCACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCCCGGGGGGATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGc                                                                                                                                                                                                       >  1:26758/1‑251 (MQ=255)
                                               ccGCCATACGGCAAACGGGGTCCCGGGGTGGTGCCGATATCGGACTTGTTCCCGGACTGGGGGGGTAACTTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGcac                                                                                                                                                                                          <  2:117295/250‑1 (MQ=255)
                                                 gCCATACGGCAAACGCGCGCCCGGCGTGGAGCCGATATCGGAACTGAACCCCGACTGTGCGTGTAACCTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTAGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCCTAGCGCGGGGCAATACATTCTTTCTCAAAAAGCAATTTAACCCCTTCCGCCCTCCCATCACTTGC‑AAACATTGGcaac                                                                                                                                                                                         >  2:591484/1‑248 (MQ=255)
                                                   cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑CCATGCCAAACATTGGCACaaa                                                                                                                                                                                       >  1:393597/1‑250 (MQ=255)
                                                   cATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCCAACATTGGCACaaa                                                                                                                                                                                       >  1:35307/1‑250 (MQ=255)
                                                         gCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACACACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGACAATCCCTTCTTTCTCAAACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACAGTGGCCCAACACgatc                                                                                                                                                                                 >  1:399293/1‑247 (MQ=255)
                                                                   gcCCGGCGGGGTGCCGATATCGGACTTGACCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTTGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCa                                                                                                                                                                      <  2:483957/251‑1 (MQ=255)
                                                                    cccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATc                                                                                                                                                                       >  1:259790/1‑249 (MQ=255)
                                                                     ccGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGc                                                                                                                                                                    <  1:503573/251‑1 (MQ=255)
                                                                        gCGTGGGGGCGATATCGGACCTGAGCCCCGACTGGGCGTGTAACATGGCGAGTGTGTGCCATCCATATTGCGCGGGGGTATTGGCGACAAACGCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCGTCTTTACCGTAGCGCGCGGCAATCCCTTCATTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAG‑ACATGCCAAACATTGGCACAAAGCGTTCAGGGTTAATCAGCCct                                                                                                                                                                 <  2:107184/251‑1 (MQ=255)
                                                                         cGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTAACGGGGGAATTGGCGCCAGACCCCGCCAGTTGCGCGCTGCTGCGATAAACCGCTTCATTTGCGGCGGCTTATTTAACGGAGAGCGCGGCAATCCCTTAATTCGCCCACAGCCAGTTCACCGCTTCCGGCCGCCCAT‑ACATGCCAAACATTGGAACAAAACATTCAAGGTTAATCAAcccct                                                                                                                                                                >  2:520898/1‑249 (MQ=255)
                                                                          ggggtgccgagagtgaaaCTGACCCCCGACTGGGCGTGGAACATGGCGCGTTGTTGCCGGCCATAGTTCACGGGGGTATTGGCGTCAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCAATTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttc                                                                                                                                                               <  2:44725/235‑1 (MQ=255)
                                                                           gggtgCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCcttct                                                                                                                                                              <  1:166918/250‑1 (MQ=255)
                                                                              tgCCGATATCGGAACTGATACCCGACTGTGCGTGTAACATGGCGCGTTTTTACCAGCCATATTTCACGGGGGTATTGGCGACAAACACCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGACTTCTTTACCGAAGCGCGCGGCAAACCCTTATTTATCACACAGAAAGTTAACCGCATCCGCCCGACCAT‑ACCTGCCAAACATTGGCACAAAACTGTCAGGGGTAATCAGcactacttac                                                                                                                                                           >  2:568712/1‑242 (MQ=255)
                                                                               cccGATATCGGAACCGATCCCCGTCTTGGCGTGTAACATGGCGCGTTTTTTCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGTCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCa                                                                                                                                                          <  1:40706/250‑1 (MQ=255)
                                                                                                            cGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGCTAACCCACTTCATTTGCGGCGGCTTCTTTACCGAAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACCTGCCAAACATTGGCACAAAACGTTACGGGTTCATCCGCCCTTCTTTCACCAGGAAGATATTCTCAAAGAAGGGTGAt                                                                                                                             >  2:284419/1‑251 (MQ=255)
                                                                                                             gtgtAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTTGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCATGAAGCTATTCTCAAAGAAGTTTGAtt                                                                                                                            >  2:651683/1‑251 (MQ=255)
                                                                                                                 accagggcccgTTTTTCCCAGCCATATTTCACGGGGGTATTGGCGACAATCTCCGCCAGTTGCGCGCTGATGCGATAACCCACGTCATGTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGAttgttg                                                                                                                        <  2:360620/242‑1 (MQ=255)
                                                                                                                    aTGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATa                                                                                                                     <  1:549201/251‑1 (MQ=255)
                                                                                                                                        atatTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCAt                                                                                                 <  1:137667/251‑1 (MQ=255)
                                                                                                                                          atTTCCCGGGGGTATTGGCGACAAACTCCGCCAGTGGCGCGCGGATGCGATAACCCACTTCATTGGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCg                                                                                               <  2:590987/251‑1 (MQ=255)
                                                                                                                                             tCACGTGGGTATTGGCGCAAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAAGCCATGCTTTCACACACAGCAAGGTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCGTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCgat                                                                                             <  2:232451/250‑1 (MQ=255)
                                                                                                                                                 gggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATgg                                                                                        >  1:301759/1‑251 (MQ=255)
                                                                                                                                                  ggggTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGc                                                                                       <  1:651683/251‑1 (MQ=255)
                                                                                                                                                         tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGatgatg                                                                                         <  1:440842/242‑1 (MQ=255)
                                                                                                                                                         tGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGAATCGatgatg                                                                                         >  2:440842/1‑242 (MQ=255)
                                                                                                                                                             cacaaacatcGCCAGTTGCTCTCTGATGCTTTAACCCCTTTCATTTGCGGCGCCTTCTTTACCGTACATCGCGGCAAACTCTTCTTTTTCTCTCGTCAAGGTTACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGTTTAATCAGCCCTGCTGTTACCAGGAAGCTATTCTTAAAGAAGTGTGATTGTTGCTAGTGGAATTCTCACCGCGCAGCGATGATGGCGAGctgttgct                                                                            <  2:595524/242‑1 (MQ=255)
                                                                                                                                                             atcAAACTCCCCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTGCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATGCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                                                                            <  2:637630/249‑1 (MQ=255)
                                                                                                                                                              aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                           >  1:300515/1‑251 (MQ=255)
                                                                                                                                                              aCAAACTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgctg                                                                           >  1:91505/1‑251 (MQ=255)
                                                                                                                                                                 acctacGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGACGGCTTCTTTCCCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACACTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCCCCAGGAAGCTATTCTCACAGAAGTGTGATTGTTGATAGATGGATTAACACCGCGCATCGATGATGGCAAGCTGCTGCAGGcc                                                                        >  2:526363/6‑250 (MQ=255)
                                                                                                                                                                 aaCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATctgctgct                                                                            >  1:637630/1‑247 (MQ=255)
                                                                                                                                                                  aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGc                                                                         >  1:589550/1‑249 (MQ=255)
                                                                                                                                                                  aCTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCag                                                                       <  1:532155/251‑1 (MQ=255)
                                                                                                                                                                   cTCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGAGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGAACAACACGTTCAAGGTTAATCAACCCTTCTTTCACCCGGAAGCTATTCTCACAGCAGTGTGATTGTTGCTAGGGGAATTCACACCGCGCATCGATGATGGCGTTCGTCTGCTGGccgg                                                                      >  2:352538/1‑248 (MQ=255)
                                                                                                                                                                    tCCGCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCagt                                                                      >  1:577185/1‑250 (MQ=255)
                                                                                                                                                                       gCCAGTTGCGCGCTGATGCGATAACCCACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTg                                                                  <  2:333292/251‑1 (MQ=255)
                                                                                                                                                                                         cgcaaacccccTTCATTTTCGGCGGTTGCTGTACCGTATCGAGCGGCCAGCCCTTCTTTCTTACACAGCTAGGTAACCGCTGTCGCCAGTCCAT‑ACATTCCAAACATTGGCACAAAACGTTGAGGGTGAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCATTAGTGCGGTATAGTGCGCGTaaa                                                <  2:26758/241‑1 (MQ=255)
                                                                                                                                                                                          gATAACCCACTTCATTTGCGGCTGTTTCTTTACCGTAGCGCGGGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATGTGCACAAAACGTTCAGGGTTAAGCAGCCCTGCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGTCAGTAGTGCGGTAGAGTGCGCGTAaag                                               <  2:393597/251‑1 (MQ=255)
                                                                                                                                                                                                 cacgttaTTTGCGGCGGCCTATTTACCGTTGCGCGCGTCAATACCTTCTTCCTCACACAGCAAGTTAACCGCTTCCGCCAGTCCAT‑ACGTGCCAAACAGTGTCACAAAACGGTCAGGGTTAATTAGCCCTTCGGTCTCCAGGAAGCTATTCTCAAAGAAGTTTGATTGGTGATAGAGGAATTCACACCGCGGGTCGTTGTGGGGGATCTGCTGCTTTCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCAt                                        <  1:282048/245‑1 (MQ=255)
                                                                                                                                                                                                     tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCa                                     <  1:224248/250‑1 (MQ=255)
                                                                                                                                                                                                      cATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGc                                   <  2:167390/251‑1 (MQ=255)
                                                                                                                                                                                                       atgtGCGGCGGCTTCTTTACCGGAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCg                                  <  2:494443/248‑1 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGtgtg                                                                                                                               <  1:492482/153‑1 (MQ=255)
                                                                                                                                                                                                            cggcggCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGtgtg                                                                                                                               >  2:492482/1‑153 (MQ=255)
                                                                                                                                                                                                                ggcTTCTTTACCTGAGCCGGCGTCAATCCCTTCTTTCTCTCACGGCAAGTTAGCCGCGTCCGCCAGCCCAT‑ACATGCCAAACATTGTCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATGGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTgc                         <  2:548194/251‑1 (MQ=255)
                                                                                                                                                                                                                ggcTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGCAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCGGTAGTGCGGTAGAGTGCTCGTAAAGAAGTCATCCAGCGATTCGCTgc                         >  2:69787/1‑251 (MQ=255)
                                                                                                                                                                                                                 gcTTCTTTACCGTAGCGCCCGGCAATCCCTTCTTTTTCACACATCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAATGCTCGTAAATAAGTCATCCAGCTATTCGCTg                          <  2:300515/249‑1 (MQ=255)
                                                                                                                                                                                                                       ttACCGTAGCGCGCGGCAATCCCTTCTGTCTCCCTCAGCAAGTTAACCGCGTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGTCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCg                    <  2:586103/249‑1 (MQ=255)
                                                                                                                                                                                                                          ccGTAGCGCGCGGCAATCCCCTCTATACCACACAGCAAGTTAACCGCTGCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCGCAAAGAAGTGTGATTGGTGATAGAGGAAGTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCGGCGCTCGGCAAt               <  2:129894/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATgg             >  1:359829/1‑251 (MQ=255)
                                                                                                                                                                                                                             tAGCGCACAGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAATATTGGCACAAAACGTTCAGGGTTAAACAGCCCTTCTTTAAACAGGAAGCAATTCACAAAGCAGTGCGATTGTTGATAGAGGTATTCACACCGCACCTAGATGATGGCGATCTGCTGCTGGCAGTTGTGCGAGAGAGGACGCGGAAAGAAGTCATCCAGCGGGTCGCGGCGCTTGGCatgggc            >  2:150847/1‑246 (MQ=255)
                                                                                                                                                                                                                                   gcgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCATTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAtt                                                                                                             >  2:273815/1‑148 (MQ=255)
                                                                                                                                                                                                                                   gcgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCATTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAAtt                                                                                                             <  1:273815/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGGGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAgg     >  2:420296/1‑251 (MQ=255)
                                                                                                                                                                                                                                    cgGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAgg     >  1:102972/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaTCCCTTCTTTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTa  >  1:224590/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              
CATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTTTTGCCAGCCATATTTCACGGGGGTATTGGCGACAAACTCCGCCAGTTGCGCGCTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCAT‑ACATGCCAAACATTGGCACAAAACGTTCAGGGTTAATCAGCCCTTCTTTCACCAGGAAGCTATTCTCAAAGAAGTGTGATTGTTGATAGAGGAATTCACACCGCGCATCGATGATGGCGATCTGCTGCTGGCAGTAGTGCGGTAGAGTGCGCGTAAAGAAGTCATCCAGCGATTCGCTGCGCTCGGCAATGGCTTTCAGGTTA  >  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: ‑