Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,328,783 (G)5→6 coding (1015/1650 nt) yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,328,7781.G100.0% 241.1 / NA 62P340P (CCT→CCCyfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (39/23);  total (39/23)

CGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGT  >  NC_000913/2328630‑2328921
                                                                                                                                                     |                                                                                                                                               
cGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAgg                                                                                                                                               <  1:720235/151‑1 (MQ=255)
cGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAgg                                                                                                                                               <  1:321608/151‑1 (MQ=255)
     aTCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAggggggc                                                                                                                                          <  2:626925/151‑3 (MQ=255)
       cGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAggggggcgt                                                                                                                                        <  1:619035/151‑5 (MQ=255)
           cgcgCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACt                                                                                                                                      <  2:221206/150‑1 (MQ=255)
             cgcCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACtttt                                                                                                                                   <  1:837024/151‑1 (MQ=255)
             cgcCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACtttt                                                                                                                                   <  1:788638/151‑1 (MQ=255)
              gcCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTc                                                                                                                                  >  1:418396/1‑151 (MQ=255)
               cccccGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCg                                                                                                                                 <  2:997437/151‑1 (MQ=255)
                      ggTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTcgcg                                                                                                                          >  2:568229/1‑151 (MQ=255)
                      ggTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGACTcgcg                                                                                                                          >  2:586475/1‑151 (MQ=255)
                       gTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTcgcg                                                                                                                          >  2:1179137/1‑150 (MQ=255)
                         cGCTGGCTAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGTATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGAt                                                                                                                        <  1:397200/150‑1 (MQ=255)
                         cGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATc                                                                                                                       <  1:932861/151‑1 (MQ=255)
                             ggCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCa                                                                                                                     <  1:825536/149‑1 (MQ=255)
                             ggCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCa                                                                                                                     <  2:547693/149‑1 (MQ=255)
                             ggCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGa                                                                                                                   <  2:24284/151‑1 (MQ=255)
                                gAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAAct                                                                                                                >  1:385051/1‑151 (MQ=255)
                                   gCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTg                                                                                                             <  1:774702/151‑1 (MQ=255)
                                    cAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTgg                                                                                                            >  1:178262/1‑151 (MQ=255)
                                     aGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGc                                                                                                           >  2:1190008/1‑151 (MQ=255)
                                       aaaCACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCa                                                                                                          >  1:1038244/1‑150 (MQ=255)
                                        aaCACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAAc                                                                                                        >  2:682061/1‑151 (MQ=255)
                                          cacGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGt                                                                                                      <  1:583092/151‑1 (MQ=255)
                                            cGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTAc                                                                                                    >  1:260120/1‑151 (MQ=255)
                                                gCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCtt                                                                                                <  1:864330/151‑1 (MQ=255)
                                                  gTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTc                                                                                               >  1:422783/1‑150 (MQ=255)
                                                  gTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTc                                                                                               >  1:924983/1‑150 (MQ=255)
                                                               cGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCg                                                                                 <  2:1099328/151‑1 (MQ=255)
                                                                    aTGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGc                                                                            >  1:625158/1‑151 (MQ=255)
                                                                        tCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAAtt                                                                        >  1:278446/1‑151 (MQ=255)
                                                                               aTCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTTCTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGttt                                                                 <  2:1105018/151‑1 (MQ=255)
                                                                                tCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTg                                                                >  2:745004/1‑151 (MQ=255)
                                                                                   ccTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGccc                                                             >  1:81294/1‑151 (MQ=255)
                                                                                    cTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCt                                                            >  1:998723/1‑151 (MQ=255)
                                                                                    cTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCt                                                            >  1:267125/1‑151 (MQ=255)
                                                                                       cccGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCtttt                                                         <  1:198156/151‑1 (MQ=255)
                                                                                           aTTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCt                                                     >  2:746021/1‑151 (MQ=255)
                                                                                           aTTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCt                                                     <  2:356568/151‑1 (MQ=255)
                                                                                             tGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTc                                                    <  2:338989/150‑1 (MQ=255)
                                                                                              ggCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTc                                                    >  1:696080/1‑149 (MQ=255)
                                                                                               gCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGt                                                 >  1:739129/1‑151 (MQ=255)
                                                                                                cGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCCGACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACCATCGAGATGCGAATGCAATTTCAGTTTGCCCTTTTCCCTCatt                                                 >  1:787461/1‑148 (MQ=255)
                                                                                                   tttGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGTCCTTTTCGCTCAGTAAgg                                             >  1:128581/1‑151 (MQ=255)
                                                                                                   tttGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAgg                                             >  1:1018395/1‑151 (MQ=255)
                                                                                                     tGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGgt                                            <  1:1066360/150‑1 (MQ=255)
                                                                                                      gCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCTTTTTCGCTCAGTAAGgtcg                                          >  2:636542/1‑151 (MQ=255)
                                                                                                        tGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGgtcgtc                                        >  1:475311/1‑151 (MQ=255)
                                                                                                           aaaaaaGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGtctc                                      >  1:1190826/1‑150 (MQ=255)
                                                                                                           aaaaaaGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCa                                     <  1:975366/151‑1 (MQ=255)
                                                                                                              aaaGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCg                                  >  2:956918/1‑151 (MQ=255)
                                                                                                               aaGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGt                                 <  2:1098622/151‑1 (MQ=255)
                                                                                                                    cGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTg                            >  1:740278/1‑151 (MQ=255)
                                                                                                                      aaTCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGAc                          >  1:20240/1‑151 (MQ=255)
                                                                                                                        tCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACa                         >  2:891139/1‑150 (MQ=255)
                                                                                                                             aCGGTCATGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGGCGTCTCACCGTGGCTGACAAACGt                    >  1:1071162/1‑150 (MQ=255)
                                                                                                                                   aTGGCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTgg             >  2:383992/1‑151 (MQ=255)
                                                                                                                                     ggCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGtt           >  1:837765/1‑151 (MQ=255)
                                                                                                                                      gCTTTCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGtt           >  1:599104/1‑150 (MQ=255)
                                                                                                                                        tttCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTAcc        >  1:676811/1‑151 (MQ=255)
                                                                                                                                          tCGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCa       <  2:299004/150‑1 (MQ=255)
                                                                                                                                           cGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAtt     <  2:644390/151‑1 (MQ=255)
                                                                                                                                           cGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAtt     >  2:802931/1‑151 (MQ=255)
                                                                                                                                           cGCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCAtt     >  2:944102/1‑151 (MQ=255)
                                                                                                                                            gCTGCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTa    <  1:888972/151‑1 (MQ=255)
                                                                                                                                               gCCTCAGGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGt  >  1:311237/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CGGTCATCGTTCGCGCCCCCGTGGTCGCTGGCGAAGCAGAAACACGCTGCGTGGCGCTGCTGTCGGGGATGGTCAGGCAATCGCCTTCCCGATTGGCGTTTTGCTGTAAAAAAGTGCGAATCGCCACGGTCATGGCTTTCGCTGCCTCAGGGGGCGTACTTTTCGCCTCGCGATCCAGAACTCTGGCAACGTACTCTTCACGATCGAGATGCGACTGCAATTTCAGTTTGCCCTTTTCGCTCAGTAAGGTCGTCTCACCGTGGCTGACAAACGTAATGTGGTTACCATTAGT  >  NC_000913/2328630‑2328921

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