Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 242,344 C→G 7.3% P320P (CCG→CCC)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913242,3440CG 7.3% 573.6 / 28.2 164P320P (CCG→CCC) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (66/86);  new base G (6/6);  total (72/92)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 167 total aligned reads displayed.

GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242196‑242491
                                                                                                                                                     |                                                                                                                                                   
ggCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGccc                                                                                                                                                     >  2:365988/1‑149 (MQ=255)
 gCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCg                                                                                                                                                    >  4:735815/1‑149 (MQ=255)
  cACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCgg                                                                                                                                                   >  4:351232/1‑149 (MQ=255)
   aCTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGc                                                                                                                                                  >  4:1578497/1‑149 (MQ=255)
   aCTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGc                                                                                                                                                  <  4:1014881/149‑1 (MQ=255)
   aCTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGc                                                                                                                                                  >  2:21205/1‑149 (MQ=255)
         ccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtggtg                                                                                                                                            <  2:643022/149‑1 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  3:1131463/149‑4 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  2:1314369/149‑4 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  2:21968/149‑4 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaa                                                                                                                                           <  3:1578497/149‑4 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGtggtgg                                                                                                                                           <  4:273981/149‑1 (MQ=255)
             cATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        <  3:231243/149‑1 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtg                                                                                                                                               <  4:1207176/141‑1 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTgg                                                                                                                                       <  1:1681766/149‑1 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGG‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtg                                                                                                                                               >  3:1207176/1‑141 (MQ=255)
               tCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTggg                                                                                                                                      <  1:318785/149‑1 (MQ=255)
               tCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTggg                                                                                                                                      <  3:735815/149‑1 (MQ=255)
                ccGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     <  4:1671726/149‑1 (MQ=255)
                ccGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTGTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGa                                                                                                                                     <  3:571070/149‑1 (MQ=255)
                 cGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGAt                                                                                                                                    >  3:832245/1‑149 (MQ=255)
                  gCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGCTCGGCACGAAGGCATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATc                                                                                                                                   >  4:1182136/1‑149 (MQ=255)
                  gCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        >  3:934367/1‑144 (MQ=255)
                  gCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTg                                                                                                                                        <  4:934367/144‑1 (MQ=255)
                     aGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAgc                                                                                                                                >  4:810645/1‑149 (MQ=255)
                         ttGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGcac                                                                                                                            >  3:1580022/1‑149 (MQ=255)
                         ttGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGcac                                                                                                                            >  4:1551040/1‑149 (MQ=255)
                          tGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGcaca                                                                                                                           <  3:1077864/149‑1 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                                     <  3:190161/149‑1 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                                     <  1:1660762/149‑1 (MQ=255)
                                 ccATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATg                                                                                                                    <  2:124905/149‑1 (MQ=255)
                                  cATTTTCGGTCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   <  3:314339/149‑1 (MQ=255)
                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   >  3:974645/1‑149 (MQ=255)
                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTGACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                                   <  4:1130770/149‑1 (MQ=255)
                                   cttttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                                <  1:1112768/118‑1 (MQ=255)
                                    ttttCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCt                                                                                                                 <  3:1551040/149‑1 (MQ=255)
                                        cGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAt                                                                                                             <  1:252511/149‑1 (MQ=255)
                                        cGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAt                                                                                                             <  2:377328/149‑1 (MQ=255)
                                         ggCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtc                                                                                                            <  3:679620/149‑1 (MQ=255)
                                          gCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGcaca                                                                                                                           <  1:841988/133‑1 (MQ=255)
                                                 gatgatGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                                      <  3:824375/131‑1 (MQ=255)
                                                 gatgatGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                                      >  4:824375/1‑131 (MQ=255)
                                                 gatgatGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTGCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAc                                                                                                    >  2:278952/1‑149 (MQ=255)
                                                     atgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGCTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgcc                                                                                                <  2:1264048/149‑1 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtctt                                                                                                          <  1:1546280/137‑1 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAAtctt                                                                                                          >  2:1546280/1‑137 (MQ=255)
                                                          atcgatcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTa                                                                                           >  3:688849/1‑149 (MQ=255)
                                                              atcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAttt                                                                                       >  2:29917/1‑149 (MQ=255)
                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  4:132393/149‑1 (MQ=255)
                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTcc                                                                                   <  2:90036/149‑1 (MQ=255)
                                                                  gCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGG‑ACGGTAAGGTTCA‑CGGGGAGTGGCGGCCGCCAATTTCCCAGCCCGGCCTGGGGGGTGGGATCCGCGGACAGGGTATGGCCAAATCTTCTGCCACGCCGAGTAATTTTTTCg                                                                                  >  2:103688/1‑149 (MQ=255)
                                                                         gACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                              >  4:1183681/1‑131 (MQ=255)
                                                                         gACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                              <  3:1183681/131‑1 (MQ=255)
                                                                          aCATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGga                                                                           <  4:1288826/149‑1 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                              >  1:104701/1‑128 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                         <  2:1518218/149‑1 (MQ=255)
                                                                             tCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGagag                                                                        <  1:1458690/149‑1 (MQ=255)
                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       >  1:795654/1‑149 (MQ=255)
                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       <  1:1292779/149‑1 (MQ=255)
                                                                              cTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGt                                                                       >  3:701861/1‑149 (MQ=255)
                                                                               tttACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGtt                                                                      >  4:1399904/1‑149 (MQ=255)
                                                                               tttACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGtt                                                                      <  4:479925/149‑1 (MQ=255)
                                                                               tttACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGCGTCGGAGAGtt                                                                      <  2:58318/149‑1 (MQ=255)
                                                                                   ccGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTaaa                                                                  <  1:21205/149‑1 (MQ=255)
                                                                                    cgcgcgTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAAc                                                                 <  3:1116600/149‑1 (MQ=255)
                                                                                    cgcgcgTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAAc                                                                 <  2:147732/149‑1 (MQ=255)
                                                                                    cgcgcgTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAAc                                                                 <  1:1063574/149‑1 (MQ=255)
                                                                                         gTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAg                                                            >  2:1350107/1‑149 (MQ=255)
                                                                                                 gTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCa                                                    >  2:141992/1‑149 (MQ=255)
                                                                                                      ggAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                               >  1:154486/1‑149 (MQ=255)
                                                                                                        aaCGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGc                                             <  1:1269166/149‑1 (MQ=255)
                                                                                                          cGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCaa                                           <  1:612061/149‑1 (MQ=255)
                                                                                                            gT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATc                                         <  4:509458/149‑1 (MQ=255)
                                                                                                             t‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCg                                        >  2:870431/1‑149 (MQ=255)
                                                                                                                cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                                      <  1:141992/149‑1 (MQ=255)
                                                                                                                cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                                      <  2:105423/149‑1 (MQ=255)
                                                                                                                cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                                      <  4:1117345/149‑1 (MQ=255)
                                                                                                                         ggAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTa                             >  1:285565/1‑149 (MQ=255)
                                                                                                                          gAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCa                                                             >  1:1037498/1‑116 (MQ=255)
                                                                                                                          gAAGTGGCGACGACCAAATTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCa                                                             <  2:1037498/116‑1 (MQ=255)
                                                                                                                           aaGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAAt                           <  2:190400/149‑1 (MQ=255)
                                                                                                                                gCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGc                      <  1:469814/149‑1 (MQ=255)
                                                                                                                                    cgacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAATCGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGttt                  <  2:262965/149‑1 (MQ=255)
                                                                                                                                     gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAg                                                                                                                          >  1:36623/1‑44 (MQ=255)
                                                                                                                                     gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGc                                >  3:894074/1‑134 (MQ=255)
                                                                                                                                     gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGc                                <  4:894074/134‑1 (MQ=255)
                                                                                                                                      acCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                <  3:22784/149‑1 (MQ=255)
                                                                                                                                      acCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                <  4:262152/149‑1 (MQ=255)
                                                                                                                                      acCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt                <  4:688849/149‑1 (MQ=255)
                                                                                                                                        caaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAAt                           <  1:908817/130‑1 (MQ=255)
                                                                                                                                        cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTa                        <  2:766558/139‑1 (MQ=255)
                                                                                                                                        cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc              >  2:160355/1‑149 (MQ=255)
                                                                                                                                           tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  4:701861/149‑1 (MQ=255)
                                                                                                                                            caccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAgg          >  4:1241686/3‑149 (MQ=255)
                                                                                                                                              ccACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           <  3:771481/146‑1 (MQ=255)
                                                                                                                                              ccACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg           >  4:771481/1‑146 (MQ=255)
                                                                                                                                               cACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCa       <  3:1001371/149‑1 (MQ=255)
                                                                                                                                               cACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCa       <  1:1350107/149‑1 (MQ=255)
                                                                                                                                                  gCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGAc    >  2:353236/1‑149 (MQ=255)
                                                                                                                                                  gCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGAc    <  1:456110/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGT‑ACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAGGTCAGACGC  >  NC_000913/242196‑242491

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