Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,092,780 C→A 5.8% G15V (GGC→GTC)  frvA ← putative enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7800CA 5.8% 333.2 / 6.5 103G15V (GGC→GTC) frvAputative enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (44/53);  new base A (4/2);  total (48/55)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 106 total aligned reads displayed.

GTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAA  >  NC_000913/4092632‑4092928
                                                                                                                                                     |                                                                                                                                                    
gTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGc                                                                                                                                                      <  2:413866/149‑1 (MQ=255)
  gCGGCTCGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                   <  2:1209110/149‑1 (MQ=255)
  gCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGccc                                                                                                                                                    <  2:1515372/149‑1 (MQ=255)
   cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                   <  1:1067002/149‑1 (MQ=255)
   cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                   <  1:243020/149‑1 (MQ=255)
   cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGGGTGAGTCGGGGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                   <  2:179983/149‑1 (MQ=255)
      cACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGaa                                                                                                                                                <  1:452170/149‑1 (MQ=255)
      cACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGaa                                                                                                                                                >  1:370272/1‑149 (MQ=255)
      cACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGaa                                                                                                                                                <  1:1325400/149‑1 (MQ=255)
          gCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATAtt                                                                                                                                            <  2:370272/149‑1 (MQ=255)
           cGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTc                                                                                                                                           <  2:943793/149‑1 (MQ=255)
            gACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCa                                                                                                                                          <  2:1331416/149‑1 (MQ=255)
             acacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAg                                                                                                                                         <  2:226357/149‑1 (MQ=255)
             acacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAg                                                                                                                                         <  2:1176245/149‑1 (MQ=255)
             acacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGGCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAg                                                                                                                                         <  1:1058001/149‑1 (MQ=255)
               acCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGt                                                                                                                                       >  1:1536721/1‑149 (MQ=255)
                 cGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTGCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCa                                                                                                                                     >  2:970571/1‑149 (MQ=255)
                  ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCaa                                                                                                                                    <  1:938572/149‑1 (MQ=255)
                   gAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTGTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAAt                                                                                                                                   >  1:699111/1‑149 (MQ=255)
                    aaCCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATa                                                                                                                                  <  2:223336/149‑1 (MQ=255)
                     aCCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGACCTGAATATTCAGGGCAAtgg                                                                                                                                 <  1:160801/149‑3 (MQ=255)
                          aTCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCtt                                                                                                                            >  1:1360559/1‑149 (MQ=255)
                           tCCAGTAGAGTGCATTTTTTCGAGCAGCAGCAGTGTTTGCATGAACTGGTGTGAGTCTGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTg                                                                                                                           <  1:996352/149‑1 (MQ=255)
                             cAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCa                                                                                                                         >  2:1301732/1‑149 (MQ=255)
                               gttgagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGGTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTTCAGt                                                                                                                       <  1:81223/146‑1 (MQ=255)
                                tAGAGTGCATTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGt                                                                                                                       <  1:772877/149‑1 (MQ=255)
                                     tGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGATCAATACAGCTTGCAGTAAGAGc                                                                                                                 >  1:133474/1‑149 (MQ=255)
                                       cATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTg                                                                                                               >  1:7306/1‑149 (MQ=255)
                                          ttttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGATCGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcc                                                                                                             >  1:1562656/1‑149 (MQ=255)
                                          ttttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAg                                                                                                                                         <  1:271611/121‑1 (MQ=255)
                                          ttttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcc                                                                                                             >  1:338074/1‑149 (MQ=255)
                                          ttttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCATGAATATTCAg                                                                                                                                         >  2:271611/1‑121 (MQ=255)
                                           tttttCGCGCAGCAGCAGGGGTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATAGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCa                                                                                                            <  2:159352/149‑1 (MQ=255)
                                                 gcgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGCGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATc                                                                                                      >  1:714855/1‑149 (MQ=255)
                                                  cgcAGCAGCAGGGTTTGCAGGAACTGGGGTGAGGCGGTGATAAACCCGTTAGGTAACGCTATTGACGCCACCTTCTTCAGACCGGAATAAGCGCCATTGCCCTGAATATTCAGGGCACTACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                     <  1:230209/149‑1 (MQ=255)
                                                    cagcagcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGa                                                                                                   <  1:1472806/149‑1 (MQ=255)
                                                        agcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTcc                                                                                               >  2:653915/1‑149 (MQ=255)
                                                        agcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTcc                                                                                               <  2:1396033/149‑1 (MQ=255)
                                                           agGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATa                                                                                            >  1:659755/1‑149 (MQ=255)
                                                              gTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaaa                                                                                         <  2:89384/149‑1 (MQ=255)
                                                              gTTGGCAGGAACTGGTGTGAGTCGGGGATAAAACCGTGTTGTAACGCTCTTGGCGCCAACTGCTTCAGAACGGACTATGCTCCTTTGCCCTTATTATTCAGCTCAATACAGCGTGCAGTAAGTGCTGCCATTAATCGGACTCCATaaaa                                                                                         <  2:258252/149‑1 (MQ=255)
                                                                  gCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatata                                                                                     <  1:1237600/149‑1 (MQ=255)
                                                                    aGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATAtt                                                                                   <  1:1422242/149‑1 (MQ=255)
                                                                         cTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAAt                                                                              <  2:553664/149‑1 (MQ=255)
                                                                           ggTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTatgga                                                                            <  2:1436559/149‑5 (MQ=255)
                                                                           ggTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTatgga                                                                            <  2:1166184/149‑5 (MQ=255)
                                                                            gtgtGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCaa                                                                           >  2:1527442/1‑149 (MQ=255)
                                                                            gtgtGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCaa                                                                           <  1:1527442/149‑1 (MQ=255)
                                                                             tgtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAt                                                                          <  2:402713/149‑1 (MQ=255)
                                                                                gAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGt                                                                       >  1:1188682/1‑149 (MQ=255)
                                                                                 aGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTa                                                                      >  2:1100191/1‑149 (MQ=255)
                                                                                   tCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAt                                                                    <  1:1086883/149‑1 (MQ=255)
                                                                                    cGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtg                                                                   <  1:1183556/149‑1 (MQ=255)
                                                                                        gATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtgg                                                               <  2:581757/149‑1 (MQ=255)
                                                                                        gATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtgg                                                               <  2:18972/149‑1 (MQ=255)
                                                                                         aTAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGa                                                              >  2:349008/1‑149 (MQ=255)
                                                                                            aaaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                           >  2:332059/1‑149 (MQ=255)
                                                                                             aaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGa                                                          >  1:826103/1‑149 (MQ=255)
                                                                                              aCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGAt                                                         >  1:280473/1‑149 (MQ=255)
                                                                                               ccGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGAtt                                                        >  1:934664/1‑149 (MQ=255)
                                                                                                cGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTa                                                       <  1:332059/149‑1 (MQ=255)
                                                                                                  ttttGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAg                                                     >  1:1144202/1‑149 (MQ=255)
                                                                                                   tttGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGt                                                    <  1:584228/149‑1 (MQ=255)
                                                                                                       tAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtgatg                                                <  2:1188682/149‑1 (MQ=255)
                                                                                                         aCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCATGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTAAATCAATTGTAATGGTGGAGCGATTAAGTGATgag                                              <  2:1037135/149‑1 (MQ=255)
                                                                                                          cGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATgaga                                             >  1:1080678/1‑149 (MQ=255)
                                                                                                             tatattCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGGCAATACAGCTTGCAGAATGAGCAGCCATAAATCGGACTCCATAAAATATATTTAAATCAAATGTAATGGTGAAGCGATTATCTGATGATAGGc                                          >  1:74122/6‑149 (MQ=255)
                                                                                                             tATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGc                                          >  1:1073131/1‑149 (MQ=255)
                                                                                                               ttGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGa                                        >  1:1448086/1‑149 (MQ=255)
                                                                                                                tGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGat                                       <  2:1360559/149‑1 (MQ=255)
                                                                                                                tGTCGCCAACTGCTTCAGAACGGAATAAGAGCCATTGCCCTGAATATTTAGGTCAATACAGCTTGCAGTAAGAGCTGACATTTATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCgtt                                       <  2:74122/149‑3 (MQ=255)
                                                                                                                            cTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAAc                           >  1:807536/1‑149 (MQ=255)
                                                                                                                             ttCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTATACa                          >  1:716059/1‑149 (MQ=255)
                                                                                                                             ttCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACa                          >  1:478587/1‑149 (MQ=255)
                                                                                                                               cAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGATCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                        <  2:133474/149‑1 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATg                       >  2:398163/1‑149 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATg                       >  2:29904/1‑149 (MQ=255)
                                                                                                                                 taaCGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      <  1:1038602/148‑1 (MQ=255)
                                                                                                                                 gAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      >  2:152399/1‑149 (MQ=255)
                                                                                                                                  aaCGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                     >  1:980240/1‑149 (MQ=255)
                                                                                                                                    cGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCATTTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGaaaa                   >  2:1386868/1‑149 (MQ=255)
                                                                                                                                      gAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAtt                 >  2:172351/1‑149 (MQ=255)
                                                                                                                                      gAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAtt                 <  1:1301732/149‑1 (MQ=255)
                                                                                                                                       aaTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                <  1:45578/149‑1 (MQ=255)
                                                                                                                                       aaTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                <  2:1080678/149‑1 (MQ=255)
                                                                                                                                        aTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTg               >  1:1596906/1‑149 (MQ=255)
                                                                                                                                        aTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTg               <  1:1514146/149‑1 (MQ=255)
                                                                                                                                         tAAGCGCCATTGCCGTGAATATTCAGTTCAATACAGATTTCAGTAAGTGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGAAATGGTGGATCGTTTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTga              <  2:342556/149‑1 (MQ=255)
                                                                                                                                         tAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTga              <  2:663292/149‑1 (MQ=255)
                                                                                                                                         tAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTga              <  1:1100191/149‑1 (MQ=255)
                                                                                                                                         tAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATCACAATTTga              >  2:1489468/1‑149 (MQ=255)
                                                                                                                                           aGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTgaga            >  2:441320/1‑149 (MQ=255)
                                                                                                                                            gcgcCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc           >  2:147702/1‑149 (MQ=255)
                                                                                                                                              tgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACag         >  2:1411557/4‑149 (MQ=255)
                                                                                                                                              tgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACag         >  1:1059336/4‑149 (MQ=255)
                                                                                                                                                taTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGt                               >  2:1259653/2‑125 (MQ=255)
                                                                                                                                                taTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGt                               <  1:1259653/124‑1 (MQ=255)
                                                                                                                                                  ttGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTg     >  2:1517312/1‑149 (MQ=255)
                                                                                                                                                  ttGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTg     >  1:1441588/1‑149 (MQ=255)
                                                                                                                                                     cccTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGaaa  >  2:986879/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAA  >  NC_000913/4092632‑4092928

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