Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,224 G→T 7.2% G82V (GGT→GTT)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,2240GT 7.2% 277.7 / 10.6 97G82V (GGT→GTT) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (40/50);  new base T (4/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

ATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGCGCTCT  >  NC_000913/4280076‑4280373
                                                                                                                                                    |                                                                                                                                                     
aTGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTGCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGt                                                                                                                                                      <  1:1802727/150‑1 (MQ=255)
  gCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCg                                                                                                                                                    >  2:1773413/1‑150 (MQ=255)
   cAAATCGCCATCGGTGCGCTTCTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCga                                                                                                                                                   <  1:818501/150‑1 (MQ=255)
   cAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCga                                                                                                                                                   >  2:1702159/1‑150 (MQ=255)
    aaaTCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgag                                                                                                                                                  >  2:215276/1‑150 (MQ=255)
      aTCGCCAGCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGGGAATTTCTCGAACATGGTCgagag                                                                                                                                                <  1:285014/150‑1 (MQ=255)
         gCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCgagagagt                                                                                                                                             >  2:1564426/1‑148 (MQ=255)
          ccATCGGTGCTCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCTCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGAttt                                                                                                                                            <  2:1522628/150‑1 (MQ=255)
          ccATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttc                                                                                                                                            <  2:427784/150‑2 (MQ=255)
           cATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTATGCATGTGGAGTGTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGTTCGAGAAAttca                                                                                                                                           <  1:180668/150‑3 (MQ=255)
              cGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATCTCTCGACCATGTTCGAGAAAttcacgg                                                                                                                                        <  2:1649332/150‑6 (MQ=255)
                gTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGGGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGcc                                                                                                                                      <  2:570127/150‑1 (MQ=255)
                   cgcTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTcg                                                                                                                                   >  2:1065701/1‑150 (MQ=255)
                     ctacccgCGTGCCCGACCTTTGGCCCGCATGTGGAGTTAGCTCCTGAACTCTTTTAAGCCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTcgcg                                                                                                                                 <  2:1346080/144‑1 (MQ=255)
                          ggCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggc                                                                                                                            >  2:815278/1‑150 (MQ=255)
                           gCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggcg                                                                                                                           <  2:1670025/150‑1 (MQ=255)
                           gCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCgcgctggcg                                                                                                                           <  2:1172701/150‑1 (MQ=255)
                               ggCCGACGTTTGGTTGGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGGTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGcgctggcgctgg                                                                                                                       <  1:219317/150‑1 (MQ=255)
                                  cGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtgg                                                                                                                    <  2:1691443/150‑1 (MQ=255)
                                     cGTTTGGTTTGCATGTGGAGTTTGCTCCTGAACTCTTTTTAGTCTTGTTTCTCCCGCCGTTGCTGTTCGCTGATGGCTGGAACACGCCGACCCGGGAATTTCTCGAACATGGTCGAGAGATTTTCTGCCTCGCGCTGGCGCtggtggtgg                                                                                                                 <  1:27670/150‑1 (MQ=255)
                                     cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                                 >  1:703787/1‑150 (MQ=255)
                                     cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                                 <  1:1660533/150‑1 (MQ=255)
                                      gtttggtttgCATGTGGAGTTTGACCCTGAACTCTTTTTTGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAACCGCCGACCCGTTAATTTCTCGATCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTggtggtggt                                                                                                                <  1:1554459/150‑1 (MQ=255)
                                         tggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCAc                                                                                                             >  1:1071988/1‑150 (MQ=255)
                                          ggtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCAGCCGTTTCTGTTCTCTGATGGCTGGAAAAATCCGACCCGTGATTTTATCAAACATTGTCTAGAGATTTTCGTCCTCGCGCTGGATCAGGTGATGTTCaac                                                                                                            >  1:389159/1‑148 (MQ=255)
                                           gtttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCg                                                                                                           >  1:832840/1‑150 (MQ=255)
                                            tttgCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGt                                                                                                          <  2:1522211/150‑1 (MQ=255)
                                               gCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGATCGCTAATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGt                                                                                                       <  2:636223/150‑1 (MQ=255)
                                               gCATGTGGAGTGTGATCCTGAACTCTTTTTTGTCGTGTTTATGCCGCCGTTGCGGGGCGCTGATGGCGGGAAAACGCCGACCCGTGAAGTTCGCGATCATGGTCGAGAGATGTTCGGCCTCGCGCTGGCGCTGGGGGTGGTCACCGTGGt                                                                                                       <  1:289996/150‑1 (MQ=255)
                                                cATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTc                                                                                                      >  1:939524/1‑150 (MQ=255)
                                                 aTGTGGAGGTAGATCCTGAACTCTTTCTCGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCg                                                                                                     <  2:331824/150‑1 (MQ=255)
                                                  tgtgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCgg                                                                                                    <  2:1886190/150‑1 (MQ=255)
                                                    tgGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCa                                                                                                  >  1:1934315/1‑150 (MQ=255)
                                                       aGTTTGCTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGCTGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTg                                                                                               <  1:1018959/150‑1 (MQ=255)
                                                        gTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACACGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGACGGCATTgg                                                                                              >  1:569130/1‑150 (MQ=255)
                                                          ttGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCt                                                                                            <  2:838501/150‑1 (MQ=255)
                                                          ttGATCCAGACCTCTTTTTAGGCTTGTTTATCCCGCCGTTGCTGTACGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGGGGTCACCGTGGTCGGCATTGGCt                                                                                            <  2:296547/150‑1 (MQ=255)
                                                            gATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTc                                                                                          >  1:1144858/1‑150 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  2:1019555/150‑1 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  2:699881/150‑1 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  2:692624/150‑1 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  2:1838791/150‑1 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGGGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  2:1060285/150‑1 (MQ=255)
                                                             aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTGCGCTGATGGCTGGAAAACGCCGACCCGTGAATTGCTCGAACCTGGGCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                         <  1:1440677/150‑1 (MQ=255)
                                                               ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                       >  1:1913132/1‑150 (MQ=255)
                                                               ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                       <  1:262075/150‑1 (MQ=255)
                                                               ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                       >  1:509905/1‑150 (MQ=255)
                                                                ccgAACTCTCTTAAGGCTTGTTTATCCCGCCGTTGCTGATCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTa                                                                                      <  1:1261584/148‑1 (MQ=255)
                                                                     ctctTTTTAGTCTTGTTTCTCCCGCCGTTGCTGTTCGCTGAGGGCTGGAAACCGCCGACCCGGGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGAGCTGGCGCTGGTGGGGGTCACCGTGGTCGGCATTGGCTTCCTTATTTAc                                                                                 <  1:269226/150‑1 (MQ=255)
                                                                     ctctTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTAc                                                                                 >  1:1462366/1‑150 (MQ=255)
                                                                      tctTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACt                                                                                >  2:908208/1‑150 (MQ=255)
                                                                         ttttAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTggg                                                                             <  1:1305069/150‑1 (MQ=255)
                                                                          tttAGTCTTGATTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGTATTGGCTTCCTTATTTACTGggt                                                                            <  2:1659621/150‑1 (MQ=255)
                                                                              gTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGgtggtg                                                                        >  1:1913302/1‑150 (MQ=255)
                                                                              gTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGggggtg                                                                        >  1:734464/1‑150 (MQ=255)
                                                                              gTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCGGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGgtggtg                                                                        <  1:654191/150‑1 (MQ=255)
                                                                                     ttATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGACCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCa                                                                 >  2:974851/1‑150 (MQ=255)
                                                                                         cccGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTcc                                                             >  1:1462119/1‑150 (MQ=255)
                                                                                           cgccgTTGCTGTGCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGc                                                           <  1:532968/150‑1 (MQ=255)
                                                                                             ccgTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTg                                                         <  1:817581/150‑1 (MQ=255)
                                                                                                     gTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGc                                                 >  2:1380396/1‑150 (MQ=255)
                                                                                                       tCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCt                                               >  2:561573/1‑150 (MQ=255)
                                                                                                        cGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCtt                                              <  2:1915813/150‑1 (MQ=255)
                                                                                                        cGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCtt                                              <  2:1827030/150‑1 (MQ=255)
                                                                                                           tGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTgc                                           >  2:1433588/1‑150 (MQ=255)
                                                                                                               ggCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGGCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTg                                       <  1:815278/150‑1 (MQ=255)
                                                                                                                 cTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGACATTCCGCTGATCCCTGCCTTTGCGCTggc                                     >  2:1609614/1‑150 (MQ=255)
                                                                                                                  tGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCTGCCTCGCGCAGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATACCTGCCTTTGCGCTggcg                                    >  1:1175468/1‑150 (MQ=255)
                                                                                                                   ggAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcgg                                   <  1:89405/150‑1 (MQ=255)
                                                                                                                    gAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTggcggc                                  >  1:1765965/1‑150 (MQ=255)
                                                                                                                     aaaaCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGgcggcg                                 >  2:830075/1‑150 (MQ=255)
                                                                                                                        aCGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTg                              <  1:215276/150‑1 (MQ=255)
                                                                                                                          gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCt                            <  2:647209/150‑1 (MQ=255)
                                                                                                                          gCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCt                            >  2:615792/1‑150 (MQ=255)
                                                                                                                             gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtc                         >  2:1792421/1‑150 (MQ=255)
                                                                                                                             gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtc                         <  2:1144858/150‑1 (MQ=255)
                                                                                                                             gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtc                         <  2:1553657/150‑1 (MQ=255)
                                                                                                                                ccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                      >  1:1686199/1‑150 (MQ=255)
                                                                                                                                ccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                      <  1:1168882/150‑1 (MQ=255)
                                                                                                                                ccGTGAATTTCTCGAACATGGGCGAGAGCTTTTCGGCCTCGCGCTGGCGCTGGGGGTGGGCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                      <  1:241274/150‑1 (MQ=255)
                                                                                                                                  gTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccga                    >  1:42234/1‑150 (MQ=255)
                                                                                                                                   tGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCGTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgac                   <  2:1430341/150‑1 (MQ=255)
                                                                                                                                   tGAAGTTCTGCATCATGGTCTTGAGATGTTCCGCCGCGCGCTGCCGCGGGCGGTGGTCACCGTGGTCGGCTTTTGATTCCATATTTACTGGGTCGTGCCGGGCATTCGGCTGATCCCTGCCTTCGCGCTGGCGGCGGTGCTTTCTccgac                   <  2:245633/150‑1 (MQ=255)
                                                                                                                                   aaaaTCTCTCGCCCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgac                   >  2:176213/3‑150 (MQ=255)
                                                                                                                                    gAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgacc                  >  2:1664229/1‑150 (MQ=255)
                                                                                                                                    gAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCGGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgacc                  <  2:75013/150‑1 (MQ=255)
                                                                                                                                    aaaTCTCTCGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgacc                  >  2:651624/2‑150 (MQ=255)
                                                                                                                                      atctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTccgaccgc                >  2:393073/4‑149 (MQ=255)
                                                                                                                                       tctctcGACCATGTTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGAt               >  1:1913174/3‑150 (MQ=255)
                                                                                                                                        ttCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATg              >  1:305232/1‑150 (MQ=255)
                                                                                                                                         tctcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTCATACCTGCCTTTGCTCTGGCGGCGGTGCTTTCTCCGACCGATGc             >  1:99355/1‑150 (MQ=255)
                                                                                                                                         gctcGAACCGGGGCGAGAGATGTTCGGCCTCGCGCTGGCGCTGGGGGTGGGCACCGTGGTCGGCATTGGCTTCCTTATATACTGGGTGGGGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGc             <  1:1360336/149‑1 (MQ=255)
                                                                                                                                          ctcGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCt            >  1:685810/1‑150 (MQ=255)
                                                                                                                                          ctcGAACATGGTCGAGAGATTTTCGACCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCt            >  1:825587/1‑150 (MQ=255)
                                                                                                                                              aaCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTgg        >  1:421256/1‑150 (MQ=255)
                                                                                                                                               aCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGCCCGATGCTGTGgc       >  1:677612/1‑150 (MQ=255)
                                                                                                                                               aCATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGgc       <  1:1113105/150‑1 (MQ=255)
                                                                                                                                                   ggTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGCGctc   >  1:933988/1‑150 (MQ=255)
                                                                                                                                                    gTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGAATTCAGCTGATCCCTGCCTTTGCCCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGCGctct  >  2:1405612/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGCGCTCT  >  NC_000913/4280076‑4280373

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