Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,805,081 T→C G89G (GGT→GGC proV → glycine betaine ABC transporter ATP binding subunit ProV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,805,0810TC98.5% 225.0 / ‑5.1 65G89G (GGT→GGCproVglycine betaine ABC transporter ATP binding subunit ProV
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/32);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGG  >  NC_000913/2804953‑2805218
                                                                                                                                |                                                                                                                                         
gCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGc                                                                                                                                >  1:507377/1‑140 (MQ=255)
     tCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCaaaa                                                                                                                           >  2:442230/1‑140 (MQ=255)
      cTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAat                                                                                                                          >  2:725772/1‑140 (MQ=255)
       tGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGa                                                                                                                                      <  2:398632/127‑1 (MQ=255)
       tGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGa                                                                                                                                      >  1:398632/1‑127 (MQ=255)
           cATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCg                                                                                                                     <  1:442230/140‑1 (MQ=255)
             ttGAAGAAGGCGAGATACTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGCGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGAc                                                                                                                   <  2:1037096/140‑1 (MQ=255)
                  gaagGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGa                                                                                                              <  2:811972/140‑1 (MQ=255)
                          gaTATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTg                                                                                                      >  1:12346/1‑140 (MQ=255)
                                tGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATTCGACGCCGAACTCCGTGAGGTgc                                                                                                >  1:294925/1‑140 (MQ=255)
                                    atcatGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAg                                                                                            >  1:856863/1‑140 (MQ=255)
                                      catGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaa                                                                                          >  1:903305/1‑140 (MQ=255)
                                      catGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaa                                                                                          >  1:650701/1‑140 (MQ=255)
                                       atGGGATTATCCGGCCCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaaa                                                                                         >  1:1006237/1‑140 (MQ=255)
                                        tGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGaaaa                                                                                        >  2:794060/1‑140 (MQ=255)
                                                tCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGc                                                                                >  1:637971/1‑140 (MQ=255)
                                                tCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGc                                                                                >  2:207126/1‑140 (MQ=255)
                                                tCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGc                                                                                >  2:528275/1‑140 (MQ=255)
                                                tCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGc                                                                                <  2:650701/140‑1 (MQ=255)
                                                 ccGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCg                                                                               >  2:606407/1‑140 (MQ=255)
                                                  cGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGa                                                                              <  1:84696/140‑1 (MQ=255)
                                                    gCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATg                                                                            <  2:474631/140‑1 (MQ=255)
                                                     cTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATgg                                                                           <  1:207126/140‑1 (MQ=255)
                                                     cTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATgg                                                                           <  2:527543/140‑1 (MQ=255)
                                                           tAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTcc                                                                     <  2:158252/140‑1 (MQ=255)
                                                            aaaTCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCa                                                                    <  2:203879/140‑1 (MQ=255)
                                                                  acaATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCtt                                                              <  2:856882/140‑1 (MQ=255)
                                                                  acaATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCtt                                                              <  2:85783/140‑1 (MQ=255)
                                                                  acaATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCtt                                                              <  1:246772/140‑1 (MQ=255)
                                                                    aaTGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTg                                                            >  1:599113/1‑140 (MQ=255)
                                                                    aaTGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTg                                                            >  1:193693/1‑140 (MQ=255)
                                                                         tACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTa                                                       >  1:265274/1‑140 (MQ=255)
                                                                         tACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTa                                                       >  1:102870/1‑140 (MQ=255)
                                                                          aCGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTaa                                                      <  1:888958/140‑1 (MQ=255)
                                                                            gCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATg                                                    <  2:265274/140‑1 (MQ=255)
                                                                            gCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATg                                                    >  2:965724/1‑140 (MQ=255)
                                                                            gCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATg                                                    <  1:880759/140‑1 (MQ=255)
                                                                              cTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGcc                                                  >  2:327526/1‑140 (MQ=255)
                                                                                tctcAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGc                                                >  2:20323/1‑140 (MQ=255)
                                                                                tctcAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGc                                                >  1:645726/1‑140 (MQ=255)
                                                                                   cAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATg                                                    >  1:312614/1‑133 (MQ=255)
                                                                                   cAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATg                                                    <  2:312614/133‑1 (MQ=255)
                                                                                      tCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGa                                          >  2:101824/1‑140 (MQ=255)
                                                                                         ccTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCg                                       <  1:526736/140‑1 (MQ=255)
                                                                                          cTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGt                                      >  1:369565/1‑140 (MQ=255)
                                                                                           tGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTg                                     <  2:507377/140‑1 (MQ=255)
                                                                                           tGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTg                                     <  2:586857/140‑1 (MQ=255)
                                                                                           tGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTg                                     <  1:725772/140‑1 (MQ=255)
                                                                                           tGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTg                                     <  2:193693/140‑1 (MQ=255)
                                                                                            gATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGc                                    >  1:1118162/1‑140 (MQ=255)
                                                                                                gAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGa                                <  2:369565/140‑1 (MQ=255)
                                                                                                gAACCCACCCGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGa                                <  2:12346/140‑1 (MQ=255)
                                                                                                  acccacccGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACa                              >  2:326036/1‑140 (MQ=255)
                                                                                                      acccGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATAc                          >  2:301929/1‑140 (MQ=255)
                                                                                                       cccGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACt                         >  1:499218/1‑140 (MQ=255)
                                                                                                       cccGCGGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACt                         >  2:311962/1‑140 (MQ=255)
                                                                                                        ccGCGGGCAAGCGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGGGAGGAGCGGAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGGGCTGGTCAATAcgg                        >  2:776244/1‑138 (MQ=255)
                                                                                                          gcgGGCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCg                      <  1:327526/140‑1 (MQ=255)
                                                                                                            gggCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    <  1:528275/140‑1 (MQ=255)
                                                                                                            gggCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    <  2:1118162/140‑1 (MQ=255)
                                                                                                            gggCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    <  2:1006237/140‑1 (MQ=255)
                                                                                                            gggCAAGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    <  2:856863/140‑1 (MQ=255)
                                                                                                                 aGTGCTGATTGATGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTa               <  2:1164830/140‑1 (MQ=255)
                                                                                                                        attgatGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    <  1:403430/128‑1 (MQ=255)
                                                                                                                        attgatGGCGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGtt                    >  2:403430/1‑128 (MQ=255)
                                                                                                                              ggTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCgg  >  2:923972/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                         
GCCAGTCTGGCCATTGAAGAAGGCGAGATATTTGTCATCATGGGATTATCCGGCTCGGGTAAATCCACAATGGTACGCCTTCTCAATCGCCTGATTGAACCCACCCGCGGGCAAGTGCTGATTGATGGTGTGGATATTGCCAAAATATCCGACGCCGAACTCCGTGAGGTGCGCAGAAAAAAGATTGCGATGGTCTTCCAGTCCTTTGCCTTAATGCCGCATATGACCGTGCTGGACAATACTGCGTTCGGTATGGAATTGGCCGG  >  NC_000913/2804953‑2805218

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: