Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,421,856 G→C 6.1% G50A (GGG→GCG) ‡ yfcH → putative NAD‑dependent epimerase YfcH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,421,8560GC 6.1% 322.1 / 9.6 98G50A (GGG→GCG) ‡yfcHputative NAD‑dependent epimerase YfcH
Reads supporting (aligned to +/- strand):  ref base G (34/58);  new base C (3/3);  total (37/61)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

TGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGATTTGATT  >  NC_000913/2421709‑2422004
                                                                                                                                                   |                                                                                                                                                    
tGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGc                                                                                                                                                    >  1:1487350/1‑150 (MQ=255)
   aCATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTg                                                                                                                                                 >  2:253639/1‑150 (MQ=255)
    cATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGc                                                                                                                                                <  2:1305180/150‑1 (MQ=255)
    cATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGc                                                                                                                                                <  2:825803/150‑1 (MQ=255)
     aTAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGcc                                                                                                                                               >  1:801789/1‑150 (MQ=255)
      tAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCg                                                                                                                                              >  1:74457/1‑150 (MQ=255)
       aGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGa                                                                                                                                             <  1:859105/150‑1 (MQ=255)
         tGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATc                                                                                                                                           <  2:1445395/150‑1 (MQ=255)
         tGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGCCCTTGccatag                                                                                                                                           <  2:272390/150‑5 (MQ=255)
          gATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCa                                                                                                                                          <  1:863387/150‑1 (MQ=255)
           aTCACCGGAGGGACAGGATTAATTGGTCGCTATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGCCCTTGccatagtg                                                                                                                                         <  2:1462046/150‑7 (MQ=255)
               ccGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCa                                                                                                                                     >  2:480605/1‑150 (MQ=255)
               ccGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCa                                                                                                                                     <  2:522470/150‑1 (MQ=255)
                 ggacggacAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAAc                                                                                                                                   >  2:1531180/1‑150 (MQ=255)
                 ggAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAAc                                                                                                                                   <  1:1181849/150‑1 (MQ=255)
                 ggAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAAc                                                                                                                                   <  1:144803/150‑1 (MQ=255)
                  gAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAAcc                                                                                                                                  <  1:719400/150‑1 (MQ=255)
                  gAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAAcc                                                                                                                                  <  1:851618/150‑1 (MQ=255)
                        cAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACg                                                                                                                            <  2:26703/150‑1 (MQ=255)
                          ggATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGc                                                                                                                          >  2:183840/1‑150 (MQ=255)
                          ggATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGc                                                                                                                          >  1:683685/1‑150 (MQ=255)
                              tAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTg                                                                                                                      <  1:1014521/150‑1 (MQ=255)
                                   ggTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCg                                                                                                                 >  1:95432/1‑150 (MQ=255)
                                    gTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTCACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCgg                                                                                                                >  1:438813/1‑150 (MQ=255)
                                    gTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCgg                                                                                                                >  1:1036580/1‑150 (MQ=255)
                                    gTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGCCGTTGaggttg                                                                                                                <  1:946608/150‑6 (MQ=255)
                                     tCGGCATTTGATTCCGCGTTTGGTGGAGCTGGGGGATCAAATTACGGTAGTGACGGGTAACCCGCAGAAAGCCAGTACCGTTCTGGGCCGTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCTACGGGGTTGATGCGGt                                                                                                               <  2:897701/150‑1 (MQ=255)
                                     tCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGt                                                                                                               <  1:1531180/150‑1 (MQ=255)
                                     tCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGCCGTTGaggttgc                                                                                                               <  1:1038378/150‑7 (MQ=255)
                                       gCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTaa                                                                                                             >  1:515611/1‑150 (MQ=255)
                                         cATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATc                                                                                                           <  1:702209/150‑1 (MQ=255)
                                             tGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAAcc                                                                                                       >  1:985750/1‑150 (MQ=255)
                                               aTTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTg                                                                                                     >  2:133570/1‑150 (MQ=255)
                                                ttCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTgg                                                                                                    >  1:1363277/1‑150 (MQ=255)
                                                 tCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGc                                                                                                   <  1:508692/150‑1 (MQ=255)
                                                   cgcgTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCg                                                                                                 <  1:1507923/150‑1 (MQ=255)
                                                    gcgTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCgg                                                                                                <  2:552239/150‑1 (MQ=255)
                                                    gcgTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCgg                                                                                                <  2:712602/150‑1 (MQ=255)
                                                       tttGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGgaga                                                                                             <  2:1011150/150‑1 (MQ=255)
                                                         tGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAAc                                                                                           >  1:20530/1‑150 (MQ=255)
                                                          gCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAAcc                                                                                          <  2:1363277/150‑1 (MQ=255)
                                                             ggAGCTGGGCCATCAAATTACGGTATTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACAGTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCTGAGAACCGAt                                                                                       <  2:1215712/150‑1 (MQ=255)
                                                              gAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGAtt                                                                                      <  1:134011/150‑1 (MQ=255)
                                                              gAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGAtt                                                                                      <  1:1061204/150‑1 (MQ=255)
                                                                 cTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCt                                                                                   <  1:64188/150‑1 (MQ=255)
                                                                  tGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCTTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGAAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTg                                                                                  <  2:148063/150‑1 (MQ=255)
                                                                  tGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTg                                                                                  <  2:1237408/150‑1 (MQ=255)
                                                                    ggCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGAt                                                                                <  2:312189/150‑1 (MQ=255)
                                                                     gCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATa                                                                               <  2:43521/150‑1 (MQ=255)
                                                                     gCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATa                                                                               <  1:800463/150‑1 (MQ=255)
                                                                       cATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATaaa                                                                             <  1:991798/150‑1 (MQ=255)
                                                                          cAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGc                                                                          <  1:239412/150‑1 (MQ=255)
                                                                           aaaTTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCt                                                                         <  1:1054349/150‑1 (MQ=255)
                                                                             aTTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTgg                                                                       >  2:945547/1‑150 (MQ=255)
                                                                               tACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGAc                                                                     <  2:140405/150‑1 (MQ=255)
                                                                               tACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGCCCTTGccatag                                                                                                                                           <  2:1122279/80‑5 (MQ=255)
                                                                               tACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGCCCTTGccatag                                                                                                                                           >  1:1122279/1‑76 (MQ=255)
                                                                                aCGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACt                                                                    <  1:253639/150‑1 (MQ=255)
                                                                                aCGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACt                                                                    <  1:727422/150‑1 (MQ=255)
                                                                                  ggTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCa                                                                  <  2:1452334/150‑1 (MQ=255)
                                                                                           gcgTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCaaaa                                                         >  2:386310/1‑150 (MQ=255)
                                                                                               aaCCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAagag                                                     <  2:381114/150‑1 (MQ=255)
                                                                                                aCCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCACCCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGc                                                    >  1:983465/1‑150 (MQ=255)
                                                                                                aCCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGc                                                    >  1:742863/1‑150 (MQ=255)
                                                                                                 cccGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCg                                                   >  1:120510/1‑150 (MQ=255)
                                                                                                 cccGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCg                                                   >  1:537517/1‑150 (MQ=255)
                                                                                                  ccGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTgg                                                                       <  2:871379/129‑1 (MQ=255)
                                                                                                  ccGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTgg                                                                       >  1:871379/1‑129 (MQ=255)
                                                                                                  ccGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGt                                                  >  1:993258/1‑150 (MQ=255)
                                                                                                     cAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGtctc                                               >  1:326301/1‑150 (MQ=255)
                                                                                                       gAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTg                                             <  2:499580/150‑1 (MQ=255)
                                                                                                       gAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTg                                             <  2:494448/150‑1 (MQ=255)
                                                                                                       gAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTg                                             >  1:452894/1‑150 (MQ=255)
                                                                                                       gAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTg                                             <  1:133570/150‑1 (MQ=255)
                                                                                                       gAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTg                                             <  1:945547/150‑1 (MQ=255)
                                                                                                        aaaGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGc                                            <  2:985750/150‑1 (MQ=255)
                                                                                                               gTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGcc                                     >  1:593322/1‑150 (MQ=255)
                                                                                                                   cGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTg                                 <  1:1095466/150‑1 (MQ=255)
                                                                                                                    gTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTgg                                >  1:1163970/1‑150 (MQ=255)
                                                                                                                    gTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTgg                                <  2:438813/150‑1 (MQ=255)
                                                                                                                       ctcGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAat                             <  1:1129099/150‑1 (MQ=255)
                                                                                                                         cGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACTCTGGACTCACTAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTTGAatat                           <  2:557188/150‑1 (MQ=255)
                                                                                                                          ggCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATc                          <  2:1096188/150‑1 (MQ=255)
                                                                                                                          ggCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATc                          >  2:1208914/1‑150 (MQ=255)
                                                                                                                          ggCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATc                          >  2:506097/1‑150 (MQ=255)
                                                                                                                              cTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGc                      <  1:480605/150‑1 (MQ=255)
                                                                                                                                    tGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAAc                <  1:386310/150‑1 (MQ=255)
                                                                                                                                       cacTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTgg             >  1:1357224/1‑150 (MQ=255)
                                                                                                                                        acTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGt            <  2:1487350/150‑1 (MQ=255)
                                                                                                                                         gatcggCAAGCCCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTc           >  2:855247/5‑150 (MQ=255)
                                                                                                                                         gatcggCAAGCCCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTc           >  2:1117955/5‑150 (MQ=255)
                                                                                                                                         cTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTc           >  2:1370254/1‑150 (MQ=255)
                                                                                                                                             ggCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGAtt       <  1:183840/150‑1 (MQ=255)
                                                                                                                                              gCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGtc                                                 <  1:970428/107‑1 (MQ=255)
                                                                                                                                              gCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGtc                                                 >  2:970428/1‑107 (MQ=255)
                                                                                                                                              gCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGAttt      <  2:326301/150‑1 (MQ=255)
                                                                                                                                               cAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGATTTg     >  2:1200592/1‑150 (MQ=255)
                                                                                                                                                aaGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGATTTGa    <  2:120510/150‑1 (MQ=255)
                                                                                                                                                  gggCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGATTTGAtt  >  2:299178/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGAACATAGTGATCACCGGAGGGACAGGATTAATTGGTCGCCATTTGATTCCGCGTTTGCTGGAGCTGGGCCATCAAATTACGGTAGTGACGCGTAACCCGCAGAAAGCCAGTTCCGTTCTCGGCCCTCGGGTGACACTATGGCAAGGGCTTGCCGATCAAAGCAACCTCAACGGCGTTGATGCGGTAATCAACCTGGCCGGAGAACCGATTGCTGATAAACGCTGGACTCACGAGCAAAAAGAGCGTCTCTGCCAAAGCCGCTGGAATATCACGCAAAAACTGGTCGATTTGATT  >  NC_000913/2421709‑2422004

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: