Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 36,358 C→A 9.0% A489A (GCG→GCT caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091336,3580CA 9.0% 665.2 / 41.7 198A489A (GCG→GCTcaiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (95/85);  new base A (10/8);  total (105/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

Only 100 of 198 total aligned reads displayed.

GCTTTTGGACGATCAAGGGTAATTTCCAGAATTGATCCATTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGACATTCTCGCCGCCACGTTTAATCATATTGCAGCGGCGATCGACGAAATAAAAAAAGTCCTCTTCGTCGCGGTATCCG  >  NC_000913/36095‑36610
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
gCTTTTGGACGATCAAGGGTAATTTCCAGAATTGATCCATTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCACTTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTg                                                                                                                                                                                                                                                   >  4:289158/1‑275 (MQ=255)
                      tttCCAGAATTGATCCATTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTCTTAGCTGGCGGGGTAGTTCCCTGTTGTTTCATTTATTTCCGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTTATCTCACGATGATAGCGCACTTTCAATCTCGCCTTATTTTGTTCGCAGAACCGGCAAAAtcac                                                                                                                                                                                                                             >  4:357476/1‑272 (MQ=255)
                                  aTCCATTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTGCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                                                                                                                                 >  4:170490/1‑275 (MQ=255)
                                                               cTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAACTGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTTCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATACGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCATCACCACAAAta                                                                                                                                                                                     >  2:274856/1‑273 (MQ=255)
                                                                         ttCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                                                                                                                                 >  4:223079/1‑236 (MQ=255)
                                                                             ccATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGttt                                                                                                                                                                                                             <  2:223079/236‑1 (MQ=255)
                                                                                     gCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGc                                                                                                                                                                                               <  2:205653/242‑1 (MQ=255)
                                                                                       tttttCAGGGACGAAAAAAGCCCTGCAAAAAATGCATATTGTTTTAGAGTTTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCTAGCCATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  4:72680/275‑1 (MQ=255)
                                                                                                   gacAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTGCGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATAcc                                                                                                                                                <  2:170490/275‑1 (MQ=255)
                                                                                                                   aaaaTGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                                                                                                                                 <  2:305918/274‑1 (MQ=255)
                                                                                                                          atatTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCatat                                                                                                                                                                                                                                                       >  6:1013029/1‑149 (MQ=255)
                                                                                                                             ttGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATAtttt                                                                                                                                                                                                                                                    >  7:1033345/1‑149 (MQ=255)
                                                                                                                             ttGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATAtttt                                                                                                                                                                                                                                                    <  5:418541/149‑1 (MQ=255)
                                                                                                                               gTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGt                                                                                                                                                                                                                                                  <  5:144995/149‑1 (MQ=255)
                                                                                                                                ttttAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGtt                                                                                                                                                                                                                                                 >  7:184901/1‑149 (MQ=255)
                                                                                                                                       gtgtGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                                                                                                                                                                                                          <  8:967769/149‑1 (MQ=255)
                                                                                                                                         gtgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGc                                                                                                                                                                                                                                        <  6:720645/149‑1 (MQ=255)
                                                                                                                                          tgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCg                                                                                                                                                                                                                                       >  7:540330/1‑149 (MQ=255)
                                                                                                                                          tgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGa                                                                                                          >  2:103943/1‑274 (MQ=255)
                                                                                                                                           gATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCgg                                                                                                                                                                                                                                      >  6:824427/1‑149 (MQ=255)
                                                                                                                                           gATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCgg                                                                                                                                                                                                                                      >  8:121215/1‑149 (MQ=255)
                                                                                                                                            attattAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGa                                                                                                                                                                                                                                     >  7:298761/1‑149 (MQ=255)
                                                                                                                                              tattaGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGc                                                                                                                                                                                               <  4:683049/185‑1 (MQ=255)
                                                                                                                                                         aGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTc                                                                                                                                                                                                                        >  5:689640/1‑149 (MQ=255)
                                                                                                                                                         aGGGTAGTTCCCTGCTGTTCCATTTATTTCAGATTCCTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTc                                                                                                                                                                                                                        >  7:90783/1‑149 (MQ=255)
                                                                                                                                                            gTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCt                                                                                                                                                                                                                     <  7:121215/149‑1 (MQ=255)
                                                                                                                                                            gTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCt                                                                                                                                                                                                                     >  8:722723/1‑149 (MQ=255)
                                                                                                                                                                ttCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                                                                                                                                 >  8:408551/1‑149 (MQ=255)
                                                                                                                                                                ttCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                                                                                                                                 >  7:581048/1‑149 (MQ=255)
                                                                                                                                                                  cccTGCTGTTTCATTTATTTCATATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGt                                                                                                                                                                                                               >  5:669823/1‑149 (MQ=255)
                                                                                                                                                                    ctgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGttt                                                                                                                                                                                                             >  5:441034/1‑149 (MQ=255)
                                                                                                                                                                     tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTc                                                                                                                                                                                                            <  7:568212/149‑1 (MQ=255)
                                                                                                                                                                     tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGACTCTTTAATACCCACAACCACGATGTCCTGAATTTTCGTGTGCGCCGCGATAATTTTTTCCAGCTCCCCGCAggcaa                                                                               >  4:519101/1‑271 (MQ=255)
                                                                                                                                                                     tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGcgcg                                                                                                              >  4:444044/1‑243 (MQ=255)
                                                                                                                                                                     tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGAc                                                                              >  4:478781/1‑275 (MQ=255)
                                                                                                                                                                     tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGtt                                                                                                                                                                                                                                                 <  8:273970/112‑1 (MQ=255)
                                                                                                                                                                      gctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCa                                                                                                                                                                                                           >  6:607992/1‑149 (MQ=255)
                                                                                                                                                                         gTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  2:688572/193‑1 (MQ=255)
                                                                                                                                                                         gTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGa                                                                                                          <  2:444044/243‑1 (MQ=255)
                                                                                                                                                                          tttCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttcat                                                                                                                                                                                                    >  2:677254/1‑152 (MQ=255)
                                                                                                                                                                            tCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttca                                                                                                                                                                                                     <  8:1033345/149‑1 (MQ=255)
                                                                                                                                                                            tCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttca                                                                                                                                                                                                     >  6:34069/1‑149 (MQ=255)
                                                                                                                                                                                         aTTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGc                                                                                                                                                                                                                                        <  6:425383/101‑1 (MQ=255)
                                                                                                                                                                                          ttctttctAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACaaa                                                                                                                                                                                       >  6:748038/1‑149 (MQ=255)
                                                                                                                                                                                                ctAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt                                                                                                                                                                                 <  7:753995/149‑1 (MQ=255)
                                                                                                                                                                                                 tAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTg                                                                                                                                                                                >  8:853158/1‑149 (MQ=255)
                                                                                                                                                                                                   attattTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGAt                                                                                                                                                                              <  8:532481/149‑1 (MQ=255)
                                                                                                                                                                                                    ttattTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATg                                                                                                                                                                             >  8:288924/1‑149 (MQ=255)
                                                                                                                                                                                                       ttttCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTc                                                                                                                                                                                                                    <  4:2976/107‑1 (MQ=255)
                                                                                                                                                                                                       ttttCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGtt                                                                                                                                                                                                              >  1:50461/1‑113 (MQ=255)
                                                                                                                                                                                                       ttttCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCACTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCt                                                                                                                                                                          <  7:286866/149‑1 (MQ=255)
                                                                                                                                                                                                           ccccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                                                                                                                                                                                                          <  3:50461/113‑1 (MQ=255)
                                                                                                                                                                                                             ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGTAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                                                                                                                                    >  7:734227/1‑149 (MQ=255)
                                                                                                                                                                                                             ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                                                                                                                                    >  7:1012051/1‑149 (MQ=255)
                                                                                                                                                                                                             ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                                                                                                                                    >  5:50302/1‑149 (MQ=255)
                                                                                                                                                                                                              cGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcgc                                                                                                                                                                   >  5:12827/1‑149 (MQ=255)
                                                                                                                                                                                                              cGAGCAATGACGTGACAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCCGAAGCGGAAAAGTTCCTGTTCGCTCAATGTGTCATCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcgc                                                                                                                                                                   <  7:669285/149‑1 (MQ=255)
                                                                                                                                                                                                                aGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGCATTTTCGTGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGACATTCTCGCCGCCACGTTTAATCATATTGCACCGCCGATcgac                                    >  1:137786/1‑274 (MQ=255)
                                                                                                                                                                                                                     ttACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            >  6:114149/1‑149 (MQ=255)
                                                                                                                                                                                                                     ttACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATAc                                                                                                                                                 >  2:291384/1‑160 (MQ=255)
                                                                                                                                                                                                                     ttACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGCCTTCATCGcgaatcga                                                                                                                                                            >  8:465534/1‑149 (MQ=255)
                                                                                                                                                                                                                      tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcgaa                                                                                                                                                           <  6:478722/149‑1 (MQ=255)
                                                                                                                                                                                                                          tGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCtt                                                                                                                                                       <  6:160793/149‑1 (MQ=255)
                                                                                                                                                                                                                             cAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCGCCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                                                                                                                                    <  7:853775/149‑1 (MQ=255)
                                                                                                                                                                                                                             cAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                                                                                                                                    <  8:569999/149‑1 (MQ=255)
                                                                                                                                                                                                                              aGATCTTTTCTGATCTCCAGATAAGAGGGAACTTTGAATTTAGCCATATTTTGTTCGCAGAAGCGTAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATCGCTTCATCGCTAATCGAATCTCTAAt                                                                                                                                                   >  6:495704/1‑149 (MQ=255)
                                                                                                                                                                                                                               gATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                                                                                                                                  <  7:319461/149‑1 (MQ=255)
                                                                                                                                                                                                                               gATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                                                                                                                                  <  6:669823/149‑1 (MQ=255)
                                                                                                                                                                                                                                      tCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACaa                                                                                                                                           >  6:681464/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aTCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTTATGGCTTCATCGcgaa                                                                                                                                                                >  3:104140/1‑124 (MQ=255)
                                                                                                                                                                                                                                           tctcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGACATTCTCGCCGCCACGTTTAATCATATTGCAGCGGCGATCGACGAAATAAAAAAAGTCCTCTTCGTcgcg         <  2:498714/274‑1 (MQ=255)
                                                                                                                                                                                                                                             tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGc                                                                                                                                                                           >  7:489774/1‑110 (MQ=255)
                                                                                                                                                                                                                                             tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGc                                                                                                                                                                           <  8:489774/110‑1 (MQ=255)
                                                                                                                                                                                                                                             tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                                                                                                                                  <  7:641564/135‑1 (MQ=255)
                                                                                                                                                                                                                                              ccAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  1:104140/124‑1 (MQ=255)
                                                                                                                                                                                                                                                aGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                                                                                                                                 <  6:518563/149‑1 (MQ=255)
                                                                                                                                                                                                                                                aGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                                                                                                                                 <  5:748038/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTcc                                                                                                                                >  5:855258/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     agagGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGaa                                                                                                                            >  8:671805/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     agagGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGaa                                                                                                                            >  7:17965/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       agGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCtt                                                                                                                                                                                                                         >  2:502785/1‑54 (MQ=255)
                                                                                                                                                                                                                                                         ggCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                                                                                                                              >  2:629641/1‑47 (MQ=255)
                                                                                                                                                                                                                                                         ggCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaa                                                                                                                                                                >  4:133587/1‑109 (MQ=255)
                                                                                                                                                                                                                                                         ggCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaa                                                                                                                                                                >  2:530977/1‑109 (MQ=255)
                                                                                                                                                                                                                                                          gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCt                                                                                                                               <  4:669456/141‑1 (MQ=255)
                                                                                                                                                                                                                                                          gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTc                                                                                                                       <  6:912229/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCt                                                                                                                                                                                                                     <  4:502785/54‑1 (MQ=255)
                                                                                                                                                                                                                                                           cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  2:217602/111‑1 (MQ=255)
                                                                                                                                                                                                                                                           cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAAc                                                                                                                                          >  2:715275/1‑129 (MQ=255)
                                                                                                                                                                                                                                                            aCTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg                                                                                                                     <  7:430490/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg                                                                                                                     >  6:1064078/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            aCTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg                                                                                                                     <  5:1013029/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg                                                                                                                     <  5:478448/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             cTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  4:530977/109‑1 (MQ=255)
                                                                                                                                                                                                                                                             cTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                                                                                                                            <  2:133587/109‑1 (MQ=255)
                                                                                                                                                                                                                                                             cTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                                                                                                                                 <  2:396238/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                tAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                                                                                                                                    <  4:21295/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaTTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGACATTCTcgccgc                                                                   <  4:694775/193‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaTTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCAGCGCGAAGCGAATCTTTAATACCCACAACCACGATTTCCTGAATTTTCGGGTtcgc                                                                                                               <  8:1061066/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    tttCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCgcgg                                                                                                             >  8:788393/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                            
GCTTTTGGACGATCAAGGGTAATTTCCAGAATTGATCCATTGCGGGTCAGATGTAATGATTCACTCATTCCTTTTCTCCATTTTTGCTTTTTCAGGGACGACAACATCCCTGCAAAAAATGCATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGGCGATAATATTTTCCAGCTCCACGCAGGAGACATTCTCGCCGCCACGTTTAATCATATTGCAGCGGCGATCGACGAAATAAAAAAAGTCCTCTTCGTCGCGGTATCCG  >  NC_000913/36095‑36610

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: