Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,377,792 A→T 9.5% intergenic (‑70/+19) blc ← / ← ampC outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,377,7920AT 9.5% 245.1 / 14.9 74intergenic (‑70/+19)blc/ampCouter membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base T (4/3);  total (37/37)
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.21e-01

GGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCG  >  NC_000913/4377649‑4377937
                                                                                                                                               |                                                                                                                                                 
ggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCgg                                                                                                                                              <  1:579404/149‑1 (MQ=255)
                   aGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                                                                               <  2:1457771/129‑1 (MQ=255)
                   aGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                                                                               >  1:1457771/1‑129 (MQ=255)
                     gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTa                                                                                                                         >  4:1211732/1‑149 (MQ=255)
                                 aaTGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                          >  4:1043982/1‑119 (MQ=255)
                                    gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTATTGTAGAGCGTTAAGAATCt                                                                                                          >  3:520732/1‑149 (MQ=255)
                                    gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCt                                                                                                          <  4:520732/149‑1 (MQ=255)
                                       gCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                                                                   <  4:305897/105‑1 (MQ=255)
                                       gCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                                                                   >  3:305897/1‑105 (MQ=255)
                                        cTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCa                                                                                                      >  2:397487/1‑149 (MQ=255)
                                            cGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcg                                                                                                  <  1:397487/149‑1 (MQ=255)
                                             gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAATGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                                                                                                                             <  1:806339/105‑1 (MQ=255)
                                             gCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGa                                                                                                                                             >  2:806339/1‑105 (MQ=255)
                                               tGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGgcggcg                                                                                               >  1:212034/1‑149 (MQ=255)
                                                   gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                                                                                                                                       >  1:1285954/1‑105 (MQ=255)
                                                   gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGAt                                                                                                                                       <  2:1285954/105‑1 (MQ=255)
                                                   gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                                                                            <  1:343217/132‑1 (MQ=255)
                                                   gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAAt                                                                                                            >  2:343217/1‑132 (MQ=255)
                                                   gCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGa                                                                                           >  3:251854/1‑149 (MQ=255)
                                                        gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGa                                                                                                              <  2:1142782/125‑1 (MQ=255)
                                                        gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGa                                                                                                              >  1:1142782/1‑125 (MQ=255)
                                                        gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACtctc                                                                                      <  3:1211732/149‑1 (MQ=255)
                                                             ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                          <  4:452142/92‑2 (MQ=255)
                                                             ggAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcct                                                                                                                                          >  3:452142/1‑91 (MQ=255)
                                                                           aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtt                                                                     <  2:442521/147‑1 (MQ=255)
                                                                           aaTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGtt                                                                     >  1:442521/1‑147 (MQ=255)
                                                                            aTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  <  2:212034/149‑1 (MQ=255)
                                                                                 tCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                                                                           <  3:1000589/71‑1 (MQ=255)
                                                                                 tCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAcc                                                                                                                                           >  4:1000589/1‑71 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                                                                                                                                <  4:779356/79‑1 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                                                                                                                                >  3:779356/1‑79 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTa                                                                                                                              <  3:108609/81‑1 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                                                                >  4:16484/1‑143 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgt                                                                <  3:16484/143‑1 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCa                                                          <  4:251854/149‑1 (MQ=255)
                                                                                    ttACTGGTTTTTTTCTAACTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTa                                                                                                                              >  4:108609/1‑81 (MQ=255)
                                                                                       cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                                                                >  3:137088/1‑108 (MQ=255)
                                                                                       cTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAggcggc                                                                                                <  4:137088/108‑1 (MQ=255)
                                                                                          gTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGt                                                                                                                   >  4:455587/1‑86 (MQ=255)
                                                                                          gTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGt                                                                                                                   <  3:455587/86‑1 (MQ=255)
                                                                                                      gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                                                                           >  2:626217/1‑82 (MQ=255)
                                                                                                      gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATc                                                                                                           <  1:626217/82‑1 (MQ=255)
                                                                                                           gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                                                               <  1:723591/105‑1 (MQ=255)
                                                                                                           gATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                                                               >  2:723591/1‑105 (MQ=255)
                                                                                                            aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                                                                              <  3:958189/105‑1 (MQ=255)
                                                                                                            aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAtt                                                                              >  4:958189/1‑105 (MQ=255)
                                                                                                            aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCttt                                    <  4:225730/147‑1 (MQ=255)
                                                                                                            aTGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCttt                                    >  3:225730/1‑147 (MQ=255)
                                                                                                                 aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttt                                                                    >  1:631020/1‑110 (MQ=255)
                                                                                                                 aGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttt                                                                    <  2:631020/110‑1 (MQ=255)
                                                                                                                      ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTCTTTGTTTGCCAGCATCACGATACCCAGCTCttt                                    >  3:1456473/1‑137 (MQ=255)
                                                                                                                      ggAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTCTTTGTTTGCCAGCATCACGATACCCAGCTCttt                                    <  4:1456473/137‑1 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                                                                        >  4:663068/1‑100 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                                                                        <  3:663068/100‑1 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg                                                                 <  1:500357/107‑1 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttg                                                                 >  2:500357/1‑107 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGa                                                 <  1:861112/123‑1 (MQ=255)
                                                                                                                       gAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGa                                                 >  2:861112/1‑123 (MQ=255)
                                                                                                                        aaaaGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                                                                             <  3:860674/94‑1 (MQ=255)
                                                                                                                        aaaaGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTg                                                                             >  4:860674/1‑94 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  >  4:637294/1‑101 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  <  3:637294/101‑1 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  <  1:913183/101‑1 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  >  2:913183/1‑101 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATCGGGATAGTTtttgtt                                                               <  2:1169637/104‑1 (MQ=255)
                                                                                                                            gCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATCGGGATAGTTtttgtt                                                               >  1:1169637/1‑104 (MQ=255)
                                                                                                                                agaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAcc                                             <  4:151361/118‑1 (MQ=255)
                                                                                                                                agaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAcc                                             >  3:151361/1‑118 (MQ=255)
                                                                                                                                 gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                                                               <  4:37994/83‑1 (MQ=255)
                                                                                                                                 gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAt                                                                               >  3:37994/1‑83 (MQ=255)
                                                                                                                                  atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTggg                                                                           <  4:743936/82‑1 (MQ=255)
                                                                                                                                  atcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTggg                                                                           >  3:743936/5‑86 (MQ=255)
                                                                                                                                      ggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTGTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCt                                >  3:1013968/1‑125 (MQ=255)
                                                                                                                                      ggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCAGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCt                                <  4:1013968/125‑1 (MQ=255)
                                                                                                                                        tCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  <  4:410922/89‑1 (MQ=255)
                                                                                                                                        tCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGttttt                                                                  >  3:410922/1‑89 (MQ=255)
                                                                                                                                            aaaaTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccg  >  4:272089/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCG  >  NC_000913/4377649‑4377937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: