Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 999,628 C→G 5.9% G165A (GGC→GCC)  PP_0861 ← outer membrane ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947999,6280CG 5.9% 323.6 / 10.4 103G165A (GGC→GCC) PP_0861outer membrane ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (48/48);  new base G (3/3);  total (51/52)
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.31e-01

Only 100 of 109 total aligned reads displayed.

TCAGGCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCG  >  NC_002947/999479‑999770
                                                                                                                                                     |                                                                                                                                              
tCAGGCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGgccgcgc                                                                                                                                                <  2:285144/150‑1 (MQ=255)
   ggCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGgccgcggcgc                                                                                                                                             >  2:366696/1‑150 (MQ=255)
           gCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGGCGCGGccggcg                                                                                                                                     <  1:366696/150‑5 (MQ=255)
           gCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCg                                                                                                                                     >  2:2076442/1‑150 (MQ=255)
              aCGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGGCGCGGccggcggtc                                                                                                                                  <  2:1606013/150‑8 (MQ=255)
                 gTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGTTTCCGGACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGcc                                                                                                                               >  2:1230272/1‑150 (MQ=255)
                 gTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGACGGcc                                                                                                                               <  2:1844781/150‑1 (MQ=255)
                     cGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGt                                                                                                                           <  1:2076442/150‑1 (MQ=255)
                      gCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGtt                                                                                                                          >  1:1682873/1‑150 (MQ=255)
                       cTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTc                                                                                                                         >  2:1911644/1‑150 (MQ=255)
                       cTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCAAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTc                                                                                                                         >  1:1167905/1‑150 (MQ=255)
                        tGAACTGGTAGTTGCCATCGAAGGTCAAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCg                                                                                                                        >  1:315726/1‑150 (MQ=255)
                          aaCTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGc                                                                                                                      >  1:1271481/1‑150 (MQ=255)
                               gTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCACtt                                                                                                                 >  1:1169948/1‑150 (MQ=255)
                                tAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCACCAGGTCCAGGGTCCGGCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                <  1:1686442/150‑1 (MQ=255)
                                 aGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTgg                                                                                                               <  1:1381448/150‑1 (MQ=255)
                                  gTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGTTACCGATGGCCGAGTTCGGCCCCTTGGc                                                                                                              <  1:1230272/150‑1 (MQ=255)
                                  gTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGc                                                                                                              <  2:832885/150‑1 (MQ=255)
                                   ttGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGcc                                                                                                             >  2:1389364/1‑150 (MQ=255)
                                      ccATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACt                                                                                                          >  2:290271/1‑150 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGccgcgccgcg                                                                                                                                            <  1:1265644/115‑1 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGccgcgccgcg                                                                                                                                            >  2:1265644/1‑115 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCCCCGATGGCCGAGTTCGGCCCCTTGGCCACtt                                                                                                         >  2:862330/1‑150 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGc                                                                                                              >  1:1282459/1‑145 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGc                                                                                                              <  2:1282459/145‑1 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACtt                                                                                                         <  2:2193823/150‑1 (MQ=255)
                                       cATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACtt                                                                                                         >  1:16628/1‑150 (MQ=255)
                                         tCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCa                                                                                                       >  2:838283/1‑150 (MQ=255)
                                          cGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCaa                                                                                                      <  1:838283/150‑1 (MQ=255)
                                           gAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCCCTTCAAc                                                                                                     >  1:976151/1‑150 (MQ=255)
                                           gAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAAc                                                                                                     >  1:1611076/1‑150 (MQ=255)
                                               gtgtAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGAt                                                                                                 <  1:1627186/150‑1 (MQ=255)
                                                tgtAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGAtt                                                                                                <  2:1030104/150‑1 (MQ=255)
                                                  tAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCg                                                                                              <  1:814276/150‑1 (MQ=255)
                                                    gcgcTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGACCACTTCAACAGATTCGAt                                                                                            <  1:761991/150‑1 (MQ=255)
                                                              tGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGAtt                                                                                  >  2:1201191/1‑150 (MQ=255)
                                                              tGGTCGCTGCCCCAGGTCAAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGAtt                                                                                  <  2:1167905/150‑1 (MQ=255)
                                                                gTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTc                                                                                >  1:1999439/1‑150 (MQ=255)
                                                                 tCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTcg                                                                               <  1:1182145/150‑1 (MQ=255)
                                                                  cgatGCGCCACGTCGAGGCGCCGCCCAGCGGGTTCCCCAGATGCGCGGGCTTGCTTACCAGGTTGATGGGACCGCCGGCGGCGCCGCGGCCTCCGATGGCCGCGTTAGGCCCCTTGGCCACTTCAACAGGTTCGATGGCGAAGATTTcgc                                                                              <  2:372566/147‑1 (MQ=255)
                                                                   gCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTcgcg                                                                             >  1:1408527/1‑150 (MQ=255)
                                                                         ccAGGTCCAGGCGCCGTCCAGCGAGTTGCCCATATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTAAGCGCTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTg                                                                       >  1:23064/1‑150 (MQ=255)
                                                                         ccAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTg                                                                       >  1:1029822/1‑150 (MQ=255)
                                                                         ccAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTAACAAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTg                                                                       >  1:1812011/1‑150 (MQ=255)
                                                                              tCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTgg                                                                                                               >  1:370597/1‑105 (MQ=255)
                                                                              tCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTgg                                                                                                               <  2:370597/105‑1 (MQ=255)
                                                                              tCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGgcgc                                                                  <  1:1546329/150‑1 (MQ=255)
                                                                                   gcgcCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAAGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAgtg                                                             <  1:1911644/150‑1 (MQ=255)
                                                                                    cgcCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAgtgt                                                            <  1:2015746/150‑1 (MQ=255)
                                                                                      ccGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                >  2:534346/1‑96 (MQ=255)
                                                                                       cGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcac                                                         <  2:834739/150‑1 (MQ=255)
                                                                                          ccAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcacgca                                                      >  2:578751/1‑150 (MQ=255)
                                                                                           cAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTcacgcac                                                     <  1:1389364/150‑1 (MQ=255)
                                                                                            aGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATgg                                                                                                                                 <  1:2000462/73‑1 (MQ=255)
                                                                                            aGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATgg                                                                                                                                 >  2:2000462/1‑73 (MQ=255)
                                                                                             gCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                <  1:1911641/89‑1 (MQ=255)
                                                                                             gCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTg                                                                                                                >  2:1911641/1‑89 (MQ=255)
                                                                                                 gTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTc                                               <  2:1408527/150‑1 (MQ=255)
                                                                                                   tGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCa                                             <  2:1999439/150‑1 (MQ=255)
                                                                                                    gCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAg                                            <  2:1877574/150‑1 (MQ=255)
                                                                                                      ccAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTc                                                                                                        >  2:1223265/1‑88 (MQ=255)
                                                                                                      ccAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTc                                                                                                        <  1:1223265/88‑1 (MQ=255)
                                                                                                       cAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTa                                         >  2:2034289/1‑150 (MQ=255)
                                                                                                           tGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTa                                     <  2:1169948/150‑1 (MQ=255)
                                                                                                                    ttgcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTgg                            >  2:1557670/1‑150 (MQ=255)
                                                                                                                    ttgcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTgg                            >  1:900852/1‑150 (MQ=255)
                                                                                                                     tgcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTggg                           <  2:257137/150‑1 (MQ=255)
                                                                                                                      gcttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGc                          >  1:318358/1‑150 (MQ=255)
                                                                                                                        ttACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGt                        >  1:2326404/1‑150 (MQ=255)
                                                                                                                           ccAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGa                     >  2:1273237/1‑150 (MQ=255)
                                                                                                                            cAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGAtt                                                                                  >  2:1332269/1‑88 (MQ=255)
                                                                                                                            cAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGaa                    <  1:1932620/150‑1 (MQ=255)
                                                                                                                            cAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGaa                    >  2:268633/1‑150 (MQ=255)
                                                                                                                            cAGGTTGATGGTACCGCCAGCCACGAAGACCACCTCGATGGCCTTGTTCGGGCTTTTTGCCACTTAAAAAACGTTGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCCaa                    <  1:521172/150‑1 (MQ=255)
                                                                                                                            cAGGTTGATGGTACCGCAGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGAtt                                                                                  <  1:1332269/88‑1 (MQ=255)
                                                                                                                              ggTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTAGCCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGc                  <  1:1227704/150‑1 (MQ=255)
                                                                                                                                ttGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGTCGGCGAATCcg                >  2:59064/1‑150 (MQ=255)
                                                                                                                                 tGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCTCTTGACAGCATGAATAAAGCCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCcgc               <  1:862330/150‑1 (MQ=255)
                                                                                                                                 tGATGGTACCGCCGGCCGCGCCGCGGCAACGGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCCATGGCGAAGATTTCGCGGGTTTGGGCGCCAGNGTCACCNACGCCGTCCAGGNAAGTATCGCCCTGGGCGTCGAAGCcgc               <  1:2284976/150‑1 (MQ=255)
                                                                                                                                   aTGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCgg             <  1:290271/150‑1 (MQ=255)
                                                                                                                                     ggTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGTCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGAt           >  2:2081176/1‑150 (MQ=255)
                                                                                                                                       tACCGCCGGCCGCGCCGCGGCCACCGATGGCCTAGTTCGGcccc                                                                                                                   <  2:573503/44‑1 (MQ=255)
                                                                                                                                       tACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGcccc                                                                                                                   >  1:573503/1‑44 (MQ=255)
                                                                                                                                         tcggtggccgCGGCGCGGCCGGCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTc                       <  1:339577/129‑1 (MQ=255)
                                                                                                                                         tcggtggccgCGGCGCGGCCGGCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTc                       >  2:339577/6‑134 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAg                                                                                     <  1:493492/72‑1 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAg                                                                                     >  2:493492/1‑72 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGt                                          <  2:818188/115‑1 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGt                                          >  1:818188/1‑115 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGc                                 <  2:289075/124‑1 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAAc       <  2:1682873/150‑1 (MQ=255)
                                                                                                                                         ccgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCAGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGc                                 >  1:289075/1‑124 (MQ=255)
                                                                                                                                          cgccgGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACg      >  1:1216801/1‑150 (MQ=255)
                                                                                                                                           gccggccgCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACgg     >  1:1269340/1‑150 (MQ=255)
                                                                                                                                            gtggccgCGGCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACggg    >  1:2258112/3‑150 (MQ=255)
                                                                                                                                             tggccgCGGCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGc   >  2:47872/2‑150 (MQ=255)
                                                                                                                                             cggccgCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGc   >  2:343528/1‑150 (MQ=255)
                                                                                                                                              ggccgCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCg  <  1:2081176/150‑1 (MQ=255)
                                                                                                                                               gccgcgccgcGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTggg                           <  2:2220029/124‑1 (MQ=255)
                                                                                                                                               gccgcgccgcGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTggg                           >  1:2220029/1‑124 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TCAGGCGACCGGCAACGGTATCGCTGAACTGGTAGTTGCCATCGAAGGTGTAGCGCTGGGTCTGGTCGCTGCCCCAGGTCCAGGCGCCGTCCAGCGAGTTGCCCAGATGCGCGCGCTTGCTTACCAGGTTGATGGTACCGCCGGCCGCGCCGCGGCCACCGATGGCCGAGTTCGGCCCCTTGGCCACTTCAACAGATTCGATGGCGAAGATTTCGCGGGTTTGGGCGCCAGTGTCACGCACGCCGTCCAGGTAAGTATCGCCCTGGGCGTCGAAGCCGCGGATGAACGGGCG  >  NC_002947/999479‑999770

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: