Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,037 T→C Q27Q (CAA→CAG)  lplA ← lipoate‑‑protein ligase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,0370TC94.4% 150.4 / ‑5.2 55Q27Q (CAA→CAG) lplAlipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (21/30);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGA  >  NC_000913/4623800‑4624278
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
gAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATac                                                                                                                                                                                                                                       <  1:77477/251‑1 (MQ=255)
    acaATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATacaccc                                                                                                                                                                                                                                   >  2:107087/1‑249 (MQ=255)
    acaATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGGTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGCCCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACCctc                                                                                                                                                                                                                                   >  2:266737/1‑251 (MQ=255)
            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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  
GAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACA‑GAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATGAGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGG‑AAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGA  >  NC_000913/4623800‑4624278

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: ‑