Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,454,594 T→C 10.7% intergenic (+16/+17) 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,5940TC10.7% 231.5 / 16.5 84intergenic (+16/+17)yjfF/fbpputative sugar ABC transporter permease/fructose‑1,6‑bisphosphatase I
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (2/7);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCTCG  >  NC_000913/4454446‑4454741
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                               <  3:227581/149‑1 (MQ=255)
      cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                               >  3:264484/1‑149 (MQ=255)
      cTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGagag                                                                                                                                               >  2:468983/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                                                                                                                                           <  1:178731/149‑2 (MQ=255)
          tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                           <  2:23000/149‑2 (MQ=255)
          tggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgt                                                                                                                                           <  4:635871/149‑2 (MQ=255)
            gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                         <  2:46700/149‑1 (MQ=255)
            gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                         <  3:173654/149‑1 (MQ=255)
            gACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg                                                                                                                                         <  4:256891/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                                                                                                                      <  1:468983/149‑1 (MQ=255)
                               aTTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGaa                                                                                                                      <  4:524981/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                                                                                                                                          >  4:209399/1‑113 (MQ=255)
                                             ttttATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCACACGGAGAGgcgtt                                                                                                                                          <  3:209399/115‑3 (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                                                        >  1:591205/1‑149 (MQ=255)
                                                                                             gagTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGgc                                                        >  3:331765/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                                                                >  1:863479/1‑122 (MQ=255)
                                                                                                                gTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGaa                                                                <  2:863479/122‑1 (MQ=255)
                                                                                                                  cAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGga                                   <  3:186732/149‑1 (MQ=255)
                                                                                                                  cAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGGGCAGGGTTTCCGGga                                   >  1:436209/1‑149 (MQ=255)
                                                                                                                   aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGCGAGGCGCAGGAGATTACGCGTCCGGGAACTCACTGATAAAGCGTTCGACATCTTCAACCATATCGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                  <  2:100089/149‑1 (MQ=255)
                                                                                                                   aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                  <  1:263836/149‑1 (MQ=255)
                                                                                                                   aaCATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGgat                                  >  3:260414/1‑149 (MQ=255)
                                                                                                                           gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                          <  4:778263/69‑1 (MQ=255)
                                                                                                                           gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATaaa                                                                                                          >  3:778263/1‑69 (MQ=255)
                                                                                                                           gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                          >  2:597266/1‑149 (MQ=255)
                                                                                                                           gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                          <  2:736970/149‑1 (MQ=255)
                                                                                                                           gcagcaATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                          >  3:50379/1‑149 (MQ=255)
                                                                                                                                    aaaaCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGct                 >  1:626907/1‑149 (MQ=255)
                                                                                                                                    aaaaCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGct                 <  1:401730/149‑1 (MQ=255)
                                                                                                                                       aCGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCtt              <  4:260414/149‑1 (MQ=255)
                                                                                                                                         gCCTCTCCACGCGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                          >  4:809276/1‑135 (MQ=255)
                                                                                                                                         gCCTCTCCACGCGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCa                          <  3:809276/135‑1 (MQ=255)
                                                                                                                                          ccTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGc           <  4:85484/149‑1 (MQ=255)
                                                                                                                                           ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGcc          >  1:373333/1‑149 (MQ=255)
                                                                                                                                           ctctCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGcc          <  1:453175/149‑1 (MQ=255)
                                                                                                                                            tctcCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCa         >  3:158056/1‑149 (MQ=255)
                                                                                                                                              ttcGTGTGGAGAGGCGCAGGAGAGTACGCGTCCGGGATCTCACGGATAAAGCGTTCGACATCTTCAACCTTATGGTCGTTGCCGACAAAGATAGACCGGCGCTGGTGCAGGATTACCGGGCGGATATCCAGAATACGGTCATTGCCCtc       <  2:436209/147‑1 (MQ=255)
                                                                                                                                                 gtgtGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtcgct    <  1:487391/149‑1 (MQ=255)
                                                                                                                                                   gtgGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCAtcgctcg  <  2:626907/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAGGAGATTACGCGTCCGGGAACTCACGGATAAAGCGTTCGACATCTTCAACCATATGGTCGTTGCCGACAAAGAATGAACGGCGCTGGTGCAGGGTTTCCGGGATGATATCCAGAATACGCTCTTTGCCATCGCTCG  >  NC_000913/4454446‑4454741

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: