Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,178,758 G→C 14.1% intergenic (+57/+37) guaA → / ← PP_1033 glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,178,7580GC14.1% 207.5 / 26.9 82intergenic (+57/+37)guaA/PP_1033glutamine‑hydrolyzing GMP synthase/sulfatase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (34/33);  new base C (6/5);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTT  >  NC_002947/1178610‑1178902
                                                                                                                                                      |                                                                                                                                                
tCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGa                                                                                                                                                 >  1:785763/1‑150 (MQ=255)
     ggTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCg                                                                                                                                            >  1:954374/1‑150 (MQ=255)
             catcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                    >  2:1244917/1‑150 (MQ=255)
             catcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCt                                                                                                                                    <  2:108495/150‑1 (MQ=255)
              atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                   <  2:776331/150‑1 (MQ=255)
              atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                   >  2:945640/1‑150 (MQ=255)
              atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                   >  2:692580/1‑150 (MQ=255)
              atcaACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                   <  2:377789/150‑1 (MQ=255)
                 aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGtt                                                                                                                                >  1:604872/1‑150 (MQ=255)
                 aaCGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAACCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGtt                                                                                                                                >  1:1619995/1‑150 (MQ=255)
                   cGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGCTACACCAGAAAGCCCGCATGATGACGGGCTTTCTGTTTCt                                                                                                                              >  2:1722086/1‑150 (MQ=255)
                    gAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCtt                                                                                                                             <  1:2058701/150‑1 (MQ=255)
                     aaaTCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTg                                                                                                                            <  1:945640/150‑1 (MQ=255)
                            cGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTgttcgttc                                                                                                                     <  1:2314554/150‑1 (MQ=255)
                             ggTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCa                                                                                                                    >  2:476548/1‑150 (MQ=255)
                               tATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGa                                                                                                                  <  2:104438/150‑1 (MQ=255)
                                   tcGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCAAGGCCCACCAGGACGGATACACCAGAAAGCCCGCCTGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGc                                                                                                              >  1:2222963/1‑150 (MQ=255)
                                      cgcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           <  2:832396/150‑1 (MQ=255)
                                       gcgTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAt                                                                                                          >  2:913257/1‑150 (MQ=255)
                                            aCCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCGGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                     <  2:2040543/150‑1 (MQ=255)
                                            aCCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTAttcttc                                                                                                     <  1:913257/150‑1 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                    >  2:1277927/1‑117 (MQ=255)
                                              cTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCt                                                                                                                                    <  1:1277927/117‑1 (MQ=255)
                                                       gcggAGCAAGGCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTgg                                                                                          <  2:1338910/147‑1 (MQ=255)
                                                           aGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                      <  1:2233272/150‑1 (MQ=255)
                                                              acgccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGa                                                                                                                  >  2:2001632/3‑119 (MQ=255)
                                                              acgccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGa                                                                                                                  <  1:2001632/117‑1 (MQ=255)
                                                                gccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                   <  1:672093/100‑1 (MQ=38)
                                                                gccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCAATCATGCGGGCTTTCTg                                                                                                                                   >  2:672093/1‑100 (MQ=38)
                                                                 ccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                   >  2:584563/1‑99 (MQ=255)
                                                                 ccgccTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTg                                                                                                                                   <  1:584563/99‑1 (MQ=255)
                                                                  cgccTGCCACTATCGAGTGGGAGTGAGCGGGTAGCCAGCCCAGGCTGGGACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCt                                                                               <  2:1155085/150‑1 (MQ=255)
                                                                       gCCACTATCGAGTGGGAGTGATCGGGTCGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTg                                                                          <  2:107280/150‑1 (MQ=255)
                                                                         cACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                    >  2:1520597/1‑122 (MQ=255)
                                                                         cACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATtcttct                                                                                                    <  1:1520597/122‑1 (MQ=255)
                                                                           cTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAg                                                                      <  2:785763/150‑1 (MQ=255)
                                                                               cGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGt                                                                  >  2:295308/1‑150 (MQ=255)
                                                                                  gTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtat                                                               <  2:847070/150‑1 (MQ=255)
                                                                                   tGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                              <  2:954374/150‑1 (MQ=255)
                                                                                     ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCt                                                                                                             >  1:1735909/1‑101 (MQ=255)
                                                                                     ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCt                                                                                                             <  2:1735909/101‑1 (MQ=255)
                                                                                     ggAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAgg                                                            >  1:1717271/1‑150 (MQ=255)
                                                                                     ggAGTGATCGGGTAGCCAGCCCAGGCTGGAACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTgg                                                                                          >  1:1865800/1‑120 (MQ=255)
                                                                                     ggAGTGATCGGGTAGCCAGCCCAGGCTGGAACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTgg                                                                                          <  2:1865800/120‑1 (MQ=255)
                                                                                       aGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTCGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                              <  1:280907/146‑1 (MQ=255)
                                                                                       aGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTCGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTtata                                                              >  2:280907/1‑146 (MQ=255)
                                                                                           aTCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                      >  1:2390251/1‑150 (MQ=255)
                                                                                                   gccagccTGGGCTGGCTACCCGCAGGACGGATACACCAGAAAGCCCGCATGAt                                                                                                                                                >  2:1608597/1‑53 (MQ=21)
                                                                                                   gccagccTGGGCTGGCTACCCGCAGGACGGATACACCAGAAAGCCCGCATGAt                                                                                                                                                <  1:1608597/53‑1 (MQ=21)
                                                                                                       gccCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCa                                          >  1:1126793/1‑150 (MQ=255)
                                                                                                        cccAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCAc                                         <  1:1148533/150‑1 (MQ=255)
                                                                                                          tgggCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           <  1:480429/80‑1 (MQ=255)
                                                                                                          tgggCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTa                                                                                                           >  2:480429/3‑82 (MQ=255)
                                                                                                             gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                    <  1:1244917/150‑1 (MQ=255)
                                                                                                             gCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                    <  2:2222963/150‑1 (MQ=255)
                                                                                                               tGGCACCCGCAAGAAGGAAACACTAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATg                                  >  2:2021259/1‑150 (MQ=255)
                                                                                                                 gCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGaa                                <  2:302991/150‑1 (MQ=255)
                                                                                                                          aGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCAc                        >  2:1539517/1‑150 (MQ=255)
                                                                                                                            gACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTc                                                                                             <  2:859289/79‑1 (MQ=255)
                                                                                                                            gACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTc                                                                                             >  1:859289/1‑79 (MQ=255)
                                                                                                                                 aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                      <  1:664412/77‑1 (MQ=21)
                                                                                                                                 aagaacCAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCaaa                                                                                      >  2:664412/5‑81 (MQ=21)
                                                                                                                                    cacCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgcg              <  2:1619995/150‑1 (MQ=255)
                                                                                                                                    cacCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgcg              >  2:2395661/1‑150 (MQ=255)
                                                                                                                                       cAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                         >  1:484854/1‑104 (MQ=38)
                                                                                                                                       cAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGaaa                                                         <  2:484854/104‑1 (MQ=38)
                                                                                                                                       cAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGa                                                           <  2:939931/102‑1 (MQ=38)
                                                                                                                                       cAGAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGa                                                           >  1:939931/1‑102 (MQ=38)
                                                                                                                                       cAGAAAGCCCGCAATCATGCGAGCTTTCAGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTATGCGTTATAGGAAATAAACCTGCTCCTCACTGCGATGAGCTCGTCACCAGCCCCGCGTTg           >  2:295255/1‑150 (MQ=255)
                                                                                                                                         gAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg         >  1:1614924/1‑150 (MQ=255)
                                                                                                                                         gAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg         >  2:2012491/1‑150 (MQ=255)
                                                                                                                                         gAAAGCCCGCAATCATGCGGGCTTTCTGGTTCTTGTTACTTCAGATTAATTATGGTTCTGGCTTGCTGGCCAAAATCGCGTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTggg         >  1:44843/1‑150 (MQ=255)
                                                                                                                                           aaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                              <  1:228328/111‑1 (MQ=255)
                                                                                                                                           aaGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGctc                                              >  2:228328/1‑111 (MQ=255)
                                                                                                                                              cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                    <  1:1242358/118‑1 (MQ=255)
                                                                                                                                              cccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGa                                    >  2:1242358/1‑118 (MQ=255)
                                                                                                                                               ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgc               >  1:1715079/1‑138 (MQ=255)
                                                                                                                                               ccGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCAcgc               <  2:1715079/138‑1 (MQ=255)
                                                                                                                                                cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt  >  2:495683/1‑150 (MQ=255)
                                                                                                                                                cGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCtt  >  1:606718/1‑150 (MQ=255)
                                                                                                                                                  cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                      <  2:362576/96‑1 (MQ=255)
                                                                                                                                                  cATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAagca                                                      >  1:362576/1‑96 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TCAGCGGTCGTATCATCAACGAAATCGACGGTATCTCGCGCGTTACCTACGACGTGTCGAGCAAGCCGCCTGCCACTATCGAGTGGGAGTGATCGGGTAGCCAGCCCAGGCTGGCACCCGCAGGACGGATACACCAGAAAGCCCGCATGATTGCGGGCTTTCTGGTTCTTGTTCGTTCAGATTGCTTATTCTTCTGGCTTTCTGGCAAATAGCGCTTTCTGGTAGGCGTTATAGGAAAGCAGCCTGCTCCTCACTGCGATGAACTCGTCACCCTCCACGCGTTGGGGGTGCTT  >  NC_002947/1178610‑1178902

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