Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,454,591 G→A 10.6% intergenic (+13/+20) yjfF → / ← fbp putative sugar ABC transporter permease/fructose‑1,6‑bisphosphatase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,454,5910GA10.6% 254.2 / 16.4 85intergenic (+13/+20)yjfF/fbpputative sugar ABC transporter permease/fructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base G (35/41);  new base A (2/7);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCT  >  NC_000913/4454444‑4454739
                                                                                                                                                   |                                                                                                                                                    
cTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgt                                                                                                                                                     <  3:118654/149‑1 (MQ=255)
 tGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtg                                                                                                                                                    >  2:396646/1‑149 (MQ=255)
 tGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtg                                                                                                                                                    >  1:61738/1‑149 (MQ=255)
  gAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtgt                                                                                                                                                   <  2:405681/149‑1 (MQ=255)
  gAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCgtgt                                                                                                                                                   <  4:448584/149‑1 (MQ=255)
   aGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGtgtg                                                                                                                                                  <  3:146238/149‑1 (MQ=255)
        cTGGTGGACGAAAATCGCCATCGGTGTTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                             <  3:577070/149‑1 (MQ=255)
        cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                             <  1:167838/149‑1 (MQ=255)
        cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                             >  2:468983/1‑149 (MQ=255)
        cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                             <  3:227581/149‑1 (MQ=255)
        cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                             >  3:264484/1‑149 (MQ=255)
          ggtggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGgc                                                                                                                                           <  3:8859/149‑1 (MQ=255)
          ggtggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgc                                                                                                                                           <  2:6000/149‑1 (MQ=255)
           gtggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGTCTCTCCACACGGAGAGgcg                                                                                                                                          <  3:695241/149‑1 (MQ=255)
            tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                         <  4:635871/149‑2 (MQ=255)
            tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                         <  1:178731/149‑2 (MQ=255)
            tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                         <  2:23000/149‑2 (MQ=255)
              gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                       <  4:256891/149‑1 (MQ=255)
              gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                       <  2:46700/149‑1 (MQ=255)
              gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                       <  3:173654/149‑1 (MQ=255)
                 gAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGgag                                                                                                                                    >  2:453175/1‑149 (MQ=255)
                      tCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTAc                                                                                                                               >  3:785202/1‑149 (MQ=255)
                              ggTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCgg                                                                                                                       <  1:211740/149‑1 (MQ=255)
                               gTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCggg                                                                                                                      >  2:127535/1‑149 (MQ=255)
                                 aTTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGaa                                                                                                                    <  4:524981/149‑1 (MQ=255)
                                 aTTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGaa                                                                                                                    <  1:468983/149‑1 (MQ=255)
                                 aTTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGaa                                                                                                                    <  2:742443/149‑1 (MQ=255)
                                  ttttGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAAc                                                                                                                   >  1:770578/1‑149 (MQ=255)
                                  ttttGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAAc                                                                                                                   <  2:75202/149‑1 (MQ=255)
                                   tttGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACt                                                                                                                  <  3:68946/149‑1 (MQ=255)
                                         gttTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGa                                                                                                            <  4:486631/149‑1 (MQ=255)
                                             aTTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                        <  4:558078/149‑1 (MQ=255)
                                               ttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgtt                                                                                                                                        <  3:209399/115‑3 (MQ=255)
                                               ttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgtt                                                                                                                                        >  4:209399/1‑113 (MQ=255)
                                                tttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCg                                                                                                     >  3:85484/1‑149 (MQ=255)
                                                   ataGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTc                                                                                                  <  3:482663/149‑1 (MQ=255)
                                                    taGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCg                                                                                                 >  2:487391/1‑149 (MQ=255)
                                                    taGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCg                                                                                                 >  3:385649/1‑149 (MQ=255)
                                                    taGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCg                                                                                                 >  2:263836/1‑149 (MQ=255)
                                                     aGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGa                                                                                                <  2:292545/149‑1 (MQ=255)
                                                          tagaGCGTGTATTAACGGTGCTGTGGGAGAATCGTCAGAATTTGCCAGTGGCATGAGTCAACATTGCGGAGAAATAAAAACGCCTCTCAGTGTGGAGAGGCGCAGGAGATTACGCGTCAGAGCACTCACGGGTAGAGCGTGGGACAtgt                                                                                           >  2:191461/4‑147 (MQ=255)
                                                            cAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTc                                                                                         <  3:756689/149‑1 (MQ=255)
                                                                 tGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCa                                                                                    <  1:19630/149‑1 (MQ=255)
                                                                    aTTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCg                                                                                                                        >  2:435033/1‑110 (MQ=255)
                                                                    aTTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCg                                                                                                                        <  1:435033/110‑1 (MQ=255)
                                                                     ttAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGgaga                                                                                                                                   <  4:137228/98‑1 (MQ=255)
                                                                     ttAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGgaga                                                                                                                                   >  3:137228/1‑98 (MQ=255)
                                                                     ttAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATg                                                                                <  4:455097/149‑1 (MQ=255)
                                                                        aCGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTc                                                                             <  2:651116/149‑1 (MQ=255)
                                                                          ggTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGt                                                                           >  4:31517/1‑149 (MQ=255)
                                                                             gCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGc                                                                        <  1:700393/149‑1 (MQ=255)
                                                                              cTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGcc                                                                       <  2:770578/149‑1 (MQ=255)
                                                                                 tgGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGAc                                                                    >  2:119194/1‑149 (MQ=255)
                                                                                    gagaATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACaaa                                                                 >  4:741655/1‑149 (MQ=255)
                                                                                       aaTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                              <  1:823254/149‑1 (MQ=255)
                                                                                       aaTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                              <  3:227290/149‑1 (MQ=255)
                                                                                        aTCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAAt                                                             >  4:186732/1‑149 (MQ=255)
                                                                                         tCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATg                                                            >  1:736970/1‑149 (MQ=255)
                                                                                           gTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGaa                                                          <  4:190423/149‑1 (MQ=255)
                                                                                             cAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCTGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACg                                                        >  2:103527/1‑149 (MQ=255)
                                                                                               gagTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgc                                                      >  3:331765/1‑149 (MQ=255)
                                                                                               gagTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgc                                                      >  1:591205/1‑149 (MQ=255)
                                                                                                agTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgcg                                                     >  4:289549/1‑149 (MQ=255)
                                                                                                    cGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTgg                                                 >  3:419384/1‑149 (MQ=255)
                                                                                                      ccAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTg                                               >  4:546544/1‑149 (MQ=255)
                                                                                                            aCAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGt                                         >  2:788260/1‑149 (MQ=255)
                                                                                                                 agTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCg                                    <  4:548191/149‑1 (MQ=255)
                                                                                                                 agTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCATTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCg                                    >  4:780066/1‑149 (MQ=255)
                                                                                                                  gTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                              <  2:863479/122‑1 (MQ=255)
                                                                                                                  gTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                              >  1:863479/1‑122 (MQ=255)
                                                                                                                    cAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGga                                 <  3:186732/149‑1 (MQ=255)
                                                                                                                    cAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGGGCAGGGTTTCCGGga                                 >  1:436209/1‑149 (MQ=255)
                                                                                                                     aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGCGAGGCGCAGGAGATTACGCGTCCGGGAACTCACTGATAAAGCGTTCGACATCTTCAACCATATCGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                <  2:100089/149‑1 (MQ=255)
                                                                                                                     aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                >  3:260414/1‑149 (MQ=255)
                                                                                                                     aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                <  1:263836/149‑1 (MQ=255)
                                                                                                                             gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                        >  3:778263/1‑69 (MQ=255)
                                                                                                                             gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                        <  4:778263/69‑1 (MQ=255)
                                                                                                                             gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                        >  3:50379/1‑149 (MQ=255)
                                                                                                                             gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                        <  2:736970/149‑1 (MQ=255)
                                                                                                                             gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                        >  2:597266/1‑149 (MQ=255)
                                                                                                                                      aaaaCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGct               <  1:401730/149‑1 (MQ=255)
                                                                                                                                      aaaaCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGct               >  1:626907/1‑149 (MQ=255)
                                                                                                                                         aCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCtt            <  4:260414/149‑1 (MQ=255)
                                                                                                                                           gCCTCTCCACGCGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                        <  3:809276/135‑1 (MQ=255)
                                                                                                                                           gCCTCTCCACGCGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                        >  4:809276/1‑135 (MQ=255)
                                                                                                                                            ccTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGc         <  4:85484/149‑1 (MQ=255)
                                                                                                                                             ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGcc        <  1:453175/149‑1 (MQ=255)
                                                                                                                                             ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGcc        >  1:373333/1‑149 (MQ=255)
                                                                                                                                              tctcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCa       >  3:158056/1‑149 (MQ=255)
                                                                                                                                                ttcGTGTGGAGAGGCGCAGGAGAGTACGCGTCCGGGATCTCACGGATAAAGCGTTCGACATCTTCAACCTTATGGTCGTTGCCGACAAAGATAGACCGGCGCTGGTGCAGGATTACCGGGCGGATATCCAGAATACGGTCATTGCCCtc     <  2:436209/147‑1 (MQ=255)
                                                                                                                                                   gtgtGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtcgct  <  1:487391/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCT  >  NC_000913/4454444‑4454739

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