Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,280,219 A→C 8.1% E80D (GAA→GAC yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,280,2190AC 8.1% 239.6 / 10.8 96E80D (GAA→GACyjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (38/47);  new base C (5/3);  total (43/53)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280070‑4280368
                                                                                                                                                     |                                                                                                                                                     
ccgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGaa                                                                                                                                                       <  2:1212549/150‑1 (MQ=255)
 cgcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAAc                                                                                                                                                      >  1:1469244/1‑150 (MQ=255)
  gcttATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACa                                                                                                                                                     <  1:1292242/150‑1 (MQ=255)
      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:1660533/150‑1 (MQ=255)
                                           cGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCtggtggtgg                                                                                                            >  1:703787/1‑150 (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:699881/150‑1 (MQ=255)
                                                                   aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                    <  2:1019555/150‑1 (MQ=255)
                                                                   aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                    <  2:1838791/150‑1 (MQ=255)
                                                                   aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                    <  2:692624/150‑1 (MQ=255)
                                                                   aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGGGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                    <  2:1060285/150‑1 (MQ=255)
                                                                   aTCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTGCGCTGATGGCTGGAAAACGCCGACCCGTGAATTGCTCGAACCTGGGCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTcc                                                                                    <  1:1440677/150‑1 (MQ=255)
                                                                     ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                  <  1:262075/150‑1 (MQ=255)
                                                                     ccTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCtt                                                                                  >  1:1913132/1‑150 (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:1553657/150‑1 (MQ=255)
                                                                                                                                   gACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTtc                    <  2:1144858/150‑1 (MQ=255)
                                                                                                                                      ccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                 <  1:1168882/150‑1 (MQ=255)
                                                                                                                                      ccGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTcc                 >  1:1686199/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                     
CCGCTTATGCAAATCGCCATCGGTGCGCTACTGGCGTGGCCGACGTTTGGTTTGCATGTGGAGTTTGATCCTGAACTCTTTTTAGTCTTGTTTATCCCGCCGTTGCTGTTCGCTGATGGCTGGAAAACGCCGACCCGTGAATTTCTCGAACATGGTCGAGAGATTTTCGGCCTCGCGCTGGCGCTGGTGGTGGTCACCGTGGTCGGCATTGGCTTCCTTATTTACTGGGTGGTGCCGGGCATTCCGCTGATCCCTGCCTTTGCGCTGGCGGCGGTGCTTTCTCCGACCGATGCTGTGGC  >  NC_000913/4280070‑4280368

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: