Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 979,121 A→T 5.7% I190F (ATC→TTC)  iscS‑I → cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947979,1210AT 5.7% 274.8 / 6.4 88I190F (ATC→TTC) iscS‑Icysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (40/43);  new base T (2/3);  total (42/46)
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.04e-01

TGCTGCGCGACGACACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAGCGCGGC  >  NC_002947/978972‑979267
                                                                                                                                                     |                                                                                                                                                  
tgctgcGCGACGACACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGaaa                                                                                                                                                    <  2:2333369/150‑1 (MQ=255)
  ctgcGCGACGACACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATc                                                                                                                                                  <  2:455162/150‑1 (MQ=255)
      gcgACGACACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGAcc                                                                                                                                              >  2:1554223/1‑150 (MQ=255)
       cgacgaCACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCt                                                                                                                                             <  1:600992/150‑1 (MQ=255)
        gacgacACCATCCTGGCCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGATTTCGACCtt                                                                                                                                            <  1:831184/150‑2 (MQ=255)
           gacACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCaa                                                                                                                                         >  2:1368134/1‑150 (MQ=255)
              acCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAg                                                                                                                                      <  1:2361970/150‑1 (MQ=255)
              acCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAg                                                                                                                                      <  1:2358553/150‑1 (MQ=255)
              acCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAg                                                                                                                                      >  1:966352/1‑150 (MQ=255)
                cATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCt                                                                                                                                    >  1:1528515/1‑150 (MQ=255)
                         ctcGCTGATGCATGTGAACAACGACGTCGGCTCGATCACCGACATCGCCGCCACCGGTGAACTGCCCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGa                                                                                                                           <  2:491749/150‑1 (MQ=255)
                          tcGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGAc                                                                                                                          <  2:1528515/150‑1 (MQ=255)
                            gCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCTGTCTACGGGTTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCt                                                                                                                        <  1:1860784/150‑1 (MQ=255)
                            gCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCt                                                                                                                        >  1:1999258/1‑150 (MQ=255)
                              tGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGa                                                                                                                      <  2:1645987/150‑1 (MQ=255)
                                       tgAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTtct                                                                                                             >  2:716990/1‑150 (MQ=255)
                                        gAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTACCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTtctc                                                                                                            >  2:1921900/1‑150 (MQ=255)
                                           caacGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGgc                                                                                                         >  2:1185278/1‑150 (MQ=255)
                                           caacGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACAGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGAc                                                                                                                                               >  1:899219/1‑112 (MQ=255)
                                           caacGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACAGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGAc                                                                                                                                               <  2:899219/112‑1 (MQ=255)
                                               gAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGcac                                                                                                     >  1:2150314/1‑150 (MQ=255)
                                               gAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGcac                                                                                                     >  2:2122019/1‑150 (MQ=255)
                                                 aGTCGGCTCGATCAACGACATCGCCGCCATCGGTTAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACaa                                                                                                   >  1:1050869/1‑150 (MQ=255)
                                                 aGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACaa                                                                                                   <  2:2106870/150‑1 (MQ=255)
                                                   tCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAgg                                                                                                 <  2:1195899/150‑1 (MQ=255)
                                                     ggCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTc                                                                                               >  1:742348/1‑150 (MQ=255)
                                                      gCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCt                                                                                              >  1:949329/1‑150 (MQ=255)
                                                          gATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGAc                                                                                                                                               >  2:1778711/1‑97 (MQ=255)
                                                          gATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGAc                                                                                                                                               <  1:1778711/97‑1 (MQ=255)
                                                          gATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATgg                                                                                          >  1:2448235/1‑150 (MQ=255)
                                                            tCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGcc                                                                                        <  2:742348/150‑1 (MQ=255)
                                                             cAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGccc                                                                                       <  2:949329/150‑1 (MQ=255)
                                                              aaCGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCg                                                                                      >  2:149488/1‑150 (MQ=255)
                                                                cGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTa                                                                                             <  2:1267680/141‑1 (MQ=255)
                                                                cGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTa                                                                                             >  1:1267680/1‑141 (MQ=255)
                                                                 gACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCACGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGaaa                                                                                   <  2:1901040/150‑1 (MQ=255)
                                                                 gACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGaaa                                                                                   >  1:2012415/1‑150 (MQ=255)
                                                                      cgccgcCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCAt                                                                              <  1:731635/150‑1 (MQ=255)
                                                                          gccATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGgc                                                                          <  2:559882/150‑1 (MQ=255)
                                                                            cATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGgcgc                                                                        >  2:1913740/1‑150 (MQ=255)
                                                                                ggTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTg                                                                    <  1:2385911/150‑1 (MQ=255)
                                                                                   gAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTAc                                                                 <  1:1368134/150‑1 (MQ=255)
                                                                                   gAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTAc                                                                 <  2:2258171/150‑1 (MQ=255)
                                                                                   gAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTAc                                                                 <  1:1428233/150‑1 (MQ=255)
                                                                                       tGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCa                                                             <  2:2012415/150‑1 (MQ=255)
                                                                                        gACCCGTTCGCGCGGTGTACTGGTCCATGTGGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAg                                                            >  1:565678/1‑150 (MQ=255)
                                                                                         aCCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAgg                                                                                                 >  1:2066807/1‑112 (MQ=255)
                                                                                         aCCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAgg                                                                                                 <  2:2066807/112‑1 (MQ=255)
                                                                                         aCCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGc                                                           <  1:2374047/150‑1 (MQ=255)
                                                                                          cccGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGcc                                                          <  2:2362101/150‑1 (MQ=255)
                                                                                           ccGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAg                                                                                                  >  2:1646693/1‑109 (MQ=255)
                                                                                           ccGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAg                                                                                                  <  1:1646693/109‑1 (MQ=255)
                                                                                              ttCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCaa                                                      >  2:2432755/1‑150 (MQ=255)
                                                                                              ttCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCaa                                                      <  1:2443754/150‑1 (MQ=255)
                                                                                                   gcgGTGTGCTGTTCCATGCCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACATAATGTCGTTCTCGGCGCACAGGGTCTATGGCCCGAAAGGCATCGGTGCGCTGTACGTCAGCCGCAAGCcgc                                                 <  2:285494/150‑1 (MQ=255)
                                                                                                    cgGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCcgcg                                                <  1:1041691/150‑1 (MQ=255)
                                                                                                     ggTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCgt                                               <  1:1185278/150‑1 (MQ=255)
                                                                                                     ggTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCgt                                               <  1:2432755/150‑1 (MQ=255)
                                                                                                         tACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTAc                                           >  1:465931/1‑150 (MQ=255)
                                                                                                            tGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCt                                                                     <  2:960523/121‑1 (MQ=255)
                                                                                                            tGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCt                                                                     >  1:960523/1‑121 (MQ=255)
                                                                                                              ttCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCtt                                      >  1:106618/1‑150 (MQ=255)
                                                                                                                ccATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGgc                                                                          <  2:1850207/112‑1 (MQ=255)
                                                                                                                ccATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGgc                                                                          >  1:1850207/1‑112 (MQ=255)
                                                                                                                   tGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGc                                 <  1:390475/150‑1 (MQ=255)
                                                                                                                   tGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGc                                 >  2:168576/1‑150 (MQ=255)
                                                                                                                    gTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGcc                                <  1:1445455/150‑1 (MQ=255)
                                                                                                                    gTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGcc                                <  2:1871174/150‑1 (MQ=255)
                                                                                                                     tCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCca                               <  2:1999258/150‑1 (MQ=255)
                                                                                                                       gACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCcatc                             >  1:65583/1‑150 (MQ=255)
                                                                                                                        aCGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCcatca                            >  1:790221/1‑150 (MQ=255)
                                                                                                                              cccAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCgt                                               <  1:1411147/125‑1 (MQ=255)
                                                                                                                              cccAGTCGGCCGGCAAGGTCGAAATCGACCCGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCgt                                               >  2:1411147/1‑125 (MQ=255)
                                                                                                                               ccAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCAcgg                     >  2:2132002/1‑150 (MQ=255)
                                                                                                                               ccAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCAcgg                     >  2:2260715/1‑150 (MQ=255)
                                                                                                                               ccAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCAcgg                     >  1:225175/1‑150 (MQ=255)
                                                                                                                                  gTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTc                                                                                                                  >  1:322752/1‑54 (MQ=255)
                                                                                                                                  gTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTc                                                                                                                  <  2:322752/54‑1 (MQ=255)
                                                                                                                                     ggccggcAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCAcggcggcgg               >  2:244909/1‑150 (MQ=255)
                                                                                                                                           gcaGGTCGATTTCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGaa                                   <  2:561767/122‑1 (MQ=255)
                                                                                                                                           gcaGGTCGATTTCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGaa                                   >  1:561767/3‑124 (MQ=255)
                                                                                                                                           gcaGGTCGATTTCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGaa                                   <  2:558876/122‑1 (MQ=255)
                                                                                                                                           gcaGGTCGATTTCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGaa                                   >  1:558876/3‑124 (MQ=255)
                                                                                                                                           cAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGt                                                               <  1:210117/96‑1 (MQ=255)
                                                                                                                                           cAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGt                                                               >  2:210117/1‑96 (MQ=255)
                                                                                                                                              ggTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCAGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAgcg      <  1:2390387/150‑1 (MQ=255)
                                                                                                                                                tCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAGcgcg    >  1:1945038/1‑150 (MQ=255)
                                                                                                                                                  gAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAGCGCGGc  >  2:42245/1‑150 (MQ=255)
                                                                                                                                                  gAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAGCGCGGc  <  2:106618/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCTGCGCGACGACACCATCCTGGTCTCGCTGATGCATGTGAACAACGAAGTCGGCTCGATCAACGACATCGCCGCCATCGGTGAACTGACCCGTTCGCGCGGTGTACTGTTCCATGTCGACGCGGCCCAGTCGGCCGGCAAGGTCGAAATCGACCTGCAAAAGCTGAAAGTCGACCTGATGTCGTTCTCGGCGCACAAGGTCTATGGCCCGAAAGGCATCGGCGCGCTGTACGTCAGCCGCAAGCCGCGTGTACGTCTTGAAGCCATCATTCACGGCGGCGGTCATGAGCGCGGC  >  NC_002947/978972‑979267

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