Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,055 T→C E21E (GAA→GAG)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,0550TC98.2% 171.0 / ‑6.5 56E21E (GAA→GAG) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/32);  minor base G (0/1);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAA  >  NC_000913/4623812‑4624297
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
cgtggggatggttttatctgtccccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCTCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTCTCTGCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGGTCTGCGCGCGACCAATCACTAACGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACCCTCCTCCCCCg                                                                                                                                                                                                                                              <  2:324751/228‑1 (MQ=255)
           ttttATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAAc                                                                                                                                                                                                                                   >  1:178044/1‑251 (MQ=255)
              tatcgagccacggctgtcacgacaTAAAGGTAAAGAGGGTATTGCCGAGTTGGTGGAACCGCGCGACTCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTACTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGAGACCAACCACTACCGTGTCGGCATTGCGACAGAGAAACA‑GAACGCGCTGCGTGGCGGGAATCAGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCa                                                                                                                                                                                                                                 <  2:207035/228‑1 (MQ=255)
              tATCGTACCCCGCGTTGCCGGCCATAAAGGTAAAGCAGGTATGGCCGATAGTGGGGAACCCCGCTCAGCCACCGCTACCGCGCCGGGCCAGGCGGACGTTATCTTCTTCCACCCGCCGGGTATTACCCTCTTTCCACGGGTGCTGCGCGCGACCAAGCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCa                                                                                                                                                                                                                                 <  2:399333/250‑1 (MQ=255)
                 cggaCTCCGGCATGCCAGCCATAAAAGCAAAGAAGGTATTACCCGGATCGTGGAACAACGCGCCACCACAGCTACTGCGCCGTGCCACGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCAACCAATCACTCCCGTGCAGGCATTGCGCCAGAGAAACA‑GAACGCACTACGTGGCGGTCCTCTGGAGAAAAATACACCACTCCACCGCCAGGTTAACCCACggg                                                                                                                                                                                                                             >  2:25459/4‑251 (MQ=255)
                 cGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTAATGCCGAGATCGTGGAACACCGCGCCGCCACAGCAACTACGCCGGGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTATGCGCGCGACCGATCACTACCGTGTGGGCATTGAGCAAGAGAAACATCAAC‑AGCCGCGGGTCGGGCATTGGGCTAAAAAGACACTCCCCCACCGTCAGGTTAAACAACggg                                                                                                                                                                                                                             >  2:54750/1‑251 (MQ=255)
                  ggaCTCCGGCTTGCCAGACATAAAAGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACC‑GAACGCGCTGCGTGGCTGCCATTGGGCGAAAAATACACTCCTCCACCCCCAGTTTCAACCACGGGt                                                                                                                                                                                                                            >  2:349370/3‑251 (MQ=255)
                  gTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGAGTGGCGGGCATCTGGCGAAAACTACACTCCTCCACCGCCAGGTTAAACCCCgggg                                                                                                                                                                                                                            >  2:48515/1‑250 (MQ=255)
                    aCTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTc                                                                                                                                                                                                                           <  1:295980/250‑1 (MQ=255)
                           cTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGt                                                                                                                                                                                                                   >  1:238658/1‑251 (MQ=255)
                                       acaGGTAACGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCAAGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCTTCTGGCGAAAAATACCCTCCTCCACCGCCAGCTTACAACACGGGTCGTAAGAGTCAGAGATGAGc                                                                                                                                                                                                        >  2:188445/3‑250 (MQ=255)
                                          ggTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGc                                                                                                                                                                                                    <  1:251912/251‑1 (MQ=255)
                                                 cAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAAtgt                                                                                                                                                                                             <  1:107087/251‑1 (MQ=255)
                                                      aTTGCCGAGATGGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCACCCGCCGGGTATTACACTCTTTCCCCGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACa                                                                                                                                                                                        <  2:395830/251‑1 (MQ=255)
                                                             gggtCGTGTAACTCCGCGCCGCCACCGCTACTGCGCCGTTCCAGGGGGACTTTATCTTCGTCCAGCCGCCGGTTATTCCACTCTTTCCACGGGTTCTGCGCGCGACCAATCGCTACCGTGTCGGAATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCGTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCATGCGTAATGTGGACAGAACGAt                                                                                                                                                                                 <  2:469461/248‑1 (MQ=255)
                                                               aTCGTGGAACCCCGCGCCGCCCCCGCTCCTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGAttt                                                                                                                                                                               <  2:466244/251‑1 (MQ=255)
                                                                   tGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCCTTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGGGGACATAACGATTTCCtt                                                                                                                                                                           >  1:17645/1‑251 (MQ=255)
                                                                    ggAACACCGCGCCGCCACCGCTACTGCGCTGTGCCAGGCGGCCGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGGCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCttt                                                                                                                                                                          >  1:404355/1‑251 (MQ=255)
                                                                     gcacccccccccgccTCCGCTACTGCGCCGTGCCAGGCGGGCGTTATCTTCTTCCATCCGCCGGGTATTACCCTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                                                                                                                                                         <  2:360170/241‑1 (MQ=255)
                                                                     gAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTc                                                                                                                                                                         <  1:17232/251‑1 (MQ=255)
                                                                         acCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                                                                                                                                                                     <  1:415409/251‑1 (MQ=255)
                                                                          ccGCGCCGCCTCCGCTTCTGCCCCGTGCCAGGCGGACGTTATCTTCCTCCATCCGCCTGGTATTACTCTCTTTCCTCGGGTTCTGCGCGCTACCAATCACGACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACTCGCGGCGTGGCGGGCATCTGTCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgta                                                                                                                                                                    <  1:227228/251‑1 (MQ=255)
                                                                           cgcgCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                                                                                                                                                                   <  1:266737/251‑1 (MQ=255)
                                                                           cgcgCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtag                                                                                                                                                                   <  1:172547/251‑1 (MQ=255)
                                                                            gcgcCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                                                                                                                                                  <  1:238672/251‑1 (MQ=255)
                                                                             cgccgcCACCGCTACTGCGCCGTGCCAGGCGGACGTTATATTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGGTAAACCACGGGTCGTAAGAGTCAGGGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgtagt                                                                                                                                                                  >  1:175223/1‑250 (MQ=255)
                                                                              gccgccACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTaa                                                                                                                                                                <  1:85598/251‑1 (MQ=255)
                                                                                         gacggcgccgtgccAGGCGGACGTTACCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGCCCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActct                                                                                                                                                     <  2:178044/247‑1 (MQ=255)
                                                                                                  gtgccAGGCGGACGTTATCTGCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCg                                                                                                                                            <  2:78397/249‑1 (MQ=255)
                                                                                                               gTTAGCTTCTGCCATCCGCCGTGTATTACACTCTGTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATTAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTATTAATGGGTCACTCTTTTTTCTCGCTTtcttcctctt                                                                                                                               <  1:221015/251‑1 (MQ=255)
                                                                                                                 tATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCCCGACCACTCACTACCGTCTCGGCATTGCGCCAGAGAAACA‑GCACGCGCTGCGTGGCGGGCATACGGCGAAAAATACACCCCTCCACCGCCCGGATAAACCACGGGTCGTCAAAGTCAGAGATGGGACGGCGTAATGTGGACATAACGGTTTCCTTGCTTGTTGGAATGGGGTACTCTTTTTTCTCGCTTtcttcctcttcc                                                                                                                             >  2:168960/1‑251 (MQ=255)
                                                                                                                  aTCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCg                                                                                                                            >  1:259360/1‑251 (MQ=255)
                                                                                                                              ccgccgGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGc                                                                                                                >  1:312610/1‑251 (MQ=255)
                                                                                                                              cagccgGGTCTTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAAAGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACACCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGAAATAAAGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCAGGCACCGGTTTGc                                                                                                                >  1:312603/3‑251 (MQ=255)
                                                                                                                               caccgGGGATTACACTCTTGCCACGGGTTCTGCGCGAGACCACTCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGATGCGTGGCGGGCCTCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGCAAGAGGCAGAGCTGGGCAGGCGTCATGAGGCCATAACGATTTCCTTTCTTGGAGATATGAGCCCCTCTTTTTTCTCGATTTCTTCCTCTTCACGCACCGGTTTGCt                                                                                                               >  2:315481/3‑251 (MQ=255)
                                                                                                                                  cgGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGc                                                                                                            >  1:273322/1‑251 (MQ=255)
                                                                                                                                      tATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGtt                                                                                                        >  1:112182/1‑251 (MQ=255)
                                                                                                                                         tACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAc                                                                                                     >  2:467263/1‑251 (MQ=255)
                                                                                                                                          acacTCTGTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCTGGCCTCTGGCGAAAAATCCACTCCTCCACCGCCCGGTTAACCCACGGGTCGTAATAGTCAGCGATACGCAGGCGTACTGTGATCATAACGATTTCCTTTCTTGTAGTAATGAGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGTCACCGATTTGCTGGCGGtttta                                                                                                     >  2:383939/1‑247 (MQ=255)
                                                                                                                                                 ttCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGTCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAgg                                                                                             <  1:200945/250‑1 (MQ=255)
                                                                                                                                                   ccACGGGTGCTGCGCGAGTCCAATCTCGACCGTGGCGGCTGTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCGGGCGAAAAATTCACTCCTCCACCGCCAGGTTAAACCACGGTTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCGTGTAGTAATGGGTCACTCGTTTTTCGCGCTTTCTTCCTCTTCCGGCACCGGTTTGCAGGCGGTTAACAGG‑AAGGGc                                                                                           <  2:284822/250‑1 (MQ=255)
                                                                                                                                                      cGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGAtt                                                                                       <  1:196724/251‑1 (MQ=255)
                                                                                                                                                               cgcgcgACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGACAGc                                                                              >  1:336491/1‑251 (MQ=255)
                                                                                                                                                                                       ggCATTGCGCCGGAGAAAAA‑GACCGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACCCCTCCACCGCCAGGTTAAACCACGGGACGGAAGAGCCAGAGCTGAGCAGGCGTCATGTGGCCGCAACGCTTTCCTTTCTTGTAGTAATGGGACCCTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGATTGCTGGCGGTTTACAGGGAAGGGCG‑TTGCTTCCAACggggaccttttccccttccacacg                                                      >  2:327081/1‑230 (MQ=255)
                                                                                                                                                                                         cATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTg                                                    <  1:431199/251‑1 (MQ=255)
                                                                                                                                                                                            ggcgcCAGAGAAACA‑GAACGCGCGGCGTGGCGGGCATCTGGCGAAAAATCCACTCCCCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCg                                                  <  2:220223/249‑1 (MQ=255)
                                                                                                                                                                                                ccAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCa                                             <  1:335622/251‑1 (MQ=255)
                                                                                                                                                                                                 cAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGGTTGCTGGCGGGTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTCCCCTGTAGCAGCGTGCGGGTCAg                                            <  1:51664/251‑1 (MQ=255)
                                                                                                                                                                                                         cA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACATG‑AAGGGCGCTTGCTGCCCACGGGTACGTTTACCCTGTAACAGCGCGCggggcatgactacgc                                    >  2:29766/1‑239 (MQ=255)
                                                                                                                                                                                                          a‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTacgcc                                   <  1:188386/251‑4 (MQ=255)
                                                                                                                                                                                                             aaCGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAgg                                                                                             <  1:105309/191‑1 (MQ=255)
                                                                                                                                                                                                                         gtcGGGCATCTGGCTAAAAATACACCCCTCCACCGCCAGGGTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCTCTTTCTTCCTCTTCCGGCACCGGTGTGCTGGCGGTTAACAGG‑AATGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGTGCAGCACTACGCCGATGGCCAGTgaaa                     <  2:324382/249‑3 (MQ=255)
                                                                                                                                                                                                                                     gCGAAAATTACCCTCCTCCTCCGCCCGGTGAACCCCCGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTTGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTTGCTATCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTAACCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGCAAGCAGCAACATCa         <  2:472950/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAgg       <  1:146428/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTaa  <  1:424207/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAAAAGG‑AAAGGCGGATGCTGCCAGCGGGTAAGTTTACCCTGTAGCAGCGGGCGGGGCAACACTAAGCCGGTTGCCCGTGCAAGAAGCAACATAAAGCGTa   >  1:12078/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
CGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAA  >  NC_000913/4623812‑4624297

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