Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,611,186 C→A 6.4% Q161K (CAG→AAG)  hyfI → hydrogenase 4 catalytic subunit HyfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,611,1860CA 6.4% 284.9 / 4.7 84Q161K (CAG→AAG) hyfIhydrogenase 4 catalytic subunit HyfI
Reads supporting (aligned to +/- strand):  ref base C (46/32);  new base A (4/2);  total (50/34)
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

GCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCG  >  NC_000913/2611037‑2611334
                                                                                                                                                     |                                                                                                                                                    
gCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAAc                                                                                                                                                      <  2:1495112/150‑1 (MQ=255)
 cGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACa                                                                                                                                                     <  1:1614095/150‑1 (MQ=255)
  ggTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAAAag                                                                                                                                                    >  1:685018/1‑150 (MQ=255)
   gtgtCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACaga                                                                                                                                                   <  1:575339/150‑1 (MQ=255)
    tgtCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACagaa                                                                                                                                                  <  2:124926/150‑1 (MQ=255)
       cggcggcggTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGAt                                                                                                                                               >  2:972924/1‑150 (MQ=255)
        ggcggcggTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGAtt                                                                                                                                              >  2:826652/1‑150 (MQ=255)
            gcggTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACg                                                                                                                                          >  2:705543/1‑150 (MQ=255)
               gTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCtg                                                                                                                                       >  2:1344867/1‑150 (MQ=255)
                   tttCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGa                                                                                                                                   <  1:666786/150‑1 (MQ=255)
                    ttCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGAt                                                                                                                                  >  1:864336/1‑150 (MQ=255)
                       cACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTAt                                                                                                                               >  1:1194933/1‑150 (MQ=255)
                       cACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTAt                                                                                                                               >  2:1018873/1‑150 (MQ=255)
                          gATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATcgc                                                                                                                            <  1:1490041/150‑1 (MQ=255)
                           aTCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATcgcg                                                                                                                           <  2:163588/150‑1 (MQ=255)
                              tctACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATc                                                                                                                        >  1:510916/1‑150 (MQ=255)
                              tctACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATc                                                                                                                        <  2:1452422/150‑1 (MQ=255)
                              tctACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTGTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATc                                                                                                                        <  2:1540870/150‑1 (MQ=255)
                               ctACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATcc                                                                                                                       >  2:1102781/1‑150 (MQ=255)
                                tACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATccc                                                                                                                      >  2:1041414/1‑150 (MQ=255)
                                  cAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCAc                                                                                                                    >  2:1656545/1‑150 (MQ=255)
                                    gCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCg                                                                                                                  >  1:1286520/1‑150 (MQ=255)
                                     cGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCgg                                                                                                                 >  1:243301/1‑150 (MQ=255)
                                       tCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCgggg                                                                                                               <  1:1115790/150‑1 (MQ=255)
                                         tGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTg                                                                                                             <  2:4834/150‑1 (MQ=255)
                                           gggCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCCCCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGAc                                                                                                           >  1:1571881/1‑150 (MQ=255)
                                           gggCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCCCCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGAc                                                                                                           >  2:825387/1‑150 (MQ=255)
                                            ggCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACt                                                                                                          >  2:460079/1‑150 (MQ=255)
                                            ggCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACt                                                                                                          >  2:1276406/1‑150 (MQ=255)
                                            ggCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCCCGGTTTGCTGCATCAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACt                                                                                                          <  2:1248368/150‑1 (MQ=255)
                                               ggTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATg                                                                                                       <  1:705543/150‑1 (MQ=255)
                                                        acCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTtg                                                                                              <  2:1071600/150‑1 (MQ=255)
                                                         ccATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCACCCGTtgt                                                                                             >  1:314708/1‑150 (MQ=255)
                                                         ccATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTtgt                                                                                             >  2:95281/1‑150 (MQ=255)
                                                          cATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTtgtg                                                                                            <  2:1194933/150‑1 (MQ=255)
                                                               tCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGc                                                                                       <  1:918009/150‑1 (MQ=255)
                                                                 cccATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAg                                                                                     >  1:313637/1‑150 (MQ=255)
                                                                  ccATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGa                                                                                    >  2:453903/1‑150 (MQ=255)
                                                                   cATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGAt                                                                                   >  2:1611526/1‑150 (MQ=255)
                                                                     ttGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATcc                                                                                 >  1:505328/1‑150 (MQ=255)
                                                                       gATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCg                                                                               >  2:1064117/1‑150 (MQ=255)
                                                                       gATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCg                                                                               <  1:1344867/150‑1 (MQ=255)
                                                                        aTGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGc                                                                              <  2:1121184/150‑1 (MQ=255)
                                                                        aTGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGc                                                                              <  1:1627724/150‑1 (MQ=255)
                                                                         tGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTCTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGcc                                                                             >  2:247669/1‑150 (MQ=255)
                                                                             tGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATca                                                                         >  1:101048/1‑150 (MQ=255)
                                                                              ggATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATcac                                                                        >  2:664957/1‑150 (MQ=255)
                                                                              ggATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATcac                                                                        >  2:400256/1‑150 (MQ=255)
                                                                                aTCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAg                                                                      <  1:1018873/150‑1 (MQ=255)
                                                                                  ccccGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCg                                                                    <  2:627796/150‑1 (MQ=255)
                                                                                  ccccGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCg                                                                    >  2:1577780/1‑150 (MQ=255)
                                                                                  ccccGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCg                                                                    <  2:220828/150‑1 (MQ=255)
                                                                                   cccGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGt                                                                   >  2:52242/1‑150 (MQ=255)
                                                                                      ggCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATc                                                                <  1:95281/150‑1 (MQ=255)
                                                                                          gCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCa                                                            >  1:301167/1‑150 (MQ=255)
                                                                                           cccGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCAt                                                           >  2:1239293/1‑150 (MQ=255)
                                                                                             cGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGCCTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTg                                                         >  1:1663695/1‑150 (MQ=255)
                                                                                             cGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTg                                                         <  1:400256/150‑1 (MQ=255)
                                                                                                 aaCACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAAAAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCg                                                     <  2:685018/150‑1 (MQ=255)
                                                                                                  acacCGGCCGCCACCATTCACGGTTTCGCCGTGGCGATCGGTTTGCTGCAAAAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCAATGAGCga                                                    >  2:250413/1‑150 (MQ=255)
                                                                                                   cacCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCgag                                                   >  2:736676/1‑150 (MQ=255)
                                                                                                       ggCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAgc                                               <  1:1128296/150‑1 (MQ=255)
                                                                                                            ccaccaTTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCgcggc                                          >  1:1450100/1‑150 (MQ=255)
                                                                                                            ccaccaTTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCgcggc                                          >  1:1444151/1‑150 (MQ=255)
                                                                                                                  ttCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTgg                                    >  1:741991/1‑150 (MQ=255)
                                                                                                                    cACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCg                                  >  1:1570369/1‑150 (MQ=255)
                                                                                                                    cACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCg                                  >  2:1168411/1‑150 (MQ=255)
                                                                                                                    cACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCACCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCg                                  >  2:1663771/1‑150 (MQ=255)
                                                                                                                     aCGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCgg                                 <  1:1276406/150‑1 (MQ=255)
                                                                                                                      cGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCggg                                <  2:510916/150‑1 (MQ=255)
                                                                                                                       ggTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGc                               <  1:247669/150‑1 (MQ=255)
                                                                                                                          ttCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTAt                            >  2:397780/1‑150 (MQ=255)
                                                                                                                           tCGCCGTGGCGCTCGGTTTGCTGCAAAAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATc                           >  2:459466/1‑150 (MQ=255)
                                                                                                                            cGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCg                          <  1:736676/150‑1 (MQ=255)
                                                                                                                                  ggCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGTTCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGGGGGCTATCGTCAggg                    <  2:1464536/150‑1 (MQ=255)
                                                                                                                                    cgcTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGc                  >  2:663695/1‑150 (MQ=255)
                                                                                                                                     gcTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCg                 <  2:864336/150‑1 (MQ=255)
                                                                                                                                       tCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCgag               >  1:329835/1‑150 (MQ=255)
                                                                                                                                           tttGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAAt           <  2:1286520/150‑1 (MQ=255)
                                                                                                                                              gctgcAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTg        >  1:697343/1‑150 (MQ=255)
                                                                                                                                               ctgcAAAAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGc       <  1:459466/150‑1 (MQ=255)
                                                                                                                                                 acAAAAGAAAAATCACGATGTGGATTATCACGATCCAAAAGCGTTTACTAAGCAAACGTAGACGACTCATATAAAACCAACACAGCGTATCTCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGa     >  2:499340/2‑150 (MQ=255)
                                                                                                                                                   aaCAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATc   >  2:1291336/1‑150 (MQ=255)
                                                                                                                                                    aCAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATcc  >  2:1372962/1‑149 (MQ=255)
                                                                                                                                                    aCAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCg  <  1:1291336/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGGTGTCGGCGGCGGTATTTTCCACGATCTCTACAGCGTCTGGGGCGGTAGCGACACCATTGTCCCCATTGATGTTTGGATCCCCGGCTGCCCGCCAACACCGGCCGCCACCATTCACGGTTTCGCCGTGGCGCTCGGTTTGCTGCAACAGAAGATTCACGCTGTGGATTATCGCGATCCCACCGGGGTGACTATGCAACCGTTGTGGCCGCAGATCCCGCCATCACAGCGTATCGCCATTGAGCGAGAAGCGCGGCGGCTGGCGGGCTATCGTCAGGGGCGAGAAATTTGCGATCG  >  NC_000913/2611037‑2611334

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