Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,845,838 G→A 10.7% W29* (TGG→TGA)  pfeS‑II → histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,845,8380GA10.7% 181.9 / 10.0 73W29* (TGG→TGA) pfeS‑IIhistidine kinase
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (4/4);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAGCGCGA  >  NC_002947/1845689‑1845984
                                                                                                                                                     |                                                                                                                                                  
cGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGCGCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTgg                                                                                                                                                    <  1:159671/150‑1 (MQ=255)
   cAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                                 <  1:988372/150‑1 (MQ=255)
     gATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCg                                                                                                                                               <  1:1822959/150‑1 (MQ=255)
       tCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGCCgat                                                                                                                                             <  1:595597/150‑3 (MQ=255)
           cacaGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTgg                                                                                                                                                    <  2:522879/139‑1 (MQ=255)
           cACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTgg                                                                                                                                                    >  1:522879/1‑139 (MQ=255)
           cACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                                         <  1:2229878/150‑1 (MQ=255)
                tgtgGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGaaa                                                                                                                                    <  1:875377/150‑1 (MQ=255)
                     ggCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGcc                                                                                                                                                <  2:309360/133‑1 (MQ=255)
                     ggCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGcc                                                                                                                                                >  1:309360/1‑133 (MQ=255)
                      gCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGc                                                                                                                                                 >  1:1524282/1‑131 (MQ=255)
                      gCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGc                                                                                                                                                 <  2:1524282/131‑1 (MQ=255)
                      gCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGa                                                                                                                              <  1:728827/150‑1 (MQ=255)
                        aaGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAAc                                                                                                                                   <  1:1881820/143‑1 (MQ=255)
                        aaGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGACTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAAc                                                                                                                                   >  2:1881820/1‑143 (MQ=255)
                           ggTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGcc                                                                                                                         <  2:1145027/150‑1 (MQ=255)
                              tACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCATCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTc                                                                                                                      <  2:1220116/150‑1 (MQ=255)
                                 gTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCctc                                                                                                                   <  2:2355523/150‑1 (MQ=255)
                                          gCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAACCCCAGAATGCCTTCCTCTCGGAGCCg                                                                                                          <  2:826989/150‑1 (MQ=255)
                                              gCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCAc                                                                                                      >  2:894449/1‑150 (MQ=255)
                                                  ggCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCa                                                                                                  <  1:1081693/150‑1 (MQ=255)
                                                   gCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAg                                                                                                 >  1:643766/1‑150 (MQ=255)
                                                     ggAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCCGAc                                                                                               <  2:643766/150‑1 (MQ=255)
                                                           aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCg                                                                                         <  1:1705035/150‑1 (MQ=255)
                                                           aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCg                                                                                         >  2:1979339/1‑150 (MQ=255)
                                                           aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCg                                                                                         >  2:1313409/1‑150 (MQ=255)
                                                           aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTAGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATTAAGACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCg                                                                                         >  2:273918/1‑150 (MQ=255)
                                                             ccATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCggg                                                                                       <  2:777530/150‑1 (MQ=255)
                                                              cATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCgggg                                                                                      >  1:767746/1‑150 (MQ=255)
                                                              cATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCgggg                                                                                      >  1:320047/1‑150 (MQ=255)
                                                                tGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAAt                                                                                                                            <  2:1629307/110‑1 (MQ=255)
                                                                tGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAAt                                                                                                                            >  1:1629307/1‑110 (MQ=255)
                                                                  cTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                                  <  1:210830/150‑1 (MQ=255)
                                                                    cGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgc                                                                                <  1:61943/150‑1 (MQ=255)
                                                                        cGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgcc                                                                            >  2:1492903/1‑150 (MQ=255)
                                                                        cGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTATGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTGGGAGCCGGCACGCCAGTCCTTGCGGGGTTATgccgcc                                                                            >  1:1691871/1‑150 (MQ=255)
                                                                           cATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGaa                                                                         <  2:767746/150‑1 (MQ=255)
                                                                             ttCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGc                                                                       >  1:442501/1‑150 (MQ=255)
                                                                                    ttCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAg                                                                >  2:1689457/1‑150 (MQ=255)
                                                                                    ctCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTACTGCCAATAGGGCGCAGTTGTACCTGTTGCCGTCACATGGAAACCCAGACTGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTCCGGGGTTATGCCGCCGCAGCTGAACAg                                                                >  1:1855708/2‑150 (MQ=255)
                                                                                               ggCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGcc                                                                                                                         <  1:426326/82‑1 (MQ=255)
                                                                                               ggCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGcc                                                                                                                         >  2:426326/1‑82 (MQ=255)
                                                                                                gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                                    >  2:2196471/1‑150 (MQ=255)
                                                                                                  cATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTgg                                                                                                                                                    >  2:2138910/1‑52 (MQ=255)
                                                                                                  cATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTgg                                                                                                                                                    <  1:2138910/52‑1 (MQ=255)
                                                                                                      ttGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcggg                                              >  2:300688/1‑150 (MQ=255)
                                                                                                        gCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggcgggcg                                            <  2:320047/150‑1 (MQ=255)
                                                                                                           ggtggGCTTCTGCCTGCTGATTATCGGCCTCAGTTATACCTGGGGCCCGCACATGGAAACCGCGATTGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCGCAACAGGCCTGGCGCCGTGGCgggcgggc                                         <  2:654497/150‑1 (MQ=255)
                                                                                                            gtggGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCa                                        >  1:1547190/1‑150 (MQ=255)
                                                                                                               ggCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGATGCTTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAggg                                     <  1:1979339/150‑1 (MQ=255)
                                                                                                                gCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGc                                    >  1:1700274/1‑150 (MQ=255)
                                                                                                                  ttCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGc                                    >  1:382039/1‑148 (MQ=255)
                                                                                                                  ttCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGc                                    <  2:382039/148‑1 (MQ=255)
                                                                                                                  ttCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTg                                  <  1:1070581/150‑1 (MQ=255)
                                                                                                                  ttCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTg                                  >  2:1734004/1‑150 (MQ=255)
                                                                                                                   tCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTgg                                 <  1:300688/150‑1 (MQ=255)
                                                                                                                      gcctgcTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCg                                                                                                          >  2:1776714/1‑74 (MQ=255)
                                                                                                                      gcctgcTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCg                                                                                                          <  1:1776714/74‑1 (MQ=255)
                                                                                                                         tgctgATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGACAGCAAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgt                           >  1:2097096/1‑150 (MQ=255)
                                                                                                                          gctgATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAgtg                          >  2:2208172/1‑150 (MQ=255)
                                                                                                                              atgatCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCa                                                                                                                                <  2:1782080/44‑1 (MQ=255)
                                                                                                                              atgatCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCa                                                                                                                                >  1:1782080/1‑44 (MQ=255)
                                                                                                                                  tCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGATGTTGCGTGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcgg                  >  1:542358/1‑150 (MQ=255)
                                                                                                                                   cGGCCCCAGGTATAACTGAGGCCGGCACCGCGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTggcggc                 >  1:1268968/1‑150 (MQ=255)
                                                                                                                                    ggCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACAAGTGGGGGGCGGcc                >  2:460568/1‑150 (MQ=255)
                                                                                                                                     gCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCa               >  2:422634/1‑150 (MQ=255)
                                                                                                                                      ccTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAt              <  2:1247246/150‑1 (MQ=255)
                                                                                                                                       cTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATg             >  1:346020/1‑150 (MQ=255)
                                                                                                                                       cTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATg             >  1:112986/1‑150 (MQ=255)
                                                                                                                                        tCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGTAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGc            <  1:1145417/150‑1 (MQ=255)
                                                                                                                                          aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCAt          <  1:1734004/150‑1 (MQ=255)
                                                                                                                                           gTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATc         <  1:911447/150‑1 (MQ=255)
                                                                                                                                               taCCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAgcg     >  1:1040838/1‑150 (MQ=255)
                                                                                                                                                tccTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAgcgc    <  2:2313901/149‑1 (MQ=255)
                                                                                                                                                  cTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAg                                       >  1:2144052/1‑113 (MQ=255)
                                                                                                                                                  cTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAg                                       <  2:2144052/113‑1 (MQ=255)
                                                                                                                                                  cTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAGCGCGa  >  2:1692660/1‑150 (MQ=255)
                                                                                                                                                  cTGAGGCCGGCCCATGCACACCCAGAATGCCTTCCTCTCGGATCCGGCACGCCACCGCGCTCGCGGTGATGCCGCCGAACCTGAACAAGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGCCCAGCGGGTGGCGGCCATGTATCAGCGCGa  >  2:37005/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGCATCAGCGCGA  >  NC_002947/1845689‑1845984

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