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→GGG) caiCputative 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

CATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  NC_000913/37022‑37477
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
cATTCGGTAACGGTGGCGCGGTACTTCTGGACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCttt                                                                                                                                                                                                                              <  1:930604/239‑1 (MQ=255)
cATTCGGTAACGGGCGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCCCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCttt                                                                                                                                                                                                                              <  1:494336/239‑1 (MQ=255)
      gTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCCCACTGCCAGGCGGAGGAATAT‑CCAGCGAAGCACAGGCTGTAATGGGTAATACCCACA‑CCTTtcggtcg                                                                                                                                                                                                                        >  1:169933/1‑239 (MQ=255)
               gcgcGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCCATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCATACCAGGCGGAGGAATAA‑CAAGCCAAGGGCAGGTTGTACTGGGTAAGCACCACA‑CCTTTCGGTAGGGGggcggg                                                                                                                                                                                                               >  2:27775/1‑238 (MQ=255)
                     tACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAAAGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGTAGGTGGTGCCGGa                                                                                                                                                                                                         >  1:424182/1‑239 (MQ=255)
                         tCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGCCGTCATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATAT‑CCAACAAAGCGGAGGGTGTAATGGGGAATCACCACC‑CCCTTCCGTCGGGGGGGGGGGGCGgggggg                                                                                                                                                                                                     >  2:769431/1‑237 (MQ=255)
                                  gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCCTCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCCGCGAAGCGCAGGTTGTAATGGGTAATCCCCACC‑CCTTTCGGTCGGGGGGGGGGGGCGGGGGGGaaagaaagt                                                                                                                                                                                            >  2:608679/1‑232 (MQ=255)
                                  gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGCAGACGTCATCGTCCCGCAGTGCACACTGCCAGGCGGGGTAATAT‑CCAGCGAAGCCCAGGTTTCAATGGGTAATCACCACA‑CCTTTCGGTCGGGGGGTGGGGCCGGGGGTGAAGGGAAtt                                                                                                                                                                                            >  1:518271/1‑239 (MQ=255)
                                   tCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCGGCGAAGCGCAGGTTGTAATGGGTAATCACCAAA‑CCTTTCGTTCGGGAGGTGGTGCCGGGGGTGAAGGGaagt                                                                                                                                                                                            >  1:625925/1‑236 (MQ=255)
                                             gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGAGCACACTGCCAGGCGGGGTAATAT‑CAAGCGAAGCGCAGGTTGTAATGGGGAATCACCCCA‑CCTTTCGTTCGGGAGGGGGTGCCGGGGGGGAAGAGAATTTCCGCCGTATc                                                                                                                                                                                 >  2:304429/1‑239 (MQ=255)
                                                       tttttctctcagcccaccaaaaGGTG‑GCCCGTGAAAAAACCCCCGCCATCCACCAA‑GTACCCGGGAAAACGATATGAAACGCAGGCATGACCGGCGGGTAGACGTGTGCTTCCAGGAGTGCACTCGTCCGGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCCCCCCC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGAt                                                                                                                                                                       <  2:266498/220‑1 (MQ=255)
                                                          gtttCTTCACCCGCACCAAGGGGTGCCCAGGCACAAAACGCCGCCACGCCCACA‑GTACCCGGTCAATCGATATAAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCCTGCCGGGTGGAGTAATAT‑CCAGCGACGCGCAGGTTGTAATGGGTAATCACGCCC‑CCTTTCGGTCGGGAGGTGTTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                     <  2:616530/238‑1 (MQ=255)
                                                            aaagaaacgggCAAAAGGGGG‑GCCCGGGAAGAAAACCCCGCCATCGCCGGA‑GTACCCTGGCAAACGATAGGAAACGCCGGACTTCCCGTCAGGTAGCCGTAATCGTCACGCAGTGCACCCTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcgg                                                                                                                                                                  <  2:839989/229‑1 (MQ=255)
                                                              cgcgcccAGAACAAGGGTG‑GCCCCGGCAGAAAACGCCGCAATCGCCGCA‑GTACGCGGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                                <  1:1031168/234‑1 (MQ=255)
                                                                 gACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGGGGGGGTGCCGGCGGTGAAGAGAATTTCCGCCGTATctgcagcgggatgcggcgggg                                                                                                                                                             >  2:136348/1‑219 (MQ=255)
                                                                  cccacccacaaGGTG‑GCCCCGGCACAAACCCGCGCTATCGCCCCA‑GTACACGGGCAATCGATATGAAACGCAGGCATTGCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACC‑CCTTTCGTTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGc                                                                                                                                                            <  2:1059517/230‑1 (MQ=255)
                                                                   ccAGCACAAAGGGG‑GCCCCGGCAAAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:387450/239‑1 (MQ=255)
                                                                   ccAGCAAAAAGGTG‑GCCCCGGCAGAAAACCCCGCCATCGCCGCA‑GTACCCTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:955911/239‑1 (MQ=255)
                                                                    aacagcaAAGGTG‑GCCCCGGCAAAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                           <  1:627381/233‑1 (MQ=255)
                                                                       cacaAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCCCACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGGAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGGGGGGAAGAGAATTTCCGCCGTATCGTAAGTCGATAGcggcgggggaaaga                                                                                                                                                       >  2:111712/1‑231 (MQ=255)
                                                                         aaaaGGTG‑GCCCCGGCAAAAAACCCCGTCATCGCCGCATGAACA‑TGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGGAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcac                                                                                                                                                     <  1:303888/238‑1 (MQ=255)
                                                                          aaaGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACCCTGGCAATCGATATGAAACTCAGGCATTGCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcaca                                                                                                                                                    <  1:114309/239‑1 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGcggcgggggctagcctggtccttt                                                                                                                                             >  1:839989/1‑221 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAGGAGAATTTCCCCCGTCTCGGCGGTCGATACCGGCGGTggcatggacaaaggggg                                                                                                                                             >  2:172424/1‑223 (MQ=255)
                                                                                  gCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCCCAGGTTGTAATGGGTAATCACCAAA‑CCTTTCGGTCGGGAGGGGGGGCCGGGGGTGAGGAGAATTTCGCCCGTATCGTCAGGCGATACaggcgggggagagaaaaagagggg                                                                                                                                             >  2:923436/1‑221 (MQ=255)
                                                                                                    gccACCGCCGCA‑GTACACTGGCAATCGATAGGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                                           <  1:179346/239‑1 (MQ=255)
                                                                                                            acaggggtccttGGCAAGTGAAATGAACCGAAGGCATTACCGTCAGGTAGACGTCCTCGTCCCGCAGTGCACACTGCAGGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTGGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATCTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                                                      <  2:399166/226‑1 (MQ=255)
                                                                                                                    ccatGGCAATGGTTTGGACACGCAGGCCTTGCCGTCAGGTAAACTACGTCGTCACGGGGTGCCCCCCGCCAGGCGGAAT‑ATATGCCAGCGAAGCGGAGGTTGTAATGGGTAATCACCCCA‑CCTTTCGGTCGGGAGGTGGTGCGGGAGGTGAAGAGATTTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACt                                                                                                            <  2:350889/236‑1 (MQ=255)
                                                                                                                     acTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATATAACCAGAGGGACGGGTGTTGAttttttacgtagggaaaagaacac                                                                                                           >  2:921033/1‑220 (MQ=255)
                                                                                                                            aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGGATAGCCCAAGGTGGAGGGGTGTTGGTTTTTCCATTTAGTAAAAGAACt                                                                                                            >  1:399166/1‑231 (MQ=255)
                                                                                                                            aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCAAA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGGTAGCGGCGGTGGATAGGACAAGGTGGCCGGGTGTTGATTTTTCAGTTGGGTAAACGACATCACCCca                                                                                                     >  1:350889/1‑238 (MQ=255)
                                                                                                                                                    cggtaGGTAGACGTCAGCTGCACGCAGTGCCCGCGTCCAGGCGGAGTAATCT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTg                                                                          <  1:27775/235‑1 (MQ=255)
                                                                                                                                                       ccggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGCCAGTCGATAGCGCCGGTGCATAGCACAAAGCAGAGGGTTGTTGGTTTTTTCTTTTAGTCAACACAACCACGCCAGCATCAACTGGAGAGTG‑CA‑CATATGTCa                                                                       >  2:438748/3‑239 (MQ=255)
                                                                                                                                                             gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGGCCACGGTGGCAGGTTGTTGATTTTTCAGTTTGGTAAACGAACTCCGGCCCCCATCCGCGGGA‑AGTGCCC‑CATCGTTCA‑GGCaaa                                                                >  1:16213/1‑239 (MQ=255)
                                                                                                                                                                 tCATCGTCACGCAGTGCAAACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACG‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                                                <  1:1046163/190‑1 (MQ=255)
                                                                                                                                                                  cATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAA‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGc                                                           <  1:136348/239‑1 (MQ=255)
                                                                                                                                                                   aTCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGGACAAGGTGGCAGGTTGTTGGTTTTTCAGTTGAGTTAACGAACTAAAGCCCTCATACACGGGA‑AGTGCCC‑CCACTTTCA‑GGGACATTTGc                                                           >  1:1051007/1‑238 (MQ=255)
                                                                                                                                                                    tCGGCAGGCGGTGCCCCCTGCCGGGGGGGTTAATAT‑CCGGGGAGGGGCGGGTTGTAAGGGGTAATCCCCACC‑CCTTTCGGGGGGGGGGGGGTGCCGGGGGTGAAAGAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCg                                                         <  2:1051007/239‑1 (MQ=255)
                                                                                                                                                                    ccGTCAAGCAGTGCACCCTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCA‑CACACCCTTTCGGTGGGGGGGTGGTGCCGGAGGTGAAAAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCg                                                         <  1:261840/238‑1 (MQ=255)
                                                                                                                                                                           gCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGa                                                  <  1:556189/239‑1 (MQ=255)
                                                                                                                                                                             agggCACACGGCCAGGCGGAGTAATAG‑CCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAg                                                 <  1:172424/235‑1 (MQ=255)
                                                                                                                                                                               tGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGGATAGGACAAAGGCGCAGGTTGTGGATTTTTTAGTTGAGTAAAAGAACACCCACCCACATCACAGGGA‑AGTGCCA‑CACCTGTCACCGCAAAtggacgccaatgcaagg                                               >  2:66362/1‑224 (MQ=255)
                                                                                                                                                                                 cccacggacAGGCGGAGTAGTGC‑CCAGGGAAGCGAAGGTTGTAATGGGTAATCCCCACC‑CCTTTCGGTCGGGAGGTGGGGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGAGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTgg                                              <  1:662662/229‑1 (MQ=255)
                                                                                                                                                                                 cacaCTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCCGGTTGTTGATTTTTCCGTTTAGTAAACGAAATCACGCCATCATCAGCGGGA‑AAGTCCA‑CATCTTGCA‑GGCAAATTTGCCGCAATTTAGGGGc                                             >  1:539661/1‑238 (MQ=255)
                                                                                                                                                                                  ccacTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTGCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCa                                            <  1:111712/237‑1 (MQ=255)
                                                                                                                                                                                    acTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGGACTCACGCCATCATCAGCGGGA‑AGGGCCC‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCACCt                                         >  1:891239/1‑239 (MQ=255)
                                                                                                                                                                                       gCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCACAAGGTGGCCGGTTGTTGATTTTTCAGTTGAGTAAACGAACCCACCCCCTCATCAACGGGA‑AGTGCCA‑CAACTGTCA‑GGCAAATGGTCCGCAATTTAGTGGGAtttt                                        >  1:989892/1‑237 (MQ=255)
                                                                                                                                                                                         cAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGGTTTTTCAGTTGAGTAAACGAACTCACGCCATCATCCGCGGGA‑AGGGCCC‑CATCTGTCA‑GGCAAAATGGCCGCAATTGAGTGGCATCTTCct                                     >  1:530944/1‑238 (MQ=255)
                                                                                                                                                                                          cggCGGGGTTATTC‑CCAGGGGAGCCGAGTTTTTAAGGGGAAAACCCCACC‑CCTTTCGGGCGGGGGGTGGTGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCCCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCctgc                                   <  2:625925/238‑1 (MQ=255)
                                                                                                                                                                                              ggAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCCCCACC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa                                <  1:380586/238‑1 (MQ=255)
                                                                                                                                                                                                aGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAACGAACAAGGGGGGAGGTTTTTGATTTTTTAGTTGGGGTAACGGACTCACGCCCACATCAACGGGA‑AGTGGCA‑CATATGTGA‑GGCCAATATGACCCAATTTAGTgggcgcttccactgcaaca                             >  2:165501/1‑222 (MQ=255)
                                                                                                                                                                                                aGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATGACACAAAGTGGCAGGGTTTTGTTTTTTTAGTTTGGGAAACGCAAACACCCCCACCACATCGGGG‑AATGCCACCATCTGTCA‑‑GCAAATATGCCCCCATTAAgggcgagtctccctgtggatc                             >  2:1026795/1‑219 (MQ=255)
                                                                                                                                                                                                  ttttttccatgggAAGCGCGGTTTGTTATGGGTAATCCCCTCC‑CCTTGGGGTCGGGAGGGGGGTCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                             <  2:539661/226‑1 (MQ=255)
                                                                                                                                                                                                   aaTGT‑CCAGGGAAGCGCAGGTTGTAATGGGTAATCCCCTCC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  1:822647/238‑1 (MQ=255)
                                                                                                                                                                                                   aaTAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  1:893829/238‑1 (MQ=255)
                                                                                                                                                                                                    aattgcgggCGAAGGGCAGGTTTTAATGGTTAATCCCCCCC‑CCGTTGGGTCGGGGGGGTGGGCGGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                           <  2:518271/230‑1 (MQ=255)
                                                                                                                                                                                                     tat‑CCAGGAAAGCGCAGGTTTAAATGGGTAATCCCCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  1:1094892/239‑1 (MQ=255)
                                                                                                                                                                                                      tgcccacggaaggccAGGTTGTAAGGGGTAACCCACACT‑CCTTTCGGTCGGGAGGTTGTGCCGGAGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCGCGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  1:902526/225‑1 (MQ=255)
                                                                                                                                                                                                      gccggtgggggtttagaagagTAAGGGTTACTCCCCACC‑CCTTTCGGGGGGGGGGGGGTGGCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                        <  2:181813/219‑1 (MQ=255)
                                                                                                                                                                                                           gacaagcggcgtgTTGTAATGGGTAAGCCCCCAG‑CCTTGCGGGGGGGAGGTGGTGCGGGAGGTGAAGTGAATTTCCGCCGTTTCGTCAGTCGGTAGCGGCGGAGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  2:118473/227‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAt                                                                                                                                                                                  <  1:35411/83‑1 (MQ=255)
                                                                                                                                                                                                           aGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAt                                                                                                                                                                                  >  2:35411/1‑83 (MQ=255)
                                                                                                                                                                                                             cGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGAACAAGGGGGCAGGTTGTTGTTTTTTCAGTTGGGTAAAAGAAATCACGCCACCATCAGCGGGG‑AATGCCC‑CATCTTTAC‑GGCAAATGTGCCGCAATTTAGTGGGCACTTCCTGCTGAAATTGTTTATACATa                 >  1:370105/1‑239 (MQ=255)
                                                                                                                                                                                                              gaatgtgaGGTTGTAATGGGTAATCCCCACA‑CCTTTCGGTCGGGAGGGGGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:864227/229‑1 (MQ=255)
                                                                                                                                                                                                              gAAGGGCGGGTTGTTATGGTTAAGCCCCACC‑CCTTTCGTTCGGGAGGGGTGGCCGGGGGTGAAGAGAATTTCCGCCGTTTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  2:424182/236‑1 (MQ=255)
                                                                                                                                                                                                              gAAGCGTGGGTTGTAATGGGTAATCCCCACC‑CCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATa                 <  1:769431/238‑1 (MQ=255)
                                                                                                                                                                                                              gAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:304429/236‑1 (MQ=255)
                                                                                                                                                                                                              aaaGCGTAGGTTGTAATGGGTAATCACCAAC‑CCGTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                   <  1:954565/235‑1 (MQ=255)
                                                                                                                                                                                                                 gcggAGGTTGGAAGGGGTAATCACCGCA‑CCTTTGGGTCGGGGGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGa              <  1:119182/238‑1 (MQ=255)
                                                                                                                                                                                                                   tcagtgtggaaggggTAATCCCCACC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATa            <  1:401866/226‑1 (MQ=255)
                                                                                                                                                                                                                    cAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAAAGAACAAAGGGGCAGGTTGTTGATTTTTTAGTTGAGTAAAAGAAACCCCGCCAGCATCCACCGGG‑AGTTCCC‑CCTCCGTCA‑GGCCAATGAGCCCCCATTTGTGGGCCGCCTTCTGCTGAAACTGGTAGTAAATAGGAAAGa          >  2:631234/1‑239 (MQ=255)
                                                                                                                                                                                                                            aaTGGGTAATCCCCACA‑CCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCa  <  1:246419/239‑1 (MQ=255)
                                                                                                                                                                                                                            aaTGGGTAAGCCCCACC‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCa  <  1:173809/239‑1 (MQ=255)
                                                                                                                                                                                                                                      caccacA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAg                                                 <  1:42088/182‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
CATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTG‑GCCCCGGCAGAAAACGCCGCCATCGCCGCA‑GTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATAT‑CCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACA‑CCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGA‑AGTGCCA‑CATCTGTCA‑GGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCA  >  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: ‑