Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,250 A→C 18.6% V192G (GTG→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,2500AC18.6% 132.3 / 2.7 73V192G (GTG→GGG) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (27/26);  new base C (2/18);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  NC_000913/37015‑37477
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
cGGAATCCATTCGGTAACGGTGGCGCGGTACCTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                      <  1:528398/239‑1 (MQ=255)
cGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                      <  1:696273/239‑1 (MQ=255)
  gAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGCCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGTTGGAAAGGGTAACcaccac                                                                                                                                                                                                                                     >  1:476242/1‑238 (MQ=255)
       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)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  NC_000913/37015‑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: