Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,255 A→C 22.4% G190G (GGT→GGG caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,2550AC22.4% 117.5 / 5.3 74G190G (GGT→GGGcaiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base C (2/17);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  NC_000913/37022‑37477
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
cATTCGGTAACGGTGGCGCGGTACTTCTGGACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCttt                                                                                                                                                                                                                              <  1:930604/239‑1 (MQ=255)
cATTCGGTAACGGGCGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCCCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCttt                                                                                                                                                                                                                              <  1:494336/239‑1 (MQ=255)
      gTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCCCACTGCCAGGCGGAGGAATATCCAGCGAAGCACAGGCTGTAATGGGTAATACCCACACCTTtcggtcg                                                                                                                                                                                                                        >  1:169933/1‑239 (MQ=255)
               gcgcGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCCATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCATACCAGGCGGAGGAATAACAAGCCAAGGGCAGGTTGTACTGGGTAAGCACCACACCTTTCGGTAGGGGggcggg                                                                                                                                                                                                               >  2:27775/1‑238 (MQ=255)
                     tACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAAAGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGTAGGTGGTGCCGGa                                                                                                                                                                                                         >  1:424182/1‑239 (MQ=255)
                         tCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGCCGTCATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATATCCAACAAAGCGGAGGGTGTAATGGGGAATCACCACCCCCTTCCGTCGGGGGGGGGGGGCGgggggg                                                                                                                                                                                                     >  2:769431/1‑237 (MQ=255)
                                  gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCCTCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATATCCCGCGAAGCGCAGGTTGTAATGGGTAATCCCCACCCCTTTCGGTCGGGGGGGGGGGGCGGGGGGGaaagaaagt                                                                                                                                                                                            >  2:608679/1‑232 (MQ=255)
                                  gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGCAGACGTCATCGTCCCGCAGTGCACACTGCCAGGCGGGGTAATATCCAGCGAAGCCCAGGTTTCAATGGGTAATCACCACACCTTTCGGTCGGGGGGTGGGGCCGGGGGTGAAGGGAAtt                                                                                                                                                                                            >  1:518271/1‑239 (MQ=255)
                                   tCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCGGCGAAGCGCAGGTTGTAATGGGTAATCACCAAACCTTTCGTTCGGGAGGTGGTGCCGGGGGTGAAGGGaagt                                                                                                                                                                                            >  1:625925/1‑236 (MQ=255)
                                             gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGAGCACACTGCCAGGCGGGGTAATATCAAGCGAAGCGCAGGTTGTAATGGGGAATCACCCCACCTTTCGTTCGGGAGGGGGTGCCGGGGGGGAAGAGAATTTCCGCCGTATc                                                                                                                                                                                 >  2:304429/1‑239 (MQ=255)
                                                       tttttctctcagcccaccaaaaGGTGGCCCGTGAAAAAACCCCCGCCATCCACCAAGTACCCGGGAAAACGATATGAAACGCAGGCATGACCGGCGGGTAGACGTGTGCTTCCAGGAGTGCACTCGTCCGGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCCCCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGAt                                                                                                                                                                       <  2:266498/220‑1 (MQ=255)
                                                          gtttCTTCACCCGCACCAAGGGGTGCCCAGGCACAAAACGCCGCCACGCCCACAGTACCCGGTCAATCGATATAAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCCTGCCGGGTGGAGTAATATCCAGCGACGCGCAGGTTGTAATGGGTAATCACGCCCCCTTTCGGTCGGGAGGTGTTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                     <  2:616530/238‑1 (MQ=255)
                                                            aaagaaacgggCAAAAGGGGGGCCCGGGAAGAAAACCCCGCCATCGCCGGAGTACCCTGGCAAACGATAGGAAACGCCGGACTTCCCGTCAGGTAGCCGTAATCGTCACGCAGTGCACCCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcgg                                                                                                                                                                  <  2:839989/229‑1 (MQ=255)
                                                              cgcgcccAGAACAAGGGTGGCCCCGGCAGAAAACGCCGCAATCGCCGCAGTACGCGGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                                <  1:1031168/234‑1 (MQ=255)
                                                                 gACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGGGGTGCCGGCGGTGAAGAGAATTTCCGCCGTATctgcagcgggatgcggcgggg                                                                                                                                                             >  2:136348/1‑219 (MQ=255)
                                                                  cccacccacaaGGTGGCCCCGGCACAAACCCGCGCTATCGCCCCAGTACACGGGCAATCGATATGAAACGCAGGCATTGCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGTTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGc                                                                                                                                                            <  2:1059517/230‑1 (MQ=255)
                                                                   ccAGCACAAAGGGGGCCCCGGCAAAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:387450/239‑1 (MQ=255)
                                                                   ccAGCAAAAAGGTGGCCCCGGCAGAAAACCCCGCCATCGCCGCAGTACCCTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:955911/239‑1 (MQ=255)
                                                                    aacagcaAAGGTGGCCCCGGCAAAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:627381/233‑1 (MQ=255)
                                                                       cacaAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCCCACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGGAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGGGGGGAAGAGAATTTCCGCCGTATCGTAAGTCGATAGcggcgggggaaaga                                                                                                                                                       >  2:111712/1‑231 (MQ=255)
                                                                         aaaaGGTGGCCCCGGCAAAAAACCCCGTCATCGCCGCATGAACATGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGGAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcac                                                                                                                                                     <  1:303888/238‑1 (MQ=255)
                                                                          aaaGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACCCTGGCAATCGATATGAAACTCAGGCATTGCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcaca                                                                                                                                                    <  1:114309/239‑1 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGcggcgggggctagcctggtccttt                                                                                                                                             >  1:839989/1‑221 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAGGAGAATTTCCCCCGTCTCGGCGGTCGATACCGGCGGTggcatggacaaaggggg                                                                                                                                             >  2:172424/1‑223 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCCCAGGTTGTAATGGGTAATCACCAAACCTTTCGGTCGGGAGGGGGGGCCGGGGGTGAGGAGAATTTCGCCCGTATCGTCAGGCGATACaggcgggggagagaaaaagagggg                                                                                                                                             >  2:923436/1‑221 (MQ=255)
                                                                                                    gccACCGCCGCAGTACACTGGCAATCGATAGGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                                           <  1:179346/239‑1 (MQ=255)
                                                                                                            acaggggtccttGGCAAGTGAAATGAACCGAAGGCATTACCGTCAGGTAGACGTCCTCGTCCCGCAGTGCACACTGCAGGGCGGAGTAATATCCAGCGAAGCGCAGGTGGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATCTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                                                      <  2:399166/226‑1 (MQ=255)
                                                                                                                    ccatGGCAATGGTTTGGACACGCAGGCCTTGCCGTCAGGTAAACTACGTCGTCACGGGGTGCCCCCCGCCAGGCGGAATATATGCCAGCGAAGCGGAGGTTGTAATGGGTAATCACCCCACCTTTCGGTCGGGAGGTGGTGCGGGAGGTGAAGAGATTTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACt                                                                                                            <  2:350889/236‑1 (MQ=255)
                                                                                                                     acTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATATAACCAGAGGGACGGGTGTTGAttttttacgtagggaaaagaacac                                                                                                           >  2:921033/1‑220 (MQ=255)
                                                                                                                            aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGGATAGCCCAAGGTGGAGGGGTGTTGGTTTTTCCATTTAGTAAAAGAACt                                                                                                            >  1:399166/1‑231 (MQ=255)
                                                                                                                            aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCAAACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGGTAGCGGCGGTGGATAGGACAAGGTGGCCGGGTGTTGATTTTTCAGTTGGGTAAACGACATCACCCca                                                                                                     >  1:350889/1‑238 (MQ=255)
                                                                                                                                                    cggtaGGTAGACGTCAGCTGCACGCAGTGCCCGCGTCCAGGCGGAGTAATCTCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTg                                                                          <  1:27775/235‑1 (MQ=255)
                                                                                                                                                       ccggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGCCAGTCGATAGCGCCGGTGCATAGCACAAAGCAGAGGGTTGTTGGTTTTTTCTTTTAGTCAACACAACCACGCCAGCATCAACTGGAGAGTGCACATATGTCa                                                                       >  2:438748/3‑239 (MQ=255)
                                                                                                                                                             gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGGCCACGGTGGCAGGTTGTTGATTTTTCAGTTTGGTAAACGAACTCCGGCCCCCATCCGCGGGAAGTGCCCCATCGTTCAGGCaaa                                                                >  1:16213/1‑239 (MQ=255)
                                                                                                                                                                 tCATCGTCACGCAGTGCAAACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACGCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                                                <  1:1046163/190‑1 (MQ=255)
                                                                                                                                                                  cATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAACCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGc                                                           <  1:136348/239‑1 (MQ=255)
                                                                                                                                                                   aTCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGGACAAGGTGGCAGGTTGTTGGTTTTTCAGTTGAGTTAACGAACTAAAGCCCTCATACACGGGAAGTGCCCCCACTTTCAGGGACATTTGc                                                           >  1:1051007/1‑238 (MQ=255)
                                                                                                                                                                    tCGGCAGGCGGTGCCCCCTGCCGGGGGGGTTAATATCCGGGGAGGGGCGGGTTGTAAGGGGTAATCCCCACCCCTTTCGGGGGGGGGGGGGTGCCGGGGGTGAAAGAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                                         <  2:1051007/239‑1 (MQ=255)
                                                                                                                                                                    ccGTCAAGCAGTGCACCCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACACACCCTTTCGGTGGGGGGGTGGTGCCGGAGGTGAAAAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                                         <  1:261840/238‑1 (MQ=255)
                                                                                                                                                                           gCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGa                                                  <  1:556189/239‑1 (MQ=255)
                                                                                                                                                                             agggCACACGGCCAGGCGGAGTAATAGCCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAg                                                 <  1:172424/235‑1 (MQ=255)
                                                                                                                                                                               tGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGGATAGGACAAAGGCGCAGGTTGTGGATTTTTTAGTTGAGTAAAAGAACACCCACCCACATCACAGGGAAGTGCCACACCTGTCACCGCAAAtggacgccaatgcaagg                                               >  2:66362/1‑224 (MQ=255)
                                                                                                                                                                                 cccacggacAGGCGGAGTAGTGCCCAGGGAAGCGAAGGTTGTAATGGGTAATCCCCACCCCTTTCGGTCGGGAGGTGGGGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGAGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTgg                                              <  1:662662/229‑1 (MQ=255)
                                                                                                                                                                                 cacaCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCCGGTTGTTGATTTTTCCGTTTAGTAAACGAAATCACGCCATCATCAGCGGGAAAGTCCACATCTTGCAGGCAAATTTGCCGCAATTTAGGGGc                                             >  1:539661/1‑238 (MQ=255)
                                                                                                                                                                                  ccacTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTGCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCa                                            <  1:111712/237‑1 (MQ=255)
                                                                                                                                                                                    acTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGGACTCACGCCATCATCAGCGGGAAGGGCCCCATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCACCt                                         >  1:891239/1‑239 (MQ=255)
                                                                                                                                                                                       gCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCACAAGGTGGCCGGTTGTTGATTTTTCAGTTGAGTAAACGAACCCACCCCCTCATCAACGGGAAGTGCCACAACTGTCAGGCAAATGGTCCGCAATTTAGTGGGAtttt                                        >  1:989892/1‑237 (MQ=255)
                                                                                                                                                                                         cAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGGTTTTTCAGTTGAGTAAACGAACTCACGCCATCATCCGCGGGAAGGGCCCCATCTGTCAGGCAAAATGGCCGCAATTGAGTGGCATCTTCct                                     >  1:530944/1‑238 (MQ=255)
                                                                                                                                                                                          cggCGGGGTTATTCCCAGGGGAGCCGAGTTTTTAAGGGGAAAACCCCACCCCTTTCGGGCGGGGGGTGGTGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCCCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgc                                   <  2:625925/238‑1 (MQ=255)
                                                                                                                                                                                              ggAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa                                <  1:380586/238‑1 (MQ=255)
                                                                                                                                                                                                aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAACGAACAAGGGGGGAGGTTTTTGATTTTTTAGTTGGGGTAACGGACTCACGCCCACATCAACGGGAAGTGGCACATATGTGAGGCCAATATGACCCAATTTAGTgggcgcttccactgcaaca                             >  2:165501/1‑222 (MQ=255)
                                                                                                                                                                                                aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATGACACAAAGTGGCAGGGTTTTGTTTTTTTAGTTTGGGAAACGCAAACACCCCCACCACATCGGGGAATGCCACCATCTGTCAGCAAATATGCCCCCATTAAgggcgagtctccctgtggatc                             >  2:1026795/1‑219 (MQ=255)
                                                                                                                                                                                                  ttttttccatgggAAGCGCGGTTTGTTATGGGTAATCCCCTCCCCTTGGGGTCGGGAGGGGGGTCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                             <  2:539661/226‑1 (MQ=255)
                                                                                                                                                                                                   aaTGTCCAGGGAAGCGCAGGTTGTAATGGGTAATCCCCTCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  1:822647/238‑1 (MQ=255)
                                                                                                                                                                                                   aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  1:893829/238‑1 (MQ=255)
                                                                                                                                                                                                    aattgcgggCGAAGGGCAGGTTTTAATGGTTAATCCCCCCCCCGTTGGGTCGGGGGGGTGGGCGGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  2:518271/230‑1 (MQ=255)
                                                                                                                                                                                                     tatCCAGGAAAGCGCAGGTTTAAATGGGTAATCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  1:1094892/239‑1 (MQ=255)
                                                                                                                                                                                                      tgcccacggaaggccAGGTTGTAAGGGGTAACCCACACTCCTTTCGGTCGGGAGGTTGTGCCGGAGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCGCGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  1:902526/225‑1 (MQ=255)
                                                                                                                                                                                                      gccggtgggggtttagaagagTAAGGGTTACTCCCCACCCCTTTCGGGGGGGGGGGGGTGGCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  2:181813/219‑1 (MQ=255)
                                                                                                                                                                                                           gacaagcggcgtgTTGTAATGGGTAAGCCCCCAGCCTTGCGGGGGGGAGGTGGTGCGGGAGGTGAAGTGAATTTCCGCCGTTTCGTCAGTCGGTAGCGGCGGAGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  2:118473/227‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAt                                                                                                                                                                                  <  1:35411/83‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAt                                                                                                                                                                                  >  2:35411/1‑83 (MQ=255)
                                                                                                                                                                                                             cGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGAACAAGGGGGCAGGTTGTTGTTTTTTCAGTTGGGTAAAAGAAATCACGCCACCATCAGCGGGGAATGCCCCATCTTTACGGCAAATGTGCCGCAATTTAGTGGGCACTTCCTGCTGAAATTGTTTATACATa                 >  1:370105/1‑239 (MQ=255)
                                                                                                                                                                                                              gaatgtgaGGTTGTAATGGGTAATCCCCACACCTTTCGGTCGGGAGGGGGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:864227/229‑1 (MQ=255)
                                                                                                                                                                                                              gAAGGGCGGGTTGTTATGGTTAAGCCCCACCCCTTTCGTTCGGGAGGGGTGGCCGGGGGTGAAGAGAATTTCCGCCGTTTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  2:424182/236‑1 (MQ=255)
                                                                                                                                                                                                              gAAGCGTGGGTTGTAATGGGTAATCCCCACCCCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATa                 <  1:769431/238‑1 (MQ=255)
                                                                                                                                                                                                              gAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:304429/236‑1 (MQ=255)
                                                                                                                                                                                                              aaaGCGTAGGTTGTAATGGGTAATCACCAACCCGTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:954565/235‑1 (MQ=255)
                                                                                                                                                                                                                 gcggAGGTTGGAAGGGGTAATCACCGCACCTTTGGGTCGGGGGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGa              <  1:119182/238‑1 (MQ=255)
                                                                                                                                                                                                                   tcagtgtggaaggggTAATCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATa            <  1:401866/226‑1 (MQ=255)
                                                                                                                                                                                                                    cAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAAAGAACAAAGGGGCAGGTTGTTGATTTTTTAGTTGAGTAAAAGAAACCCCGCCAGCATCCACCGGGAGTTCCCCCTCCGTCAGGCCAATGAGCCCCCATTTGTGGGCCGCCTTCTGCTGAAACTGGTAGTAAATAGGAAAGa          >  2:631234/1‑239 (MQ=255)
                                                                                                                                                                                                                            aaTGGGTAATCCCCACACCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCa  <  1:246419/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaTGGGTAAGCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCa  <  1:173809/239‑1 (MQ=255)
                                                                                                                                                                                                                                      caccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAg                                                 <  1:42088/182‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
CATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  NC_000913/37022‑37477

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