Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 965,520 A→G 5.1% I401V (ATC→GTC)  ycaI → ComEC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913965,5200AG 5.1% 256.0 / 3.1 78I401V (ATC→GTC) ycaIComEC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base G (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACA  >  NC_000913/965374‑965664
                                                                                                                                                  |                                                                                                                                                
cgccgGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTGTc                                                                                                                                                >  3:404034/1‑149 (MQ=255)
    ggTATTGCGTGCAGTTGTTTCCCTCATGCATCTGCAACTGGGAATCACACTTCTGCTTATGCACGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACgg                                                                                                                                            <  1:168474/149‑1 (MQ=255)
        ttGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTcc                                                                                                                                        >  4:488802/1‑149 (MQ=255)
            gTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTg                                                                                                                                    >  3:201437/1‑149 (MQ=255)
            gTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTg                                                                                                                                    >  3:515281/1‑149 (MQ=255)
            gTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTg                                                                                                                                    >  1:590401/1‑149 (MQ=255)
               cAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATc                                                                                                                                 <  2:598651/149‑1 (MQ=255)
               cAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATc                                                                                                                                 <  2:70617/149‑1 (MQ=255)
                           tCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATg                                                                                                                     >  4:248568/1‑149 (MQ=255)
                            catccatcTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATgg                                                                                                                    <  2:435749/149‑1 (MQ=255)
                                    tGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCAt                                                                                                            <  2:556602/149‑1 (MQ=255)
                                       aaCTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTGTCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTa                                                                                                         >  1:337487/1‑149 (MQ=255)
                                        aCTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTaa                                                                                                        <  4:378399/149‑1 (MQ=255)
                                          tGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGc                                                                                                      >  1:228014/1‑149 (MQ=255)
                                           gggAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCATCGCCGCGATGGTTGTGCATTTAAGCg                                                                                                     >  3:398777/1‑149 (MQ=255)
                                              aaTCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCATGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGc                                                                                                  >  1:459712/1‑149 (MQ=255)
                                               aTCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGCTGGTTGTGCATTTAAGCGGGcc                                                                                                 >  2:26557/1‑149 (MQ=255)
                                                      ttCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATc                                                                                          <  1:441642/149‑1 (MQ=255)
                                                      ttCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATc                                                                                          <  4:515281/149‑1 (MQ=255)
                                                      ttCTGCTTATGCCCATGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATc                                                                                          <  4:335245/149‑1 (MQ=255)
                                                         tGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTg                                                                                       <  4:300127/149‑1 (MQ=255)
                                                           cTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                                                                     <  4:524340/149‑1 (MQ=255)
                                                           cTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                                                                     <  3:485591/149‑1 (MQ=255)
                                                              aTGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCa                                                                                  >  4:418220/1‑149 (MQ=255)
                                                               tGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCaa                                                                                 <  1:469326/149‑1 (MQ=255)
                                                                 cccGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGTGCCGTTAATCCTGGAGCAAgg                                                                               <  4:486687/149‑1 (MQ=255)
                                                                     tGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTa                                                                           <  2:301531/149‑1 (MQ=255)
                                                                      gCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTAt                                                                          <  3:217850/149‑1 (MQ=255)
                                                                       cAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATg                                                                         >  2:362659/1‑149 (MQ=255)
                                                                        aaaTCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATgg                                                                                                                    <  2:180875/105‑1 (MQ=255)
                                                                        aaaTCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATgg                                                                                                                    >  1:180875/1‑105 (MQ=255)
                                                                           tCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGttt                                                                     <  4:333770/149‑1 (MQ=255)
                                                                            cGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTc                                                                    >  1:9173/1‑149 (MQ=255)
                                                                              tCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCtt                                                                  <  2:363730/149‑1 (MQ=255)
                                                                              tCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCtt                                                                  >  2:485543/1‑149 (MQ=255)
                                                                              tCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCtt                                                                  <  1:552076/149‑1 (MQ=255)
                                                                                atatTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCa                                                                                                                                               >  3:571048/1‑70 (MQ=255)
                                                                                atatTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCa                                                                                                                                               <  4:571048/70‑1 (MQ=255)
                                                                                atatTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGc                                                                >  2:149444/1‑149 (MQ=255)
                                                                                  atTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCg                                                              >  3:513069/1‑149 (MQ=255)
                                                                                    ttCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGAc                                                            >  3:349606/1‑149 (MQ=255)
                                                                                         ggCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTc                                                       >  1:158682/1‑149 (MQ=255)
                                                                                           cATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCtt                                                     >  2:226081/1‑149 (MQ=255)
                                                                                            aTTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCttt                                                    >  1:328576/1‑149 (MQ=255)
                                                                                            aTTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCttt                                                    >  1:438677/1‑149 (MQ=255)
                                                                                            aTTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCttt                                                    <  4:201437/149‑1 (MQ=255)
                                                                                            aTTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCttt                                                    >  1:299896/1‑149 (MQ=255)
                                                                                                    gACCTCGTTTATTGCAAATCTATTAGCAACTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTAc                                            <  1:362659/149‑1 (MQ=255)
                                                                                                        tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTGTCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACtttt                                        <  4:404034/149‑1 (MQ=255)
                                                                                                        tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACtttt                                        >  1:424919/1‑149 (MQ=255)
                                                                                                        tCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACtttt                                        >  3:55315/1‑149 (MQ=255)
                                                                                                          gTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCt                                      >  3:433325/1‑149 (MQ=255)
                                                                                                          gTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCt                                      <  3:489011/149‑1 (MQ=255)
                                                                                                              aTTGCAAATCTATTAGCAATTCCCTTGGTGACATTTGTCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTgggg                                  <  2:337487/149‑1 (MQ=255)
                                                                                                              aTTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTgggg                                  <  4:184618/149‑1 (MQ=255)
                                                                                                              aTTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTgggg                                  <  1:261235/149‑1 (MQ=255)
                                                                                                              aTTGCAAATCTATTAGCAAATCCCTTGGTGACATTTATCACGGTTCCGTTGATCATCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTgggg                                  <  4:398777/149‑1 (MQ=255)
                                                                                                                tGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGt                                <  4:470399/149‑1 (MQ=255)
                                                                                                                   aaaTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTaa                             >  4:104588/1‑149 (MQ=255)
                                                                                                                     aTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                                                                     >  4:477238/1‑91 (MQ=255)
                                                                                                                     aTCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGa                                                                                     <  3:477238/91‑1 (MQ=255)
                                                                                                                      tCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAaga                          <  2:328999/149‑1 (MQ=255)
                                                                                                                      tCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAaga                          >  1:328999/1‑149 (MQ=255)
                                                                                                                          ttAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGttt                      >  3:146262/1‑149 (MQ=255)
                                                                                                                            aGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGc                    >  1:256435/1‑149 (MQ=255)
                                                                                                                             gCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGcc                   <  3:104588/149‑1 (MQ=255)
                                                                                                                              cAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCg                  <  4:146262/149‑1 (MQ=255)
                                                                                                                                 ttCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGaa               <  3:510806/149‑1 (MQ=255)
                                                                                                                                 ttCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGaa               >  2:476044/1‑149 (MQ=255)
                                                                                                                                    ccTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAggg            >  1:278622/1‑149 (MQ=255)
                                                                                                                                     cTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGt           <  3:488802/149‑1 (MQ=255)
                                                                                                                                      ttGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTg          >  4:429210/1‑149 (MQ=255)
                                                                                                                                          tGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATc      >  4:536704/1‑149 (MQ=255)
                                                                                                                                          tGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATc      >  2:553843/1‑149 (MQ=255)
                                                                                                                                           gACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCa     >  2:37515/1‑149 (MQ=255)
                                                                                                                                              aTTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACa  <  2:278622/149‑1 (MQ=255)
                                                                                                                                              aTTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACa  >  2:383232/1‑149 (MQ=255)
                                                                                                                                              aTTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACa  >  3:130778/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGCCGGTATTGCGTGCAGTTGTTTCCCTCATCCATCTGCAACTGGGAATCACACTTCTGCTTATGCCCGTGCAAATCGTCATATTTCATGGCATTAGTCTGACCTCGTTTATTGCAAATCTATTAGCAATTCCCTTGGTGACATTTATCACGGTTCCGTTGATCCTCGCCGCGATGGTTGTGCATTTAAGCGGGCCGTTAATCCTGGAGCAAGGGTTATGGTTTCTTGCCGACCGGTCTTTGGCTTTACTTTTCTGGGGGTTAAAGAGTTTGCCGGAAGGGTGGATCAACA  >  NC_000913/965374‑965664

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