Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 242,502 C→G 8.6% V268L (GTG→CTG)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913242,5020CG 8.6% 241.7 / 11.3 70V268L (GTG→CTG) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base G (2/4);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGCTCGTC  >  NC_000913/242359‑242649
                                                                                                                                               |                                                                                                                                                   
gagcAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGgtgct                                                                                                                                                >  2:419350/4‑149 (MQ=255)
gATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGgtgct                                                                                                                                                >  2:826276/1‑149 (MQ=255)
 aTCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACctgctg                                                                                                                                               >  3:545127/1‑149 (MQ=255)
     gcgcACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGGTGCTgggga                                                                                                                                           <  4:747972/149‑4 (MQ=255)
     gcgcACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGGTGCTgggga                                                                                                                                           <  1:419350/149‑4 (MQ=255)
     gcgcACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGGTGCTgggga                                                                                                                                           <  4:29817/149‑4 (MQ=255)
     gcgcACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCAGGTGCTgggga                                                                                                                                           <  1:826276/149‑4 (MQ=255)
       gcACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTg                                                                                                                                         >  2:671890/1‑149 (MQ=255)
        cacaGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGc                                                                                                                                        <  3:176172/149‑1 (MQ=255)
           aGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCAt                                                                                                                                     >  4:1344/1‑149 (MQ=255)
                 tGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGccc                                                                                                                               <  1:109197/149‑1 (MQ=255)
                   ccTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCAt                                                                                                                             >  1:433795/1‑149 (MQ=255)
                    cTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATg                                                                                                                            >  1:676116/1‑149 (MQ=255)
                    cTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATg                                                                                                                            >  1:798344/1‑149 (MQ=255)
                    cTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATg                                                                                                                            >  4:719456/1‑149 (MQ=255)
                     tAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGc                                                                                                                           >  3:39719/1‑149 (MQ=255)
                     tAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGc                                                                                                                           >  3:809143/1‑149 (MQ=255)
                     tAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGc                                                                                                                           >  3:94043/1‑149 (MQ=255)
                        aTCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGa                                                                                                                        <  2:840227/149‑1 (MQ=255)
                           ttctGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACaa                                                                                                                     <  1:534933/149‑1 (MQ=255)
                           ttctGCACCGCCGAGTAATTTTGCCGGGTCGGAGAGTTTAATCGCCAGCCCAAGCACGGCCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCTCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACaa                                                                                                                     >  4:550519/1‑149 (MQ=255)
                            tctGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAAt                                                                                                                    >  3:285669/1‑149 (MQ=255)
                                  ccgccgAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGgt                                                                                                              <  3:703770/149‑1 (MQ=255)
                                      cgAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTcc                                                                                                          <  1:658071/149‑1 (MQ=255)
                                          tAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGaa                                                                                                      >  4:807022/1‑149 (MQ=255)
                                          tAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAACCCCGGTGTCCGGaa                                                                                                      <  4:502589/149‑1 (MQ=255)
                                           aaTTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAAt                                                                                                     >  4:691744/1‑149 (MQ=255)
                                                              gTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAAccc                                                                                  <  1:385427/149‑1 (MQ=255)
                                                                ttAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGc                                                                                <  1:353602/149‑1 (MQ=255)
                                                                 tAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCt                                                                               >  2:619840/1‑149 (MQ=255)
                                                                  aaaCGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCACCGGAACCCGCtt                                                                              >  2:527099/1‑149 (MQ=255)
                                                                       ccAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCggg                                                                         <  2:24296/149‑1 (MQ=255)
                                                                            ccAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGt                                                                    <  2:798344/149‑1 (MQ=255)
                                                                            ccAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGt                                                                    >  1:738426/1‑149 (MQ=255)
                                                                                 cACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTg                                                               <  1:584776/149‑1 (MQ=255)
                                                                                     gTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCaaa                                                           <  2:676116/149‑1 (MQ=255)
                                                                                          aaTCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAgg                                                      >  4:591502/1‑149 (MQ=255)
                                                                                              ggTGCCAGCGTTATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCCCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGAt                                                  <  4:545127/149‑1 (MQ=255)
                                                                                              ggTGCCAGCGTGATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGAt                                                  <  1:688525/149‑1 (MQ=255)
                                                                                              ggTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGAt                                                  <  2:801631/149‑1 (MQ=255)
                                                                                              ggTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGACAGTGCAAAGCAGGGGAt                                                  <  4:285669/149‑1 (MQ=255)
                                                                                                 gCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGAtctc                                               <  2:417727/149‑1 (MQ=255)
                                                                                                   cAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCCCCTGCTGGCCCTGCCATCCGCCCATGCAGACAATCCCGGTGTCCGGACTCGCGCCAGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCCAAGCAGGGGATCTccc                                             <  4:491790/149‑2 (MQ=255)
                                                                                                   cAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCt                                             >  4:217054/1‑149 (MQ=255)
                                                                                                     gCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCGGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGa                                           >  4:247487/1‑149 (MQ=255)
                                                                                                        tAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCa                                        >  1:851202/1‑149 (MQ=255)
                                                                                                           tGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAgggg                                                    >  1:214530/1‑134 (MQ=255)
                                                                                                           tGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAgggg                                                    <  2:214530/134‑1 (MQ=255)
                                                                                                                      ttCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCgg                          <  1:421558/149‑1 (MQ=255)
                                                                                                                       tCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGc                         >  3:604182/1‑149 (MQ=255)
                                                                                                                        ccAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCa                        >  4:677574/1‑149 (MQ=255)
                                                                                                                          aGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGa                      >  3:404586/1‑149 (MQ=255)
                                                                                                                           ggTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGAt                     >  2:527026/1‑149 (MQ=255)
                                                                                                                            gTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATa                    >  2:310843/1‑149 (MQ=255)
                                                                                                                            gTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATa                    <  4:272404/149‑1 (MQ=255)
                                                                                                                             tCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAg                   >  2:182587/1‑149 (MQ=255)
                                                                                                                             ccAGAAGCATCTCCAGGACCTACTGGCCCTGCCATTCGCCCATCCAGACAATCCCGGTGTCCGGAATCGCGACCGCATGGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCATGGGATCTCCTGACCTCGCGCCAGACGcagcaggtcg                   >  2:54906/2‑145 (MQ=255)
                                                                                                                               aGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTg                 <  2:61434/149‑1 (MQ=255)
                                                                                                                               aGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTg                 >  2:507216/1‑149 (MQ=255)
                                                                                                                                 aCGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATg                                                                                                                            >  3:11657/1‑40 (MQ=255)
                                                                                                                                 aCGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATg                                                                                                                            <  4:11657/40‑1 (MQ=255)
                                                                                                                                  cGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATc              <  4:177387/149‑1 (MQ=255)
                                                                                                                                      tCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTc          <  2:407904/149‑1 (MQ=255)
                                                                                                                                      tCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTc          <  2:180101/149‑1 (MQ=255)
                                                                                                                                        cccAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTg        <  2:433795/149‑1 (MQ=255)
                                                                                                                                        cccAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTg        >  2:542723/1‑149 (MQ=255)
                                                                                                                                         ccAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGc       <  4:459164/149‑1 (MQ=255)
                                                                                                                                         ccAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGc       >  1:224972/1‑149 (MQ=255)
                                                                                                                                          cAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGCt      <  4:533910/149‑1 (MQ=255)
                                                                                                                                              aCCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGCTCGTc  <  4:39719/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGCATCCCCAGCACCTGCTGGCCCTGCCATTCGCCCATGCAGACAATCCCGGTGTCCGGAATCGCGCCCGCATCGGAACCCGCTTCCGGGCTGGTCAGTGCAAAGCAGGGGATCTCCTGACCACGCGCCAGACGCGGCAGATAGTGATCTTTCTGCTCGTC  >  NC_000913/242359‑242649

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