Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 39,897 A→G Y164H (TAC→CAC)  caiA ← crotonobetainyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091339,8970AG100.0% 276.0 / NA 78Y164H (TAC→CAC) caiAcrotonobetainyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/38);  total (40/38)

TCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCAC  >  NC_000913/39761‑40030
                                                                                                                                         |                                                                                                                                     
tCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCgg                                                                                                                                    <  2:1486965/140‑1 (MQ=255)
      cagcTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCCGTCGCGCGCCATCACCACGATGTGCGGGgtgt                                                                                                                               <  1:163015/140‑1 (MQ=255)
        gcTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAgg                                                                                                                            <  2:394547/140‑1 (MQ=255)
           aTCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGgcgc                                                                                                                         >  1:1022143/1‑140 (MQ=255)
               aTACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCgctgct                                                                                                                     <  1:195659/140‑1 (MQ=255)
               aTACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCgctgct                                                                                                                     <  2:1022143/140‑1 (MQ=255)
                      gACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaata                                                                                                              >  1:202911/1‑140 (MQ=255)
                      gACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaata                                                                                                              >  2:492445/1‑140 (MQ=255)
                      gACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaata                                                                                                              >  2:392149/1‑140 (MQ=255)
                      gACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaata                                                                                                              >  1:150156/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGCTGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  2:551449/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  2:345166/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  2:1042696/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  2:214375/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  2:1388109/1‑140 (MQ=255)
                       aCCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGtaataa                                                                                                             >  1:842493/1‑140 (MQ=255)
                        ccGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATaaa                                                                                                            <  2:101141/140‑1 (MQ=255)
                        ccGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATaaa                                                                                                            <  1:345166/140‑1 (MQ=255)
                        ccGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATaaa                                                                                                            <  1:492445/140‑1 (MQ=255)
                         cGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATaaaa                                                                                                           <  1:428714/140‑1 (MQ=255)
                          gAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAac                                                                                                          >  1:556597/1‑140 (MQ=255)
                           aGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAaca                                                                                                         <  1:214375/140‑1 (MQ=255)
                           aGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAaca                                                                                                         <  1:283420/140‑1 (MQ=255)
                           aGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAaca                                                                                                         <  1:1042696/140‑1 (MQ=255)
                                 tCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtact                                                                                                   >  2:173660/1‑140 (MQ=255)
                                 tCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtact                                                                                                   <  1:1388109/140‑1 (MQ=255)
                                 tCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtact                                                                                                   <  1:26692/140‑1 (MQ=255)
                                 tCAAGTTTGGGCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtact                                                                                                   >  2:752369/1‑140 (MQ=255)
                                      tttGGTCGCTTTGATGCACGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGGCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                                              <  2:489894/140‑1 (MQ=255)
                                      tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                                              <  2:150156/140‑1 (MQ=255)
                                      tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                                              <  2:556597/140‑1 (MQ=255)
                                      tttGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCAt                                                                                              <  1:1275572/140‑1 (MQ=255)
                                          gTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGTGGTGTAGGCGCTGCTGGTAATaaaa                                                                                                           >  2:1412915/1‑123 (MQ=255)
                                          gTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATaaaa                                                                                                           <  1:1412915/123‑1 (MQ=255)
                                          gTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtac                                                                                                    <  1:485748/130‑1 (MQ=255)
                                          gTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTtac                                                                                                    >  2:485748/1‑130 (MQ=255)
                                          gTCACTTTGATGCCCGGTTTGATCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAg                                                                                          >  2:1078835/1‑140 (MQ=255)
                                                    tGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTa                                                                                <  1:12115/140‑1 (MQ=255)
                                                      cccGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTAcc                                                                              <  2:70980/140‑1 (MQ=255)
                                                         ggTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCAtt                                                                           >  1:522106/1‑140 (MQ=255)
                                                           tttGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTc                                                                         >  2:671919/1‑140 (MQ=255)
                                                               cTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACg                                                                     >  2:225182/1‑140 (MQ=255)
                                                                    atCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                                                >  2:780685/1‑140 (MQ=255)
                                                                    atCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtat                                                                >  2:131530/1‑140 (MQ=255)
                                                                     tCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGtata                                                               >  1:778479/1‑140 (MQ=255)
                                                                         caaACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTc                                                           >  2:955541/1‑140 (MQ=255)
                                                                             ccATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGttt                                                       <  1:392149/140‑1 (MQ=255)
                                                                               aTTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTc                                                     >  1:922902/1‑140 (MQ=255)
                                                                               aTTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTc                                                     <  1:780685/140‑1 (MQ=255)
                                                                                ttCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                                    >  2:832095/1‑140 (MQ=255)
                                                                                ttCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                                    >  2:996291/1‑140 (MQ=255)
                                                                                ttCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                                    >  1:555107/1‑140 (MQ=255)
                                                                                              gTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCg                                      >  1:365525/1‑140 (MQ=255)
                                                                                               tttGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                                    >  2:890840/1‑125 (MQ=255)
                                                                                               tttGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCa                                                    <  1:890840/125‑1 (MQ=255)
                                                                                               tttGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCgg                                     >  1:170042/1‑140 (MQ=255)
                                                                                                         gaAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGcc                           <  1:955541/140‑1 (MQ=255)
                                                                                                         gaAGCCCCGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGcc                           <  2:170042/140‑1 (MQ=255)
                                                                                                              ccccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGt                       >  2:1218403/1‑140 (MQ=255)
                                                                                                              ccccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGt                       <  2:777239/140‑1 (MQ=255)
                                                                                                              ccccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGt                       >  1:633416/1‑140 (MQ=255)
                                                                                                              ccccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGt                       >  2:1023812/1‑140 (MQ=255)
                                                                                                                ccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTc                     <  1:832095/140‑1 (MQ=255)
                                                                                                                ccGTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTc                     <  2:922902/140‑1 (MQ=255)
                                                                                                                  gTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCgg                   >  2:535954/1‑140 (MQ=255)
                                                                                                                  gTCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCgg                   >  2:331944/1‑140 (MQ=255)
                                                                                                                   tCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGt                  >  2:514896/1‑140 (MQ=255)
                                                                                                                   tCGCGCGCCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGt                  >  2:270508/1‑140 (MQ=255)
                                                                                                                        cgcCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCg             <  1:535954/140‑1 (MQ=255)
                                                                                                                         gcCATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGc            <  1:1078835/140‑1 (MQ=255)
                                                                                                                          ccATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCt           <  1:752369/140‑1 (MQ=255)
                                                                                                                           cATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTg          <  1:671919/140‑1 (MQ=255)
                                                                                                                           cATCACCACGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTg          <  2:633416/140‑1 (MQ=255)
                                                                                                                               accacGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGtt      <  1:893376/140‑1 (MQ=255)
                                                                                                                                ccacGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTc     <  1:270508/140‑1 (MQ=255)
                                                                                                                                   cGATGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCcac  >  2:580511/1‑140 (MQ=255)
                                                                                                                                     aTGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCt           <  2:1032226/129‑1 (MQ=255)
                                                                                                                                     aTGTGCGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCt           >  1:1032226/1‑129 (MQ=255)
                                                                                                                                         |                                                                                                                                     
TCACAGCAGCTATCCATACGCAGACCGAGCTTTTCAAGTTTGGTCACTTTGATGCCCGGTTTGCTCATATCAACAAACCATTCGGTGTAGACAGGTTTGTCCGGAGAAGCCCCGTCGCGCGCCATCACCACGATGTACGGGGTGTAGGCGCTGCTGGTAATAAAACACTTACTACCATTAAGATAAATCTTACCATTTCTACGGGTATAAGTCGTTTTCAGGCTACCCACGTCGGAGCCCGCGCCCGGTTCGGTAATCGCTGAGTTCCAC  >  NC_000913/39761‑40030

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