Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,623,539 T→C G47G (GGT→GGC yhhH → putative NTF2 fold immunity protein for polymorphic toxin RhsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,623,5390TC100.0% 318.7 / NA 82G47G (GGT→GGCyhhHputative NTF2 fold immunity protein for polymorphic toxin RhsB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (49/33);  total (49/33)

TACCTAATGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGA  >  NC_000913/3623393‑3623685
                                                                                                                                                  |                                                                                                                                                  
tACCTAATGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCaaa                                                                                                                                                 <  1:432084/150‑1 (MQ=255)
  ccTAATGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAatat                                                                                                                                              >  2:898393/1‑151 (MQ=255)
    tAATGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATAc                                                                                                                                            >  1:244642/1‑151 (MQ=255)
     aaTGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCGAATATAcc                                                                                                                                           <  2:185113/151‑1 (MQ=255)
     aaTGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATAc                                                                                                                                            <  1:467274/150‑1 (MQ=255)
        gTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGtt                                                                                                                                        <  2:792593/151‑1 (MQ=255)
        gTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGt                                                                                                                                         <  2:299962/150‑1 (MQ=255)
            gCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGa                                                                                                                                     >  1:234674/1‑150 (MQ=255)
             cAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATTCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGaaa                                                                                                                                   <  1:368968/151‑1 (MQ=255)
                    cATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGtt                                                                                                                            <  1:535181/151‑1 (MQ=255)
                    cATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGt                                                                                                                             <  1:77565/150‑1 (MQ=255)
                    cATGGTAATATGTAAAAAAATATTGTCTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATAGGCGCAGAGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGtt                                                                                                                            >  2:324017/1‑151 (MQ=255)
                       ggTAATATGTAAAAAAATATTGTTTAGGAACTGTTTCATTGTATCTTTGTTTGTTTTTACATACAAAGCATGGGCGAAGTGTAATAATAATATTAAAATAACGAGCAAGTATGAAAGTGAAGGCAAAGATTCCGTTAGAAGTGTGGTTaaa                                                                                                                         <  2:273093/151‑1 (MQ=255)
                       ggTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTaaa                                                                                                                         <  1:1051052/151‑1 (MQ=255)
                       ggTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTaaa                                                                                                                         <  2:806341/151‑1 (MQ=255)
                        gTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTaaa                                                                                                                         >  1:301070/1‑150 (MQ=255)
                          aaTATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAAt                                                                                                                      >  1:400789/1‑151 (MQ=255)
                          aaTATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAAt                                                                                                                      >  1:985048/1‑151 (MQ=255)
                           atatGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAAt                                                                                                                      >  1:339762/1‑150 (MQ=255)
                            tatGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTAAGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATa                                                                                                                     >  1:282382/1‑150 (MQ=255)
                             atGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATaaa                                                                                                                   <  2:476791/151‑1 (MQ=255)
                               gTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAg                                                                                                                  >  1:697221/1‑150 (MQ=255)
                                tAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCt                                                                                                                >  2:50938/1‑151 (MQ=255)
                                   aaaaaTATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCtat                                                                                                              >  1:648628/1‑150 (MQ=255)
                                           tGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTgg                                                                                                      >  2:487052/1‑150 (MQ=255)
                                               tAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAAtt                                                                                                  >  2:380384/1‑150 (MQ=255)
                                                aGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAg                                                                                                <  2:262756/151‑1 (MQ=255)
                                                    aCTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTg                                                                                             >  1:901929/1‑150 (MQ=255)
                                                          tCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGatat                                                                                       >  1:848420/1‑150 (MQ=255)
                                                           cATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGatatat                                                                                     >  2:214692/1‑151 (MQ=255)
                                                             ttGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGatatat                                                                                     <  2:131165/149‑1 (MQ=255)
                                                             ttGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGCGTAATACTAATATTAAAATAATGAGCAAGTATGAAAGCGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGt                                                                                   >  1:269344/1‑151 (MQ=255)
                                                                tATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATTATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTa                                                                                 >  2:276358/1‑150 (MQ=255)
                                                                tATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTa                                                                                 >  2:380198/1‑150 (MQ=255)
                                                                 aacTTTGATTGGTTTTACATACAACACAGGGGCGCCGTGTAATAATAATATTGAAGTAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTaa                                                                                >  2:83975/3‑150 (MQ=255)
                                                                  tCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTa                                                                                 >  2:560020/1‑148 (MQ=255)
                                                                     ttgtttgtTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAAt                                                                            <  2:400779/150‑1 (MQ=255)
                                                                      tgtttgttTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATc                                                                           <  1:98281/150‑1 (MQ=255)
                                                                       gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGt                                                                         >  1:349095/1‑151 (MQ=255)
                                                                       gtttgtttTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGt                                                                         <  2:514767/151‑1 (MQ=255)
                                                                              tttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGa                                                                   <  2:321042/150‑1 (MQ=255)
                                                                               ttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGAc                                                                  >  1:17562/1‑150 (MQ=255)
                                                                               ttACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGAc                                                                  <  1:851029/150‑1 (MQ=255)
                                                                                  catacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAgg                                                               >  1:240348/1‑150 (MQ=255)
                                                                                   atacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                             <  2:553308/151‑1 (MQ=255)
                                                                                   atacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGAt                                                             >  2:714957/1‑151 (MQ=255)
                                                                                    tacaACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATg                                                            >  2:720898/1‑151 (MQ=255)
                                                                                         cacaTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCa                                                       <  2:892471/151‑1 (MQ=255)
                                                                                            aTGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCaga                                                     >  1:184019/1‑150 (MQ=255)
                                                                                                   cAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAgaagaagaa                                              >  1:238083/1‑150 (MQ=255)
                                                                                                    aGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGaaa                                             >  2:828144/1‑150 (MQ=255)
                                                                                                    aGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATAACGTTAGAAATTTGGTTAAAAATAAAGATATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAgaagaagaa                                              >  2:194333/1‑149 (MQ=255)
                                                                                                      tgtAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAc                                          >  1:274845/1‑151 (MQ=255)
                                                                                                       gtAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAc                                          <  2:157297/150‑1 (MQ=255)
                                                                                                        taataataaTATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAAcc                                         >  1:32635/1‑150 (MQ=255)
                                                                                                         aataataatATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCa                                        <  2:173540/150‑1 (MQ=255)
                                                                                                         aataataatATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCa                                        <  2:360108/150‑1 (MQ=255)
                                                                                                         aataataatATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCAt                                       >  2:817749/1‑151 (MQ=255)
                                                                                                         aataataatATTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCAt                                       >  2:866570/1‑151 (MQ=255)
                                                                                                             ataataTTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACg                                    >  1:476642/1‑150 (MQ=255)
                                                                                                              taataTTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGa                                   >  1:1060964/1‑150 (MQ=255)
                                                                                                              taataTTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGa                                   >  1:85284/1‑150 (MQ=255)
                                                                                                              taataTTAAAATAATGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGa                                   >  2:896714/1‑150 (MQ=255)
                                                                                                                        ataatGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGt                        <  2:271651/151‑1 (MQ=255)
                                                                                                                        ataatGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGt                        <  1:712774/151‑1 (MQ=255)
                                                                                                                         taatGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGtt                       >  1:1047758/1‑151 (MQ=255)
                                                                                                                           atGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTa                      >  2:182736/1‑150 (MQ=255)
                                                                                                                           atGCGCAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTa                      >  1:213766/1‑150 (MQ=255)
                                                                                                                               gcAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTaacaac                 >  2:415432/1‑151 (MQ=255)
                                                                                                                                cAAGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAACCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAacaaca                <  2:185225/151‑1 (MQ=255)
                                                                                                                                  aGTATGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAacaaca                >  1:145129/1‑149 (MQ=255)
                                                                                                                                     aTGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGtt            <  2:303970/150‑1 (MQ=255)
                                                                                                                                      tGAAAGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTggg         >  1:783090/1‑151 (MQ=255)
                                                                                                                                        aaaGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt        >  2:758584/1‑151 (MQ=255)
                                                                                                                                        aaaGTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgt        <  2:542080/151‑1 (MQ=255)
                                                                                                                                           gTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgt     <  1:753993/151‑1 (MQ=255)
                                                                                                                                           gTGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgt     <  2:458544/151‑1 (MQ=255)
                                                                                                                                           gTGAAAGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgt     <  1:210679/151‑1 (MQ=255)
                                                                                                                                            tGAAGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgtt    >  1:897236/1‑151 (MQ=255)
                                                                                                                                              aaGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGgttgtt    <  2:476316/149‑1 (MQ=255)
                                                                                                                                               aGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGa  >  2:155969/1‑150 (MQ=255)
                                                                                                                                               aGGCAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATAGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGa  >  2:573567/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACCTAATGTCTGCAGAGTTCATGGTAATATGTAAAAAAATATTGTTTAGGAACTGTGTCATTGTATCTTTGTTTGTTTTTACATACAACACATGGGCGCAGTGTAATAATAATATTAAAATAATGCGCAAGTATGAAAGTGAAGGTAAATATACCGTTAGAAATTTGGTTAAAAATAAAGCTATAGCATTGGAATTAGCTGAGATATATGTTAAGAATCGTTATGGACAGGATGCCGCAGAAGAAGAAAAACCATACGAAATCACTGAGTTAACAACAAGTTGGGTTGTTGA  >  NC_000913/3623393‑3623685

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