Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,419,742 G→A 5.8% A16V (GCG→GTG)  yfcE ← phosphodiesterase activity on bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,419,7420GA 5.8% 219.3 / 2.4 69A16V (GCG→GTG) yfcEphosphodiesterase activity on bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base A (2/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCCACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCTGGACAGA  >  NC_000913/2419598‑2419885
                                                                                                                                                |                                                                                                                                               
gCCACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCACCgg                                                                                                                                             >  3:42840/1‑149 (MQ=255)
gCCACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCACCgg                                                                                                                                             >  1:61794/1‑149 (MQ=255)
 ccACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCgg                                                                                                                                             >  6:336769/1‑148 (MQ=255)
 ccACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCgg                                                                                                                                             >  6:232097/1‑148 (MQ=255)
             gCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATg                                                                                                                                 >  1:53593/1‑148 (MQ=255)
              ccGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGaa                                                                                                                               <  1:209990/149‑1 (MQ=255)
               cGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAAt                                                                                                                              <  5:317729/149‑1 (MQ=255)
                  tCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTAGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGAAACGACCCATGAAAGTc                                                                                                                           >  8:38366/1‑149 (MQ=255)
                     gCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGa                                                                                                                        >  2:274842/1‑149 (MQ=255)
                                 ggcGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGG‑AACGACCCATGAATGTCCGATGCAAACATCAGt                                                                                                           >  7:91862/1‑149 (MQ=255)
                                  gcgcGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAACCAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCa                                                                                                             <  8:166021/147‑1 (MQ=255)
                                  gcgcGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGt                                                                                                           <  4:160934/149‑1 (MQ=255)
                                    gcgTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGttt                                                                                                         <  4:10288/149‑1 (MQ=255)
                                         aCCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTtcatc                                                                                                     <  2:53593/148‑1 (MQ=255)
                                            cTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCAtcactc                                                                                                 <  7:306020/149‑1 (MQ=255)
                                             tCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTGTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCAtcactca                                                                                                >  5:281027/1‑149 (MQ=255)
                                               cGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCAtcactca                                                                                                <  5:336769/147‑1 (MQ=255)
                                                 gTAAAGCATTACGCGGGCCGTGAGTCCACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCAGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCActcc                                                                                            <  4:62870/149‑1 (MQ=255)
                                                    aaGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCtcg                                                                                         >  3:113688/1‑149 (MQ=255)
                                                     aGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCtcgt                                                                                        <  1:274842/149‑1 (MQ=255)
                                                     aGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACAGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCtcgt                                                                                        >  2:215394/1‑149 (MQ=255)
                                                      gCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCtcgtc                                                                                       <  4:113688/149‑1 (MQ=255)
                                                             gcgGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGaaaaaa                                                                                >  2:210706/1‑149 (MQ=255)
                                                             gcgGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGaaaaaa                                                                                <  5:232097/149‑1 (MQ=255)
                                                              cgGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAAc                                                                               <  3:307768/149‑1 (MQ=255)
                                                                    gTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGc                                                                         >  1:143524/1‑149 (MQ=255)
                                                                    gTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGc                                                                         >  8:30215/1‑149 (MQ=255)
                                                                     tGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTtca                                                                                                       >  3:61243/1‑118 (MQ=255)
                                                                     tGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTtca                                                                                                       <  4:61243/118‑1 (MQ=255)
                                                                            acacGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATAc                                                                 >  4:92929/1‑149 (MQ=255)
                                                                                gTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGa                                                             >  4:255059/1‑149 (MQ=255)
                                                                                   gCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTc                                                          <  8:27532/149‑1 (MQ=255)
                                                                                     cAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTg                                                        <  1:210706/149‑1 (MQ=255)
                                                                                     cAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTACCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTg                                                        >  8:70393/1‑149 (MQ=255)
                                                                                       aGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACCTCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTc                                                      <  3:33722/149‑1 (MQ=255)
                                                                                       aGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGAAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTc                                                      <  2:124052/149‑1 (MQ=255)
                                                                                          aTCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGtt                                                   <  7:30215/149‑1 (MQ=255)
                                                                                           tCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTg                                                  >  8:181080/1‑149 (MQ=255)
                                                                                           tCACCAGACACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTg                                                  <  6:247588/149‑1 (MQ=255)
                                                                                                aGCCACTGGGCACCGCTTTGGGCAAACAACTACAGACCACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATc                                             <  2:301130/149‑1 (MQ=255)
                                                                                                 gCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCa                                            <  3:92929/149‑1 (MQ=255)
                                                                                                    aCTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTg                                                               >  8:23394/1‑127 (MQ=255)
                                                                                                    aCTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTg                                                               <  7:23394/127‑1 (MQ=255)
                                                                                                     cTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGc                                          >  1:151010/1‑147 (MQ=255)
                                                                                                        ggCAGCGCGTTGGGCAAACAACTCCCGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCCAACATCAGTTTCATCACTCACTCCGCGTCGAAAAAACGTACGCTAGGAGACTGGAGTGGGTCGTTGCAATCAGCCAGag                                     >  2:159889/1‑148 (MQ=255)
                                                                                                           aCCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTCGATTCTGTCGTTGCTATCAGCCAGATTGc                                  >  7:349923/1‑149 (MQ=255)
                                                                                                             cGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGAtt                                                           <  7:264349/122‑1 (MQ=255)
                                                                                                             cGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGAtt                                                           >  8:264349/1‑122 (MQ=255)
                                                                                                             cGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCtt                                >  4:262124/1‑149 (MQ=255)
                                                                                                                     gCAAACAACTCCAGAACACGTTCCGTCTCCGGTAACGACCCATGAATTTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGaaa                        <  3:67841/149‑1 (MQ=255)
                                                                                                                     gCAAACAACTCCAGAACACGTACCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGaaa                        <  7:70393/149‑1 (MQ=255)
                                                                                                                      cAAACAACTCCAGAACACGTTCCGTCGCCGG‑AACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGt                      <  8:91862/149‑1 (MQ=255)
                                                                                                                      cAAACAACTCCAGAACACGTTCCGTCACCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAg                       <  2:61794/149‑1 (MQ=255)
                                                                                                                      cAAACAACTCCAGAACACGTTCCGTCACCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAg                       <  4:42840/149‑1 (MQ=255)
                                                                                                                         acaacTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCaaa                                                                                                                  >  3:191663/1‑55 (MQ=255)
                                                                                                                         acaacTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCaaa                                                                                                                  <  4:191663/55‑1 (MQ=255)
                                                                                                                         acaacTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGc                    <  1:98451/149‑1 (MQ=255)
                                                                                                                              tCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGt                      <  6:57800/142‑1 (MQ=255)
                                                                                                                              tCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGt                      >  5:57800/1‑142 (MQ=255)
                                                                                                                              tCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATa                >  4:129928/1‑148 (MQ=255)
                                                                                                                              tCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGAAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATa                >  1:107928/1‑148 (MQ=255)
                                                                                                                                cAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATaaaa             <  8:161693/149‑1 (MQ=255)
                                                                                                                                cAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATaaaa             >  6:55098/1‑149 (MQ=255)
                                                                                                                                cAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATaa               <  6:285187/147‑1 (MQ=255)
                                                                                                                                    acacGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCt         <  6:281027/149‑1 (MQ=255)
                                                                                                                                    acacGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCt         <  8:44950/149‑1 (MQ=255)
                                                                                                                                     ccctttCCGTCGCCGTTAACGACCCATGAATGTCCGTTGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAATCTATGATACTGGATTCTGGCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCTg        <  8:76485/145‑1 (MQ=255)
                                                                                                                                         ttCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCTGGa      >  4:244235/1‑147 (MQ=255)
                                                                                                                                            cGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCTGGACAGa  <  3:255059/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
GCCACTTCATTAAGCCGTTCAGCGACTTTGGCTGGCGCGTAACCCTCCGGTAAAGCATTACGCGGGCCGTGATTCAACACGTCGCCAAGGATCACCAGCCACTGGGCACCGCTTTGGGCAAACAACTCCAGAACACGTTCCGTCGCCGGTAACGACCCATGAATGTCCGATGCAAACATCAGTTTCATCACTCACTCCTCGTCGAAAAAACGTAAGCTATGATACTGGATTCTGTCGTTGCTATCAGCCAGATTGCTTCGCCGAAAGTGCAATAAAACGCTGGACAGA  >  NC_000913/2419598‑2419885

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