Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,859,571 T→A 6.4% intergenic (‑121/+35) PP_1662 ← / ← PP_1663 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,859,5710TA 6.4% 184.9 / 4.6 62intergenic (‑121/+35)PP_1662/PP_1663hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base A (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCTT  >  NC_002947/1859425‑1859719
                                                                                                                                                  |                                                                                                                                                    
aaGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGAGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTggg                                                                                                                                                   <  1:923402/150‑1 (MQ=255)
    tGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCggg                                                                                                                                               <  2:844705/150‑1 (MQ=255)
       cTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGTGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCt                                                                                                                                            >  2:1212615/1‑150 (MQ=255)
        tGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAGAGAAAAGCCCCGCACTGAGCGGGGCtt                                                                                                                                           >  1:326105/1‑150 (MQ=255)
          cGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCa                                                                                                                                         <  1:365752/150‑1 (MQ=255)
           gACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAg                                                                                                                                        <  2:735370/150‑1 (MQ=255)
                  tgtCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgc                                                                                                                                 <  1:1160941/150‑1 (MQ=255)
                     cGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                                                                                              <  1:718271/150‑1 (MQ=255)
                     cGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTgcagg                                                                                                                              <  1:1170788/150‑1 (MQ=255)
                         cGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGt                                                                                                                          <  2:1074251/150‑1 (MQ=255)
                            cttcttTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgcg                                                                                                                       <  1:1069259/150‑1 (MQ=255)
                                    taGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGt                                                                                                               <  1:849181/150‑1 (MQ=255)
                                         tttGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTTTTTGGGCTTTTTCGTAGCAAATAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGGTCAGTTGccc                                                                                                          >  1:1055063/1‑150 (MQ=255)
                                                gCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCg                                                                                                   >  1:735154/1‑150 (MQ=255)
                                                   gACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAGAAGCCCCGCCCAGTGCGGGGCttttctt                                                                                                                                      >  2:535850/1‑107 (MQ=38)
                                                   gACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAGAAGCCCCGCCCAGTGCGGGGCttttctt                                                                                                                                      <  1:535850/112‑6 (MQ=38)
                                                    aCTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCtt                                                                                               <  2:735154/150‑1 (MQ=255)
                                                      taagtaagTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGc                                                                                             >  1:1216922/1‑150 (MQ=255)
                                                      taagtaagTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGc                                                                                             >  1:1052993/1‑150 (MQ=255)
                                                          taagtaaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATc                                                                                         <  2:22020/150‑1 (MQ=255)
                                                            agtaaAACAGGCATGAAACAAGTTTCATGCGTCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGa                                                                                       >  1:358344/1‑150 (MQ=255)
                                                              taaAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGAcc                                                                                     >  1:127945/1‑150 (MQ=255)
                                                               aaaaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCg                                                                                    >  1:1127227/1‑150 (MQ=255)
                                                                aaaCAGGTATGAAACAAGTTTCATGAAGCTCTAATGGATGCTGCCTGGTTGGGCTTTTTCGCAGCAAAGAAAAGACCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGACCTTCACCGGCTTGCCATCGACCGt                                                                                   <  2:516538/150‑1 (MQ=255)
                                                                 aaCAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                                  >  2:70384/1‑150 (MQ=255)
                                                                   cAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                                                                                            >  2:607932/1‑122 (MQ=255)
                                                                   cAGGCATGAAACAAGTTTCATCAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGc                                                                                                            <  1:607932/122‑1 (MQ=255)
                                                                      gCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCa                                                                                                      <  2:35273/125‑1 (MQ=255)
                                                                      gCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCa                                                                                                      >  1:35273/1‑125 (MQ=255)
                                                                       cATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATcc                                                                            >  1:546163/1‑150 (MQ=255)
                                                                        aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgc                                                                                                                        >  1:843277/1‑105 (MQ=255)
                                                                        aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTgc                                                                                                                        <  2:843277/105‑1 (MQ=255)
                                                                        aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCt                                                                           >  1:514900/1‑150 (MQ=255)
                                                                        aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCt                                                                           <  2:1055063/150‑1 (MQ=255)
                                                                        aTGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCt                                                                           <  2:1216922/150‑1 (MQ=255)
                                                                             aCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                                  >  2:1045357/1‑138 (MQ=255)
                                                                             aCAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTg                                                                                  <  1:1045357/138‑1 (MQ=255)
                                                                                     cATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATcac                                                              <  2:1127227/150‑1 (MQ=255)
                                                                                        gAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                                                                                                   <  1:586045/94‑1 (MQ=255)
                                                                                        gAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGAt                                                                                                                   >  2:586045/1‑94 (MQ=255)
                                                                                        gAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                                                                                                 >  1:1019919/1‑112 (MQ=255)
                                                                                        gAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGc                                                                                                 <  2:1019919/112‑1 (MQ=255)
                                                                                        gAAGCTCGAAGGGATGCTGCCTGAGTGGGCTTTTTCGCAGCAACGATTTTCCCCTCACTTGGCTGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATa                                                           <  2:358344/150‑1 (MQ=255)
                                                                                         aaGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATAc                                                          <  2:546163/150‑1 (MQ=255)
                                                                                           gCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTc                                                        <  1:855562/150‑1 (MQ=255)
                                                                                           gCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTc                                                        >  1:582921/1‑150 (MQ=255)
                                                                                                    gATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCATGCAGCACGAACACATACTCCTTTCCAtt                                               >  2:713553/1‑149 (MQ=255)
                                                                                                          gcCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGAAGTACGTGATTTGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGtctc                                         <  2:1067785/150‑1 (MQ=255)
                                                                                                                   ggCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGgcg                                >  2:389994/1‑150 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGc                                                                                             <  1:908730/86‑1 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGc                                                                                             >  2:908730/1‑86 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATca                                                               >  2:957568/1‑116 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATca                                                               <  1:957568/116‑1 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCtt                                                     <  2:1006205/126‑1 (MQ=255)
                                                                                                                      tttttCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCtt                                                     >  1:1006205/1‑126 (MQ=255)
                                                                                                                            gcagcaAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCAc                       >  2:1233466/1‑150 (MQ=255)
                                                                                                                              agcaAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGc                         >  2:1113448/1‑146 (MQ=255)
                                                                                                                              agcaAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGc                         <  1:1113448/146‑1 (MQ=255)
                                                                                                                                     aaaaGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTa              <  2:127945/150‑1 (MQ=255)
                                                                                                                                      gaaGCCCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGtct                                          <  2:1131291/120‑1 (MQ=255)
                                                                                                                                      gaaGCCCCGCCCAGTGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGtct                                          >  1:1131291/2‑121 (MQ=255)
                                                                                                                                                 cTGGGCGGGGCTTCAGGTAGTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCtt  <  1:389994/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAGATGGCTGCGACAGTGTGTCGCACGTCTTCTTTATAGAGTTTGTTAGCAGACTAAGTAAGTAAAACAGGCATGAAACAAGTTTCATGAAGCTCGAAGGGATGCTGCCTGGGTGGGCTTTTTCGCAGCAAAGAAAAGCCCCGCACTGGGCGGGGCTTCAGGTATTGCAGGCAGTGCGCGATCAGTTGCCCTTCACCGGCTTGCCATCGACCGTGCCATCCTGCAGCACGATCACATACTCCTTGCCATCGGTCTCTACCTGGCGCAGTTGCACCAGCAGGTAGTCCCAGTTCTT  >  NC_002947/1859425‑1859719

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