Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 36,358 C→A 12.2% 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,3580CA12.2% 392.7 / 35.1 121A489A (GCG→GCTcaiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (56/50);  new base A (8/7);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 121 total aligned reads displayed.

ATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGG  >  NC_000913/36217‑36503
                                                                                                                                             |                                                                                                                                                 
atatTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCatat                                                                                                                                            >  2:1013029/1‑149 (MQ=255)
   ttGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATAtttt                                                                                                                                         >  3:1033345/1‑149 (MQ=255)
   ttGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATAtttt                                                                                                                                         <  1:418541/149‑1 (MQ=255)
     gTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGt                                                                                                                                       <  1:144995/149‑1 (MQ=255)
     gTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGt                                                                                                                                       <  1:956447/149‑1 (MQ=255)
      ttttAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGtt                                                                                                                                      >  3:184901/1‑149 (MQ=255)
          agagTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCa                                                                                                                                  <  4:52304/149‑1 (MQ=255)
             gtgtGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGaa                                                                                                                               <  4:967769/149‑1 (MQ=255)
               gtgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGc                                                                                                                             <  2:720645/149‑1 (MQ=255)
                tgATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCg                                                                                                                            >  3:540330/1‑149 (MQ=255)
                 gATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCgg                                                                                                                           >  2:824427/1‑149 (MQ=255)
                 gATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCgg                                                                                                                           >  4:121215/1‑149 (MQ=255)
                  attattAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGa                                                                                                                          >  3:298761/1‑149 (MQ=255)
                      ttaGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGaaaaa                                                                                                                      <  3:441192/149‑1 (MQ=255)
                         gCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                   >  4:830457/1‑149 (MQ=255)
                           tGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCct                                                                                                                 <  3:152340/149‑1 (MQ=255)
                             gCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCctct                                                                                                               >  4:568973/1‑149 (MQ=255)
                               aGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTc                                                                                                             >  1:689640/1‑149 (MQ=255)
                               aGGGTAGTTCCCTGCTGTTCCATTTATTTCAGATTCCTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTc                                                                                                             >  3:90783/1‑149 (MQ=255)
                                  gTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCt                                                                                                          >  4:722723/1‑149 (MQ=255)
                                     gTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCaa                                                                                                       <  4:433263/149‑1 (MQ=255)
                                      ttCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                      >  4:408551/1‑149 (MQ=255)
                                      ttCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAAt                                                                                                      >  3:581048/1‑149 (MQ=255)
                                       tCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCCTCGCTCAATg                                                                                                     >  4:469662/1‑149 (MQ=255)
                                        cccTGCTGTTTCATTTATTTCATATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGt                                                                                                    >  1:669823/1‑149 (MQ=255)
                                          ctgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGttt                                                                                                  >  1:441034/1‑149 (MQ=255)
                                           tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTc                                                                                                 <  3:568212/149‑1 (MQ=255)
                                           tgctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGtt                                                                                                                                      >  3:273970/1‑112 (MQ=255)
                                            gctgTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCa                                                                                                >  2:607992/1‑149 (MQ=255)
                                             ctgTTTCATTTATTTCAGATTCTGTCTAATTATTATCCCCGAGCAATTACGTGACAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCTATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCAc                                                                                               <  2:418250/149‑1 (MQ=255)
                                                  tCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttca                                                                                          >  3:1003754/1‑149 (MQ=255)
                                                  tCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttca                                                                                          <  4:1033345/149‑1 (MQ=255)
                                                  tCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttca                                                                                          >  2:34069/1‑149 (MQ=255)
                                                     tttatttCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttcattc                                                                                       >  4:747882/1‑149 (MQ=255)
                                                     tttatttCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttcattc                                                                                       >  3:569999/1‑149 (MQ=255)
                                                     tttatttCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCttcattc                                                                                       <  4:69247/149‑1 (MQ=255)
                                                               aTTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGc                                                                                                                             >  1:425383/1‑101 (MQ=255)
                                                               aTTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGc                                                                                                                             <  2:425383/101‑1 (MQ=255)
                                                               aTTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                   <  2:534891/111‑1 (MQ=255)
                                                               aTTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                   >  1:534891/1‑111 (MQ=255)
                                                                ttctttctAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACaaa                                                                            >  2:748038/1‑149 (MQ=255)
                                                                 tctttctAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAAt                                                                           <  3:968448/149‑1 (MQ=255)
                                                                      ctAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCttt                                                                      <  3:753995/149‑1 (MQ=255)
                                                                       tAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTg                                                                     >  4:853158/1‑149 (MQ=255)
                                                                         attattTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGAt                                                                   <  4:532481/149‑1 (MQ=255)
                                                                          ttattTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATg                                                                  >  4:288924/1‑149 (MQ=255)
                                                                                   ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGTAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                         >  3:734227/1‑149 (MQ=255)
                                                                                   ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                         >  3:1012051/1‑149 (MQ=255)
                                                                                   ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                         >  1:50302/1‑149 (MQ=255)
                                                                                   ccGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcg                                                         >  1:457228/1‑149 (MQ=255)
                                                                                    cGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcgc                                                        >  1:12827/1‑149 (MQ=255)
                                                                                    cGAGCAATGACGTGACAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCCGAAGCGGAAAAGTTCCTGTTCGCTCAATGTGTCATCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATcgc                                                        <  3:669285/149‑1 (MQ=255)
                                                                                           ttACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                 >  3:293909/1‑149 (MQ=255)
                                                                                           ttACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcga                                                 >  2:114149/1‑149 (MQ=255)
                                                                                            tACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaatcgaa                                                <  2:478722/149‑1 (MQ=255)
                                                                                               gTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGca                                                                                   >  3:88842/1‑111 (MQ=255)
                                                                                               gTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGca                                                                                   <  4:88842/111‑1 (MQ=255)
                                                                                                tGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCtt                                            <  2:160793/149‑1 (MQ=255)
                                                                                                tGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCtt                                            <  2:153588/149‑1 (MQ=255)
                                                                                                  gCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTa                                          >  3:617692/1‑149 (MQ=255)
                                                                                                  gCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTa                                          >  3:1018796/1‑149 (MQ=255)
                                                                                                  gCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTa                                          >  3:515512/1‑149 (MQ=255)
                                                                                                   cAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCGCCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                         <  3:853775/149‑1 (MQ=255)
                                                                                                   cAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                         <  2:802469/149‑1 (MQ=255)
                                                                                                   cAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTaa                                         >  3:350888/1‑149 (MQ=255)
                                                                                                    aGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAATCGGAAAAATTCCTCTTCGCTCAATGTTTCa                                                                                                <  2:831871/93‑1 (MQ=255)
                                                                                                    aGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAATCGGAAAAATTCCTCTTCGCTCAATGTTTCa                                                                                                >  1:831871/1‑93 (MQ=255)
                                                                                                    aGATCTTTTCTGATCTCCAGATAAGAGGGAACTTTGAATTTAGCCATATTTTGTTCGCAGAAGCGTAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATCGCTTCATCGCTAATCGAATCTCTAAt                                        >  2:495704/1‑149 (MQ=255)
                                                                                                     gATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                       <  3:319461/149‑1 (MQ=255)
                                                                                                     gATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                       <  2:669823/149‑1 (MQ=255)
                                                                                                         ttttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                   >  4:1078126/1‑69 (MQ=255)
                                                                                                         ttttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTc                                                                                                                   <  3:1078126/69‑1 (MQ=255)
                                                                                                         ttttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaa                                                     >  3:304163/1‑131 (MQ=255)
                                                                                                         ttttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGcgaa                                                     <  4:304163/131‑1 (MQ=255)
                                                                                                         ttttCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCca                                   >  4:284598/1‑149 (MQ=255)
                                                                                                            tCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACaa                                >  2:681464/1‑149 (MQ=255)
                                                                                                                  ctcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGa                          <  2:689640/149‑1 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGc                                                                >  3:489774/1‑110 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGc                                                                <  4:489774/110‑1 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                       <  3:641564/135‑1 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATa                                       >  4:641564/1‑135 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGAt                         >  1:1060731/1‑149 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGAt                         <  2:441034/149‑1 (MQ=255)
                                                                                                                   tcCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGAt                         <  1:133298/149‑1 (MQ=255)
                                                                                                                      aGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                      <  2:518563/149‑1 (MQ=255)
                                                                                                                      aGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTc                      <  1:748038/149‑1 (MQ=255)
                                                                                                                       gATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTcc                     <  3:747882/149‑1 (MQ=255)
                                                                                                                       gATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTcc                     <  2:514430/149‑1 (MQ=255)
                                                                                                                       gATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTcc                     >  1:855258/1‑149 (MQ=255)
                                                                                                                           agagGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGaa                 >  4:671805/1‑149 (MQ=255)
                                                                                                                           agagGGCACTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGaa                 >  3:17965/1‑149 (MQ=255)
                                                                                                                              gggCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAAttt              >  3:712777/1‑149 (MQ=255)
                                                                                                                                gCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTc            <  3:183702/149‑1 (MQ=255)
                                                                                                                                 cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCAc                            <  2:671916/132‑1 (MQ=255)
                                                                                                                                 cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCAc                            >  1:671916/1‑132 (MQ=255)
                                                                                                                                 cACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCg           <  4:219205/149‑1 (MQ=255)
                                                                                                                                  aCTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg          <  1:1013029/149‑1 (MQ=255)
                                                                                                                                  aCTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg          >  2:110328/1‑149 (MQ=255)
                                                                                                                                  aCTTTAAATTTAGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCgg          <  1:478448/149‑1 (MQ=255)
                                                                                                                                        aaTTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCAGCGCGAAGCGAATCTTTAATACCCACAACCACGATTTCCTGAATTTTCGGGTtcgc    <  4:1061066/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
ATATTGTTTTAGAGTGTGATTATTAGCTGGCAGGGTAGTTCCCTGCTGTTTCATTTATTTCAGATTCTTTCTAATTATTTTCCCCGAGCAATTACGTGGCAGATCTTTTCTGATCTCCAGATAAGAGGGCACTTTAAATTTCGCCATATTTTGTTCGCAGAAGCGGAAAAATTCCTCTTCGCTCAATGTTTCACCTTCATTCAGCACCACAAATGCTTTGATGGCTTCATCGCGAATCGAATCTTTAATACCCACAACCACGATGTCCTGAATTTTCGGGTGCGCGG  >  NC_000913/36217‑36503

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