Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,811 A→C 5.2% V254G (GTG→GGG)  yihU ← 3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,8110AC 5.2% 333.8 / 4.4 126V254G (GTG→GGG) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (55/61);  new base C (4/4);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATC  >  NC_000913/4072664‑4072957
                                                                                                                                                   |                                                                                                                                                  
ttATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCCGCGGCATTGGcaca                                                                                                                                                  >  1:1631711/1‑150 (MQ=255)
 tctCCAGTCGTTACATTTTTACTGGGGCCGTCATCCCGGCACTGCCACGGACATGTTCCATAATGGCGGACCAGGCCTGGCGACCTCGACCCGCTGCGCGCGCCTGGCTATAACCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAt                                                                                                                                                 <  2:306535/148‑1 (MQ=255)
 tATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCAGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCCGCGGCATTGGcacca                                                                                                                                                 >  1:1395998/1‑148 (MQ=255)
 tATCCAGACGGTACAGATTTACTTTTGCAGGCATCCCGGCACTGACACGGACCTGAGCCAGACTGGCGTACCCGTCCTGGCGACCGCGACCCGCTGCGCTCGCCTGGCTCTAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAt                                                                                                                                                 <  2:306534/150‑1 (MQ=255)
  aTCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATg                                                                                                                                                >  1:580429/1‑150 (MQ=255)
  aTCCAGTCGTTACATTTTTAATTTGGCAGTCAAACCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCCCTGCGCGCGCCGGGCTATAAACCTCCCGTGAGGCGGACCCCAGCGGGATTGGcacctt                                                                                                                                                >  2:1562082/1‑147 (MQ=255)
   tCCAGACGTTACATTTTTACTTTGGCAGTCCGCCCGGCACGGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGc                                                                                                                                               <  2:236070/150‑1 (MQ=255)
    ccAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCCCATGcc                                                                                                                                              >  2:49416/1‑149 (MQ=255)
    ccAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCa                                                                                                                                              >  1:685040/1‑150 (MQ=255)
     cAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGCCCGCGACCCGCTGCGCGCGCCTGGCTATACACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACCTGCAg                                                                                                                                             >  2:358002/1‑150 (MQ=255)
        tCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTg                                                                                                                                          <  2:1262944/150‑1 (MQ=255)
        gcGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGATTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTg                                                                                                                                          <  1:361733/149‑1 (MQ=255)
         cGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCAGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                                                         <  1:1673101/150‑1 (MQ=255)
         cGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGGTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTgg                                                                                                                                         <  1:556319/150‑1 (MQ=255)
         cGTTACATATTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGGCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                                                         <  1:1735339/150‑1 (MQ=255)
          gTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGt                                                                                                                                        >  2:693367/1‑150 (MQ=255)
          gTTACATTTTTACTATGGCAGTCATCCCGGCACTGACACGGACCTGGTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTTCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggg                                                                                                                                        <  1:1494170/150‑2 (MQ=255)
            tACATTTTTACTATGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggggg                                                                                                                                      <  1:220943/150‑4 (MQ=255)
              cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCTGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGc                                                                                                                                    >  2:555273/1‑150 (MQ=255)
               aTTTTTACTTTGGCAGTCCTCCCGCCCCTGCCACGGACCTGATCCAGACTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCg                                                                                                                                   <  2:786161/150‑1 (MQ=255)
                    tACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATc                                                                                                                              >  1:1681449/1‑150 (MQ=255)
                     aCTTTGGCAGTCATCCCGGCACTGACACGGACCTGGTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCa                                                                                                                             <  1:651366/150‑1 (MQ=255)
                      cTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCaa                                                                                                                            >  1:1861227/1‑150 (MQ=255)
                      cTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCaa                                                                                                                            >  2:252240/1‑150 (MQ=255)
                      cTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGACATTGGCACATGCAGCTGGTTGGCGACATCaa                                                                                                                            >  1:480240/1‑150 (MQ=255)
                       tttGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCCGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAg                                                                                                                           <  1:1199010/150‑1 (MQ=255)
                        ttGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCCGCGGCATAGGCACATGCAGCTGGTTGGCGACATCCAgg                                                                                                                          >  1:293018/1‑150 (MQ=255)
                         tGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGGTGGCGACATCAAggg                                                                                                                         >  2:1003446/1‑150 (MQ=255)
                             aGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCCGCGGCATTGGCGCATGCAGCTGGTTGGCGACATCAAGGGCGAt                                                                                                                     >  1:1301552/1‑150 (MQ=255)
                              gTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                    <  2:1560398/150‑1 (MQ=255)
                                cATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGAACATGCAGCTGGTTGGCGACATCAAGGGCGATGcc                                                                                                                  >  1:1523706/1‑150 (MQ=255)
                                 aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCa                                                                                                                 <  1:1690785/150‑1 (MQ=255)
                                 aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCa                                                                                                                 <  2:410680/150‑1 (MQ=255)
                                     cGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGAt                                                                                                             <  2:1799550/150‑1 (MQ=255)
                                        cACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCt                                                                                                          <  2:1512067/150‑1 (MQ=255)
                                        cACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCt                                                                                                          <  2:240385/150‑1 (MQ=255)
                                          cTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGCCCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTa                                                                                                        >  2:1701056/1‑150 (MQ=255)
                                          cTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTa                                                                                                        <  2:1820212/150‑1 (MQ=255)
                                             acacGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATgg                                                                                                     <  2:930647/150‑1 (MQ=255)
                                                cGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCa                                                                                                  >  1:955307/1‑150 (MQ=255)
                                                  gACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGCCCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAg                                                                                                >  2:261540/1‑150 (MQ=255)
                                                   aCCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTGATGGGCAAga                                                                                               >  2:1846846/1‑150 (MQ=255)
                                                     cTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatc                                                                                             >  2:1030449/1‑150 (MQ=255)
                                                      tGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatcg                                                                                            >  1:913722/1‑150 (MQ=255)
                                                      tGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAACATCCTTATGGGCAAgatcg                                                                                            >  1:907684/1‑150 (MQ=255)
                                                          ccAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCa                                                                                        >  2:1781247/1‑150 (MQ=255)
                                                          ccAGAATGGCGGACCACGCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCa                                                                                        <  2:1636876/150‑1 (MQ=255)
                                                           cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCAt                                                                                       <  2:598486/150‑1 (MQ=255)
                                                                tGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCGGGTTGGCGCCATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAgg                                                                                  >  2:1813540/1‑150 (MQ=255)
                                                                   cGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCgg                                                                               >  1:1283606/1‑150 (MQ=255)
                                                                    ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                              >  2:600252/1‑150 (MQ=255)
                                                                    ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                              <  2:1159019/150‑1 (MQ=255)
                                                                          gTCCTGGCGACCGCGCCCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTTTGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                        >  1:424611/1‑150 (MQ=255)
                                                                          gTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                        <  1:579979/150‑1 (MQ=255)
                                                                          gTCCGGGCGACCGCGACCCGCTGCGCGCGCCGGGCTATCAACCTCCCGTGAGGCGGCCCCCAGCGGCAGTGGCCCATGCAGCGGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                        <  1:887114/150‑1 (MQ=255)
                                                                            ccTGGCGCCCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCAGATGCAGCTGGTTGGCGACTTCAAGGGCGATGCCAAGCTCCTGATGGGAAAGCTCGATCATGATAGCGGGAGAAAGAt                                                                      >  2:174160/1‑150 (MQ=255)
                                                                            ccTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCCCATGCAGCTGGTTGGCGACATCAAGGGCTATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAt                                                                      >  2:443212/1‑150 (MQ=255)
                                                                             cTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtc                                                                     <  2:1588358/150‑1 (MQ=255)
                                                                               ggCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGCTCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtctc                                                                   >  2:1384652/1‑150 (MQ=255)
                                                                                    ccgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCTCATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTg                                                              >  2:835685/1‑150 (MQ=255)
                                                                                    ccgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTg                                                              >  1:1778295/1‑150 (MQ=255)
                                                                                     cgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGa                                                             >  2:456940/1‑150 (MQ=255)
                                                                                      gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                            <  2:20939/150‑1 (MQ=255)
                                                                                         accCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGAc                                                         >  1:1537660/1‑150 (MQ=255)
                                                                                           ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                                                       <  1:622646/150‑1 (MQ=255)
                                                                                            cGCTGCGCGCGCCTGGGTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAGCTCCGCTGAGGACttt                                                      <  1:1504803/150‑1 (MQ=255)
                                                                                             gCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttg                                                     <  2:1071508/150‑1 (MQ=255)
                                                                                                gcgcgcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgttt                                                  <  2:1411635/150‑1 (MQ=255)
                                                                                                 cgcgcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtttg                                                 <  2:1740078/150‑1 (MQ=255)
                                                                                                     cgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCa                                             <  1:1646305/150‑1 (MQ=255)
                                                                                                         tGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGaa                                         >  1:156449/1‑150 (MQ=255)
                                                                                                          ggCTATACAACTCCCGTGCGGCGGCCCTCACCGGCAGTGTCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGAGCGAGCATGAAGGCGGGAGACAGATCTCCGCGGAGGACTTTGTTTGGCCAGGAAg                                        <  2:148992/150‑1 (MQ=255)
                                                                                                            cTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtg                                      >  1:1230178/1‑150 (MQ=255)
                                                                                                               taAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCTCTGAGGACTTTGTTTGGCCAGGAAgtggtg                                   <  1:1320838/150‑1 (MQ=255)
                                                                                                                aaaCCCCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCCAGGTCGATCATGAAGGCGGGGGAAAGATCTCCGCTGAGGACTTTGGTTGGCCAGGAAGTGGTGa                                  >  1:55911/1‑150 (MQ=255)
                                                                                                                    cTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGCTCATGAAGGCGGGTGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTg                              >  1:716331/1‑150 (MQ=255)
                                                                                                                       ccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                           <  2:1356581/150‑1 (MQ=255)
                                                                                                                       ccGGGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAATATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                           <  2:1681449/150‑1 (MQ=255)
                                                                                                                          tGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGCTCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTGTGGCCAGGAAGGGGTGAAGTGACCttt                        >  1:922103/1‑150 (MQ=255)
                                                                                                                          tGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCttt                        <  2:982353/150‑1 (MQ=255)
                                                                                                                             ggcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCTGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAcc                     >  2:1261035/1‑150 (MQ=255)
                                                                                                                             ggcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAcc                     >  2:551362/1‑150 (MQ=255)
                                                                                                                             ggcggcCCCCAGCGGCAGTGGCACATGCAGCTGGTTGGCGCCATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGACAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGGGGTGAAGTGCCCTTTAcc                     >  2:1853982/1‑150 (MQ=255)
                                                                                                                              gcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcg                    <  1:1186650/150‑1 (MQ=255)
                                                                                                                              gcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcg                    <  1:1306425/150‑1 (MQ=255)
                                                                                                                                ggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggc                  <  1:1254224/150‑1 (MQ=255)
                                                                                                                                gccaaccAGCTGCATGTGCCAATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCGTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggc                  >  1:1237181/6‑150 (MQ=255)
                                                                                                                                  cccccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGCTCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcgg                >  1:1067403/1‑150 (MQ=255)
                                                                                                                                    cccAGCGGCTTTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATCGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcg              <  1:1665593/150‑1 (MQ=255)
                                                                                                                                    cccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcg              >  2:833358/1‑150 (MQ=255)
                                                                                                                                     ccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGCCCTTTACcggcggcgg             >  1:339309/1‑150 (MQ=255)
                                                                                                                                     ccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcgg             <  2:955307/150‑1 (MQ=255)
                                                                                                                                         cGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAcc         <  2:57109/150‑1 (MQ=255)
                                                                                                                                         cGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAcc         >  1:4171/1‑150 (MQ=255)
                                                                                                                                             aTTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtc     <  2:685040/150‑1 (MQ=255)
                                                                                                                                             aTTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtc     <  1:835540/150‑1 (MQ=255)
                                                                                                                                               tGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCGTTATGGGCAAGATCGATCATGATGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCTCtcat   >  1:227029/1‑150 (MQ=255)
                                                                                                                                               tGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcat   >  1:1394238/1‑150 (MQ=255)
                                                                                                                                                ggCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcatc  <  1:1030449/150‑1 (MQ=255)
                                                                                                                                                ggCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcatc  <  1:1383943/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATC  >  NC_000913/4072664‑4072957

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑